Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,065,747 (C)6→5 coding (18/546 nt) cobU ← cobinamide‑P guanylyltransferase/cobinamide kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,065,7420C.99.5% 2385.7 / 12.2 571coding (23/546 nt)cobUcobinamide‑P guanylyltransferase/cobinamide kinase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (252/311);  total (256/315)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 606 total aligned reads displayed.

GTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAA  >  NC_000913/2065603‑2065877
                                                                                                                                            |                                                                                                                                       
   cTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGcccc                                                                                                                                      >  2:657503/1‑140 (MQ=255)
    tCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGccccc                                                                                                                                     <  1:3610123/140‑1 (MQ=255)
     ccGGGCGGCCTTGCCGATGATGCTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGcccccg                                                                                                                                    <  1:3266956/140‑2 (MQ=255)
             ccTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCaaa                                                                                                                           >  2:3786885/1‑140 (MQ=255)
             ccTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCaaa                                                                                                                           >  2:6237698/1‑140 (MQ=255)
             ccTTGCCGATGATGTTCAATCCGTGCAGCCACCTCATCACCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCaaa                                                                                                                           <  2:2618977/140‑1 (MQ=255)
              cTTGCCGATGATGTTCTATCCGCGCAGCCTTCTCATCCTCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCaaaa                                                                                                                          <  1:2592847/140‑1 (MQ=255)
                   cGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAatcatc                                                                                                                     <  1:3812039/140‑1 (MQ=255)
                   cGATGATGTTCTATCCGCGCAGCCATCTCATCATCCAGGGTTTGCGAGCTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAatcatc                                                                                                                     <  2:2085962/140‑1 (MQ=255)
                       gatgTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATAc                                                                                                                 >  2:5075544/1‑140 (MQ=255)
                       gatgTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATAc                                                                                                                 >  2:5909871/1‑140 (MQ=255)
                        atgTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATAcc                                                                                                                <  2:895688/140‑1 (MQ=255)
                         tgTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCg                                                                                                               <  2:4885222/140‑1 (MQ=255)
                         tgTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCg                                                                                                               <  1:4067530/140‑1 (MQ=255)
                         tgTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCg                                                                                                               >  1:4609152/1‑140 (MQ=255)
                           tttattCCGTGCAGCCATCTAATCATCAAGAATTTGCGAGTTGGCGATATAAAAACCCTGGAAGAAGCCCCCAATAAGCGCCCCTGCGTGGCGCTTCTCCCCGCTCCGGCCGCCCCGTCACCAAAATCATCATACCGACt                                                                                                            <  2:7362397/135‑1 (MQ=255)
                              tATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTcc                                                                                                          <  1:5410816/140‑1 (MQ=255)
                              tATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTcc                                                                                                          <  2:321569/140‑1 (MQ=255)
                                tCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTg                                                                                                        >  2:6769423/1‑140 (MQ=255)
                                tCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTg                                                                                                        >  1:5325660/1‑140 (MQ=255)
                                  cGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGAt                                                                                                      >  2:6298025/1‑140 (MQ=255)
                                    tGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTc                                                                                                    <  2:1662897/140‑1 (MQ=255)
                                    tGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTc                                                                                                    >  1:6178237/1‑140 (MQ=255)
                                    tGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTc                                                                                                    >  1:5546491/1‑140 (MQ=255)
                                    tGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGATCCGTGCTCCCCCGTGACCGAAGTCATCAGACCGACTCCTGATTc                                                                                                    >  1:2313098/1‑140 (MQ=255)
                                     gCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCa                                                                                                   >  1:5203521/1‑140 (MQ=255)
                                     gCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCa                                                                                                   >  1:6902446/1‑140 (MQ=255)
                                      cAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTcc                                                                                                          >  1:4909296/1‑132 (MQ=255)
                                      cAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCaa                                                                                                  <  2:4015562/140‑1 (MQ=255)
                                      cAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCaa                                                                                                  <  1:5855689/140‑1 (MQ=255)
                                      cAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGA‑‑‑‑‑‑‑‑‑CTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTatgat                                                                                         >  1:6703436/1‑140 (MQ=255)
                                        gCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAAt                                                                                                >  2:6682737/1‑140 (MQ=255)
                                         ccATCTCATCATCAAGGATTTGAGAGGTAGCCATACACATAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCACTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATg                                                                                               <  2:2267688/140‑1 (MQ=255)
                                           aTCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTa                                                                                             <  2:1328916/140‑1 (MQ=255)
                                           aTCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTa                                                                                             <  1:7573047/140‑1 (MQ=255)
                                           aTCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTa                                                                                             <  1:7414181/140‑1 (MQ=255)
                                           aTCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTa                                                                                             <  1:224313/140‑1 (MQ=255)
                                             ctcATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTatg                                                                                           <  1:1273888/140‑1 (MQ=255)
                                               catcatcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTatgat                                                                                         <  1:1755232/140‑1 (MQ=255)
                                               catcatcaAGGATTTGAGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTatgat                                                                                         <  2:4606793/140‑1 (MQ=255)
                                                  catcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTc                                                                                                    >  2:6452868/1‑126 (MQ=255)
                                                  catcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGt                                                                                      <  1:999107/140‑1 (MQ=255)
                                                     caAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCATCACCAAAATCAACATGTCGACTCCGGATTAAAATGTATGATGGTAtt                                                                                   >  2:3662858/1‑140 (MQ=255)
                                                       aGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCa                                                                                 <  1:5024589/140‑1 (MQ=255)
                                                       aGGATTTGCGAGGTAGCAATATACAGAACCTATGAAGAGACCCCAATATGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCa                                                                                 <  2:2423890/140‑1 (MQ=255)
                                                       aGGATTTGCGAGGTACCGATATACAAAACCTGTGAAGCGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCa                                                                                 <  1:3662858/140‑1 (MQ=255)
                                                            ttGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAATTCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaa                                                                            <  1:6456585/140‑1 (MQ=255)
                                                            ttGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaa                                                                            <  1:1069377/140‑1 (MQ=255)
                                                            ttGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaa                                                                            <  2:5310543/140‑1 (MQ=255)
                                                            ttGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaa                                                                            <  2:4802612/140‑1 (MQ=255)
                                                            ttGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaa                                                                            <  1:4245642/140‑1 (MQ=255)
                                                             tGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaa                                                                            >  1:1629890/1‑139 (MQ=255)
                                                                 aGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaacagaa                                                                       <  1:6298025/140‑1 (MQ=255)
                                                                   gTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaacagaaca                                                                     >  2:6039097/1‑140 (MQ=255)
                                                                     aGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAAt                                                                   <  2:5009027/140‑1 (MQ=255)
                                                                       cGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACTAAAATCATCATACCGACTCCTAATTTAAATGTATGATGATATTCATTGAACAGAACAATgg                                                                 >  2:2329177/1‑140 (MQ=255)
                                                                            taCAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACGCCTGATTCCAATCTATTATGGTATTAATTGTACAGAACAATGGATccc                                                            >  2:2978310/1‑140 (MQ=255)
                                                                               aGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCa                                                         <  1:657503/140‑1 (MQ=255)
                                                                                 aaCCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCAcc                                                       >  2:4469461/1‑140 (MQ=255)
                                                                                 aaCCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCAcc                                                       >  2:1535688/1‑140 (MQ=255)
                                                                                 aaCCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCAcc                                                       >  1:7198633/1‑140 (MQ=255)
                                                                                 aaCCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCAGTGCACAGAAGAATGGATCCCCCAcc                                                       >  1:3431690/1‑140 (MQ=255)
                                                                                   ccTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGt                                                     >  2:779450/1‑140 (MQ=255)
                                                                                    cTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGtt                                                    >  2:7031355/1‑140 (MQ=255)
                                                                                      gtgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCACATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGc                                                  >  2:3317007/1‑140 (MQ=255)
                                                                                      gtgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGa                                                                <  1:6352686/126‑1 (MQ=255)
                                                                                      gtgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGa                                                                >  2:6352686/1‑126 (MQ=255)
                                                                                      gtgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGc                                                  >  2:5171156/1‑140 (MQ=255)
                                                                                      gtgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGc                                                  >  1:7102003/1‑140 (MQ=255)
                                                                                      gtgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGc                                                  >  2:4220285/1‑140 (MQ=255)
                                                                                      gtgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATCGAAGAGTACAATGGATCCCCCACCGTTGc                                                  >  1:2409303/1‑140 (MQ=255)
                                                                                      gtgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATCGAACAGAACAATGGATCCCCCACCGTTGc                                                  >  1:3358889/1‑140 (MQ=255)
                                                                                       tgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCt                                                 >  1:878803/1‑140 (MQ=255)
                                                                                       tgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCt                                                 >  1:1823976/1‑140 (MQ=255)
                                                                                       tgAAGAGTCCCCAATAAACGCCCCTGCGTGGTGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCt                                                 <  2:2721784/140‑1 (MQ=255)
                                                                                       tgAAGAGCCCCCAATCAGCGCCCCTCCGTCGCGACTCGTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCt                                                 <  1:2055524/140‑1 (MQ=255)
                                                                                               tCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCGTTGGACAGAACAATGGGTCCCTCACCGTTGCTGACCGCa                                          >  2:1992249/1‑140 (MQ=255)
                                                                                               tCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCGTTGCACAAAACAATGGTTCCCCAACCGTTGCTGACCGCa                                          >  2:3481242/1‑140 (MQ=255)
                                                                                               tCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCGTTGAACAGAACAATGGTTCCCCCACCGTTGCTGACCGCa                                          >  2:2837225/1‑140 (MQ=255)
                                                                                               tCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCa                                          >  1:5687097/1‑140 (MQ=255)
                                                                                               tCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAACGGATCCCCCCCCGTTGCTGACCGCa                                          >  1:2721784/1‑140 (MQ=255)
                                                                                                ccccAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCaa                                         >  2:1763260/1‑140 (MQ=255)
                                                                                                  ccAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCaaaa                                       >  1:7025852/1‑140 (MQ=255)
                                                                                                         gcgcCTCTGCGTGGCGGCTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAc                                <  1:4460390/140‑1 (MQ=255)
                                                                                                         gcgcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAc                                <  2:5396010/140‑1 (MQ=255)
                                                                                                         gcgcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAc                                >  1:392097/1‑140 (MQ=255)
                                                                                                         gcgcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAc                                <  1:5111320/140‑1 (MQ=255)
                                                                                                         gcgcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAc                                <  1:6769423/140‑1 (MQ=255)
                                                                                                         gcgcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAc                                <  1:7347959/140‑1 (MQ=255)
                                                                                                         gcgcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAc                                <  2:5687097/140‑1 (MQ=255)
                                                                                                         gcgcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAc                                <  2:5546491/140‑1 (MQ=255)
                                                                                                         gcgcCTCTGCGTGGCCACTCTTCCCACTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAc                                <  1:2817557/140‑1 (MQ=255)
                                                                                                           gcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAcgc                              <  1:481144/140‑1 (MQ=255)
                                                                                                           gcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAcgc                              <  2:7469075/140‑1 (MQ=255)
                                                                                                                     ggCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACt                    <  2:7198633/140‑1 (MQ=255)
                                                                                                                       cGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTAc                                <  1:6002825/126‑1 (MQ=255)
                                                                                                                          ctctTCCCGCTCCGTGACCCCCTGTCACCAATATCCTCATACCGACTCTTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGcgcg               <  2:3403771/140‑1 (MQ=255)
                                                                                                                               cccGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAata          <  1:2912965/140‑1 (MQ=255)
                                                                                                                                   cTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAAt      >  2:7104793/1‑140 (MQ=255)
                                                                                                                                       gTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTaaa  >  1:1381819/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                       
GTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAA  >  NC_000913/2065603‑2065877

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: