Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,817,708 A→T 5.1% intergenic (+169/‑52) yicC → / → dinD UPF0701 family protein YicC/DNA damage‑inducible protein D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,817,7080AT 5.1% 1038.5 / 28.0 316intergenic (+169/‑52)yicC/dinDUPF0701 family protein YicC/DNA damage‑inducible protein D
Reads supporting (aligned to +/- strand):  ref base A (143/157);  new base T (10/6);  total (153/163)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

Only 100 of 325 total aligned reads displayed.

CTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTA  >  NC_000913/3817559‑3817857
                                                                                                                                                     |                                                                                                                                                     
cTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATaa                                                                                                                                                       >  1:989222/1‑150 (MQ=255)
cTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAAATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATaa                                                                                                                                                       <  2:1991745/150‑1 (MQ=255)
  tttCCGCCTTCACAAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAAt                                                                                                                                                     >  1:340258/1‑150 (MQ=255)
  tttCCGCCTTCACAAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAAt                                                                                                                                                     <  1:1342370/150‑1 (MQ=255)
      cGCCTTCACAAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAg                                                                                                                                                 >  2:1591363/1‑150 (MQ=255)
       gCCTTCACAAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGt                                                                                                                                                <  2:356718/150‑1 (MQ=255)
         cTTCACAAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGGTACTCACAACTGTATTTATATACAGttg                                                                                                                                              <  2:1129931/150‑2 (MQ=255)
         cTTCACAAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTa                                                                                                                                              <  1:1826534/150‑1 (MQ=255)
           tCACAAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACa                                                                                                                                            <  1:1220145/150‑1 (MQ=255)
            cacaAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCGATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATTTATATACAGttgtga                                                                                                                                           <  2:973769/150‑5 (MQ=255)
            cacaAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAg                                                                                                                                           <  2:1444140/150‑1 (MQ=255)
            cacaAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAg                                                                                                                                           >  1:697088/1‑150 (MQ=255)
             acaAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATTTATATACAGttgtgag                                                                                                                                          <  1:1814878/150‑6 (MQ=255)
               aaaTGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGAtt                                                                                                                                        >  2:53316/1‑150 (MQ=255)
                   gccgccACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTAct                                                                                                                                    <  2:961946/150‑1 (MQ=255)
                    ccgccACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTActt                                                                                                                                   >  2:1373267/1‑150 (MQ=255)
                    ccgccACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTActt                                                                                                                                   >  1:857041/1‑150 (MQ=255)
                     cgccACTCAAACAGAGCGGCATTTTTCTTCCCCGCGACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTActtt                                                                                                                                  >  1:48302/1‑150 (MQ=255)
                        cACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTActttctt                                                                                                                               >  2:920675/1‑150 (MQ=255)
                         aCTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTActttcttt                                                                                                                              <  1:195607/150‑1 (MQ=255)
                          cTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTg                                                                                                                             <  2:448984/150‑1 (MQ=255)
                          cTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTg                                                                                                                             <  1:411472/150‑1 (MQ=255)
                          cTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTg                                                                                                                             <  2:149555/150‑1 (MQ=255)
                           tCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGc                                                                                                                            <  2:741344/150‑1 (MQ=255)
                           tCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGc                                                                                                                            <  1:1944742/150‑1 (MQ=255)
                            cAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCa                                                                                                                           >  2:678242/1‑150 (MQ=255)
                            cAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCa                                                                                                                           >  2:405150/1‑150 (MQ=255)
                             aaaCAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCaa                                                                                                                          >  2:1051097/1‑150 (MQ=255)
                                cAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTg                                                                                                                       <  1:191949/150‑1 (MQ=255)
                                 agagCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGa                                                                                                                      <  1:1967233/150‑1 (MQ=255)
                                 agagCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGa                                                                                                                      <  1:530838/150‑1 (MQ=255)
                                      ggCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATca                                                                                                                 >  2:219321/1‑150 (MQ=255)
                                        cATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTActttcttt                                                                                                                              <  2:1236688/135‑1 (MQ=255)
                                         aTTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACAt                                                                                                              <  2:1564370/150‑1 (MQ=255)
                                               cttcCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTg                                                                                                        >  2:666864/1‑150 (MQ=255)
                                                tttcccGCAACATTTAATTCTTTTTTTGGTGGCTGCTCCATTTTTTGCGTTGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTgg                                                                                                       <  2:767378/147‑1 (MQ=255)
                                                ttcCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATaa                                                                                                                                                       >  1:1132908/1‑102 (MQ=255)
                                                 tcCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCGCAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTggg                                                                                                      >  2:1926064/1‑150 (MQ=255)
                                                    ccGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCaa                                                                                                   >  2:1634759/1‑150 (MQ=255)
                                                     cGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAAt                                                                                                  <  2:1667578/150‑1 (MQ=255)
                                                      gCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATg                                                                                                 <  2:1318319/150‑1 (MQ=255)
                                                        aaCATTCAATTATGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCGCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATgaa                                                                                               <  2:610363/150‑1 (MQ=255)
                                                         aCATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATgaac                                                                                              <  1:1404453/150‑1 (MQ=255)
                                                            ttCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATgaacgaa                                                                                           <  1:670753/150‑1 (MQ=255)
                                                              cAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAAca                                                                                         <  2:1182475/150‑1 (MQ=255)
                                                               aaTTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAAcat                                                                                        >  1:1390360/1‑150 (MQ=255)
                                                                  tCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAAcatcat                                                                                     >  2:329799/1‑150 (MQ=255)
                                                                    tGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACAtcatca                                                                                   >  1:212220/1‑150 (MQ=255)
                                                                        ttGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCt                                                                               <  1:633492/150‑1 (MQ=255)
                                                                          gCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCttt                                                                             >  1:556936/1‑150 (MQ=255)
                                                                           cGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCtttt                                                                            >  1:1872762/1‑150 (MQ=255)
                                                                            gTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTg                                                                           <  1:1371268/150‑1 (MQ=255)
                                                                              gCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGaa                                                                         <  1:1211863/150‑1 (MQ=255)
                                                                               ccTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAag                                                                        >  1:1549342/1‑150 (MQ=255)
                                                                                  gCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAgaga                                                                     <  2:1167424/150‑1 (MQ=255)
                                                                                  gCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAgaga                                                                     <  2:697088/150‑1 (MQ=255)
                                                                                   cTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGAt                                                                    >  1:1279241/1‑150 (MQ=255)
                                                                                   cTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGAt                                                                    <  1:1634759/150‑1 (MQ=255)
                                                                                        gATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAAc                                                               >  2:1898383/1‑150 (MQ=255)
                                                                                        gATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAAc                                                               >  2:1765588/1‑150 (MQ=255)
                                                                                         aTTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACt                                                              <  1:400681/150‑1 (MQ=255)
                                                                                          ttttGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTg                                                             <  1:894616/150‑1 (MQ=255)
                                                                                             tGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGAtt                                                          <  2:1351049/150‑1 (MQ=255)
                                                                                                gATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAAt                                                       <  1:936354/150‑1 (MQ=255)
                                                                                                 aTGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATg                                                      <  2:1398476/150‑1 (MQ=255)
                                                                                                  tGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTTCAATTTATATCACATTGAGTGGGAAATGAAAGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGc                                                     <  2:24580/150‑1 (MQ=255)
                                                                                                  tGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGc                                                     >  1:1566842/1‑150 (MQ=255)
                                                                                                  tGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGc                                                     >  1:27754/1‑150 (MQ=255)
                                                                                                  tGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGc                                                     <  2:1872762/150‑1 (MQ=255)
                                                                                                      tttttGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAAc                                                 >  2:566588/1‑150 (MQ=255)
                                                                                                        tttGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACgg                                               <  1:1183469/150‑1 (MQ=255)
                                                                                                            cccAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCa                                           <  1:1672812/150‑1 (MQ=255)
                                                                                                             ccAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAg                                          <  1:1898383/150‑1 (MQ=255)
                                                                                                                gCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAAc                                       <  1:187677/150‑1 (MQ=255)
                                                                                                                       gagcACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTc                                <  2:857041/150‑1 (MQ=255)
                                                                                                                       gagaACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTc                                <  2:141489/150‑1 (MQ=255)
                                                                                                                        agaACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCa                               >  2:80827/1‑150 (MQ=255)
                                                                                                                           aCGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCa                            <  1:1540134/150‑1 (MQ=255)
                                                                                                                               gagaTACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGCACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGAc                        >  2:82943/1‑150 (MQ=255)
                                                                                                                                    aCTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACtt                   >  1:832454/1‑150 (MQ=255)
                                                                                                                                    aCTCACAACTGTATATAAATACAGTTACAGATTTAATTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACtt                   >  2:1802839/1‑150 (MQ=255)
                                                                                                                                     cTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTg                  >  1:1239658/1‑150 (MQ=255)
                                                                                                                                      ctgtaaCTGTATTTATATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACttcg                 >  1:238377/5‑148 (MQ=255)
                                                                                                                                      ctgtaaCTGTATTTATATACAGTTACAGATTTACTTTCTATGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTgg                 >  2:1959099/5‑150 (MQ=255)
                                                                                                                                       cacaACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCa                            >  1:1018217/1‑138 (MQ=255)
                                                                                                                                       cacaACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCa                            <  2:1018217/138‑1 (MQ=255)
                                                                                                                                       cacaACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAAc                     >  1:751794/1‑145 (MQ=255)
                                                                                                                                       cacaACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTggg                >  2:28904/1‑150 (MQ=255)
                                                                                                                                       cacaACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTggg                >  1:1475621/1‑150 (MQ=255)
                                                                                                                                       cacaACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTggg                >  2:386399/1‑150 (MQ=255)
                                                                                                                                       cacaACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTggg                >  2:1368842/1‑150 (MQ=255)
                                                                                                                                        gtaaCTGTATTTATATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGa               >  1:432292/3‑150 (MQ=255)
                                                                                                                                           aCTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAac            >  2:498189/1‑150 (MQ=255)
                                                                                                                                            cTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAact           <  1:1103712/150‑1 (MQ=255)
                                                                                                                                            cTGAATAAAAATACAGTTAAAGATTTAATTTCTTTGAAATTCATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAATAGAAAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAact           >  2:1910629/1‑150 (MQ=255)
                                                                                                                                                atataAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTg       >  1:286352/1‑150 (MQ=255)
                                                                                                                                                 tataAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTgg      >  1:1711276/1‑150 (MQ=255)
                                                                                                                                                   taAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGt    >  1:59505/1‑150 (MQ=255)
                                                                                                                                                   taAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGt    >  2:1435646/1‑150 (MQ=255)
                                                                                                                                                     aaTACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTa  >  2:121933/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTA  >  NC_000913/3817559‑3817857

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