Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,076,294 A→G 5.2% K506E (AAA→GAA)  yeeR → CP4‑44 prophage; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,076,2940AG 5.2% 239.5 / 2.9 77K506E (AAA→GAA) yeeRCP4‑44 prophage; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/43);  new base G (2/2);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

TGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTG  >  NC_000913/2076149‑2076439
                                                                                                                                                 |                                                                                                                                                 
tGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGGTGAAACCATGAGGAATAAGCTGCGCGAAGGTAAAg                                                                                                                                                >  2:63515/1‑149 (MQ=255)
  cAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATTATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAg                                                                                                                                                >  5:390385/1‑147 (MQ=255)
    gCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTaga                                                                                                                                            <  8:23000/149‑1 (MQ=255)
    gCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTaga                                                                                                                                            <  2:174464/149‑1 (MQ=255)
    gCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTCACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTa                                                                                                                                              >  8:45136/1‑147 (MQ=255)
        gTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTGAAGTagag                                                                                                                                           >  1:369240/1‑146 (MQ=255)
         tAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGTCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGTCGCTGTTGAAACAATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGgaga                                                                                                                                       <  7:133799/149‑1 (MQ=255)
          aGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTGAAGTAgagga                                                                                                                                         <  2:369240/146‑1 (MQ=255)
              tCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGa                                                                                                                                  <  1:234662/149‑1 (MQ=255)
              tCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGa                                                                                                                                  >  4:186353/1‑149 (MQ=255)
                         ttCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAg                                                                                                                       >  1:147530/1‑149 (MQ=255)
                         ttCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAg                                                                                                                       >  6:368664/1‑149 (MQ=255)
                         ttCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAg                                                                                                                       >  4:232739/1‑149 (MQ=255)
                            gACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAgcg                                                                                                                     >  3:251792/1‑148 (MQ=255)
                             aCGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCa                                                                                                                   <  7:204034/149‑1 (MQ=255)
                              cGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAg                                                                                                                  >  8:205928/1‑149 (MQ=255)
                               ggACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGc                                                                                                                 >  3:362866/1‑149 (MQ=255)
                               ggACAAAAAGTACAGTATCGTTACATGGGAAATAATGGCACCATGCAACTGGAATTCCCCCTTGAGCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAATTAGAGGAGAGATGACGCCCGGCGAGCACAGc                                                                                                                 <  2:107907/149‑1 (MQ=255)
                                 aCAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAgcg                                                                                                                     >  8:247590/1‑143 (MQ=255)
                                 aCAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAgcg                                                                                                                     <  7:247590/143‑1 (MQ=255)
                                   aaaaCGGACAGTATCGTTACATCGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGAGCAGTATCCCGGCGCTGGTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCa                                                                                                             <  3:36427/149‑1 (MQ=255)
                                    aaaCGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGATGAGAGATGACGCCCGGCGAGCGCAGCCTCAt                                                                                                            <  5:96732/149‑1 (MQ=255)
                                     aaCGGACAGTATCGTTACATGGGAAATAATTGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCAtt                                                                                                           <  2:117779/149‑1 (MQ=255)
                                     aaCGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCAtt                                                                                                           <  8:243087/149‑1 (MQ=255)
                                      aCGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTc                                                                                                          <  5:354272/149‑1 (MQ=255)
                                         gACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAAc                                                                                                       <  5:368664/149‑1 (MQ=255)
                                             gTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGCGAGGACGCCCGGCGAGCGCAGCCTCATTCAACGGGc                                                                                                   >  6:31479/1‑149 (MQ=255)
                                              tATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGCCGCCCGGCGAGCGCAGCCTCATTCAACGGGcc                                                                                                  >  5:294320/1‑149 (MQ=255)
                                              tATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGcc                                                                                                  >  4:131851/1‑149 (MQ=255)
                                                tCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCACGGCGAGCGCAGCCTCATTCAACGGGCCCt                                                                                                >  1:54541/1‑149 (MQ=255)
                                                   ttACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCAGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGaa                                                                                             <  5:118851/149‑1 (MQ=255)
                                                        tGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAAccc                                                                                        <  5:36054/149‑1 (MQ=255)
                                                         gggAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCt                                                                                       >  7:195054/1‑149 (MQ=255)
                                                          ggAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTg                                                                                      <  3:186353/149‑1 (MQ=255)
                                                           gAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTg                                                                                      >  8:4217/1‑148 (MQ=255)
                                                            aaaTAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGa                                                                                    <  6:144911/149‑1 (MQ=255)
                                                              ataatGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGAcc                                                                                  >  1:55305/1‑149 (MQ=255)
                                                                 atGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGcc                                                                               <  4:51076/149‑1 (MQ=255)
                                                                  tGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCa                                                                              >  1:62800/1‑149 (MQ=255)
                                                                    gCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATc                                                                            <  6:225935/149‑1 (MQ=255)
                                                                      cccATGCAACTGGAAGTCCACCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCAGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAAAGGGCCCTGAAAACCCTGGCCCGCCCTCtt                                                                          >  7:345200/1‑149 (MQ=255)
                                                                       ccATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTGAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCtt                                                                          >  1:378530/1‑148 (MQ=255)
                                                                        cATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCa                                                                        <  6:375939/149‑1 (MQ=255)
                                                                        cATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCa                                                                        >  3:7471/1‑149 (MQ=255)
                                                                         aTGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCAt                                                                       <  3:232739/149‑1 (MQ=255)
                                                                           gCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAACATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCCTTCAACGGGCCCTGAACACCCTGGCCCGCCATCTTCATg                                                                      >  2:26134/1‑148 (MQ=255)
                                                                                 ggAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTagag                                                                                                                                           <  2:118516/73‑1 (MQ=255)
                                                                                     gTCCCCCGTGATCAGTAGGCCGGCGCTGTTGAAACCATGAGGAATACGAGCCGCGAAGGTAATGGAGAGGAGAGATGACGCCCGGCGCGCGCCGCCTCATTCGACGGGCCCTGAACCCCCGGGCCCGCCCTCTTGGTGAACCCGGCggg                                                           >  6:316844/1‑147 (MQ=255)
                                                                                      tCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTgg                                                          <  2:147530/149‑1 (MQ=255)
                                                                                        ccccGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGcc                                                        <  4:251792/149‑1 (MQ=255)
                                                                                        ccccGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGc                                                         <  4:309408/148‑1 (MQ=255)
                                                                                        ccccGTGAGCAGTAGGCCGGCGCTGTTGAAACCATGAGGAATAAGAGCCGCGAACGTAAAGTCGAGGATAGATGACGCCCGGAGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGcc                                                        <  4:146/149‑1 (MQ=255)
                                                                                          ccGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCtt                                                      <  6:390385/149‑1 (MQ=255)
                                                                                           cGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTc                                                     <  7:202376/149‑1 (MQ=255)
                                                                                            gTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCa                                                    <  4:362866/149‑1 (MQ=255)
                                                                                              gATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCACGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCAcc                                                  >  7:232935/1‑149 (MQ=255)
                                                                                               aTCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCt                                                 <  8:195054/149‑1 (MQ=255)
                                                                                                tCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTGAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTc                                                <  2:378530/149‑1 (MQ=255)
                                                                                                  aGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCa                                              <  7:402078/149‑1 (MQ=255)
                                                                                                        ccGGCGCTGTTTAAACCTTGAGTAATATGCTCCGCGAAGGTAAAGTAGATGAGAGATGACGCCCGGCGAGCGCAGACTCATTCAACGGGCCCTGAAAACCCTGGACCTCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGt                                         <  8:345200/148‑1 (MQ=255)
                                                                                                        ccGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTg                                        >  4:183171/1‑149 (MQ=255)
                                                                                                        ccGGCGCTGTGGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCACATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGt                                         <  4:83338/148‑1 (MQ=255)
                                                                                                          ggCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTgca                                      <  4:7471/149‑1 (MQ=255)
                                                                                                                 ttGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCAtt                                                                                                           >  7:344170/1‑73 (MQ=255)
                                                                                                                 ttGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCAtt                                                                                                           <  8:344170/73‑1 (MQ=255)
                                                                                                                   gAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGaa                             >  7:278857/1‑149 (MQ=255)
                                                                                                                    aaaCCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTc                                                                                                              <  5:373430/67‑1 (MQ=255)
                                                                                                                    aaaCCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTc                                                                                                              >  6:373430/1‑67 (MQ=255)
                                                                                                                     aaCCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATg                           <  7:45136/149‑1 (MQ=255)
                                                                                                                       ccATGAGGAATAAGATCCGCGAAGGTAACGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGc                         <  5:31479/149‑1 (MQ=255)
                                                                                                                       ccATGAGGAAAAAGCTCCGCGAACGGAAAGTAGAGGAGAGATGACGCCCGGCCAGCGCACCCTCCTTCAACCGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGGGGCCTTCACCTCCACCCGAGCAGCACGGGAATGGc                         <  2:322500/149‑1 (MQ=255)
                                                                                                                        cATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCt                        <  2:55305/149‑1 (MQ=255)
                                                                                                                          tGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAAACCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGa                      >  4:267859/1‑149 (MQ=255)
                                                                                                                            aGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGAtt                    <  7:205928/149‑1 (MQ=255)
                                                                                                                               aaTAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTg                 <  3:131851/149‑1 (MQ=255)
                                                                                                                                           gAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGa     <  7:4217/149‑1 (MQ=255)
                                                                                                                                               gTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTg  >  4:216253/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTG  >  NC_000913/2076149‑2076439

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑