Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,000,244 G→C 8.7% intergenic (+99/+48) uvrD → / ← yigE ssDNA translocase and dsDNA helicase ‑ DNA helicase II/DUF2233 domain‑containing protein YigE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,000,2440GC 8.7% 272.3 / 12.3 89intergenic (+99/+48)uvrD/yigEssDNA translocase and dsDNA helicase ‑ DNA helicase II/DUF2233 domain‑containing protein YigE
Reads supporting (aligned to +/- strand):  ref base G (47/33);  new base C (4/4);  total (51/38)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

CAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGAT  >  NC_000913/4000095‑4000392
                                                                                                                                                     |                                                                                                                                                     
cAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTAcgccg                                                                                                                                                       <  1:13574/150‑1 (MQ=255)
   ggTATTAAATGGGTGGTGGCGGCATTCGCCCGGGTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTGATTTGGCGTAAAAAATCATTCAGTTTCAATCTCTTGCACCGGCATGTAGGCCGGATCGGGCGTTTCcgccgcgt                                                                                                                                                    >  2:1068090/1‑148 (MQ=255)
    gTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATc                                                                                                                                                   >  1:347132/1‑150 (MQ=255)
      aTTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCg                                                                                                                                                 <  2:1003304/150‑1 (MQ=255)
        tAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTTCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGc                                                                                                                                               >  1:1063876/1‑150 (MQ=255)
        tAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGGGTAACGTTGCCGGAGGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATGCGGCGTAAACGCCCTATCCGGc                                                                                                                                               <  2:1477624/150‑1 (MQ=255)
         aaaTGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATGCGGCGTAAACGCCCTATCCGGcc                                                                                                                                              <  1:1265037/150‑2 (MQ=255)
         aaaTGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATGCGGCGTAAACGCCCTATCCGGcc                                                                                                                                              <  1:1270840/150‑2 (MQ=255)
          aaTGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATGCGGCGTAAACGCCCTATCCGGcct                                                                                                                                             <  1:1443745/150‑3 (MQ=255)
          aaTGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCAt                                                                                                                                             >  1:903581/1‑150 (MQ=255)
          aaTGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCAt                                                                                                                                             >  2:1598263/1‑150 (MQ=255)
             ggCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTg                                                                                                                                          >  2:168812/1‑150 (MQ=255)
              gCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTgc                                                                                                                                         >  2:369739/1‑150 (MQ=255)
               cTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTgcg                                                                                                                                        <  1:887329/150‑1 (MQ=255)
                 ggtggCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGcc                                                                                                                                      >  2:668448/1‑150 (MQ=255)
                 ggtggCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGAGTTTAAGCCTCATCAGGCATATgcgcg                                                                                                                                      <  2:741997/150‑2 (MQ=255)
                      cggcATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCtt                                                                                                                                 <  2:273569/150‑1 (MQ=255)
                       ggcATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTc                                                                                                                                <  2:80489/150‑1 (MQ=255)
                         cATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAg                                                                                                                              >  1:606536/1‑150 (MQ=255)
                         cATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAg                                                                                                                              <  1:1133186/150‑1 (MQ=255)
                             cGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATc                                                                                                                          >  2:364323/1‑150 (MQ=255)
                                  ggCTGCAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGaca                                                                                                                     >  2:1239398/1‑150 (MQ=255)
                                     tGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACaa                                                                                                                  <  1:520719/150‑1 (MQ=255)
                                        aGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAAc                                                                                                               >  2:457434/1‑150 (MQ=255)
                                        aGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAAc                                                                                                               >  1:1205773/1‑150 (MQ=255)
                                         gTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACt                                                                                                              >  1:140465/1‑150 (MQ=255)
                                           cGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTAt                                                                                                            >  2:565872/1‑150 (MQ=255)
                                             gtgtAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCg                                                                                                          <  1:882324/150‑1 (MQ=255)
                                                 aaCGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGtttt                                                                                                      <  2:572218/150‑1 (MQ=255)
                                                   cGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTaa                                                                                                    >  2:528424/1‑150 (MQ=255)
                                                     ttGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAAcc                                                                                                  <  1:793023/150‑1 (MQ=255)
                                                       gCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCtt                                                                                                >  1:395696/1‑150 (MQ=255)
                                                            aTGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGc                                                                                           <  1:1195152/150‑1 (MQ=255)
                                                             tGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCt                                                                                          <  1:364323/150‑1 (MQ=255)
                                                             tGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGAATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCt                                                                                          >  2:1111280/1‑150 (MQ=255)
                                                              gCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCCTCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTc                                                                                         >  2:306620/1‑150 (MQ=255)
                                                                     tgcGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGaa                                                                                 <  2:595643/150‑1 (MQ=255)
                                                                     tgcGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGaa                                                                                 >  2:708984/1‑150 (MQ=255)
                                                                     tgcGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGaa                                                                                 >  1:1429442/1‑150 (MQ=255)
                                                                      gcgcACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGaaa                                                                                >  1:902388/1‑150 (MQ=255)
                                                                      gcgcACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGaaa                                                                                >  1:1069272/1‑150 (MQ=255)
                                                                       cgcACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAAt                                                                               >  1:288389/1‑150 (MQ=255)
                                                                        gcACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATc                                                                              <  1:1527091/150‑1 (MQ=255)
                                                                        gcACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATc                                                                              >  2:1310729/1‑150 (MQ=255)
                                                                         cACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGTCGAGCGGaaaggt                                                                             <  1:611792/150‑4 (MQ=255)
                                                                              tATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTg                                                                        <  1:31452/150‑1 (MQ=255)
                                                                              tATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTg                                                                        <  1:26441/150‑1 (MQ=255)
                                                                              tATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTg                                                                        <  1:168812/150‑1 (MQ=255)
                                                                              tATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGAGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTg                                                                        >  1:1387106/1‑150 (MQ=255)
                                                                               aTTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGa                                                                       >  2:716279/1‑150 (MQ=255)
                                                                                tttGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTgtt                                                                      <  1:674294/150‑3 (MQ=255)
                                                                                  tGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACaa                                                                    >  2:90383/1‑150 (MQ=255)
                                                                                  tGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACaa                                                                    >  2:1608106/1‑150 (MQ=255)
                                                                                    gCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACaaaa                                                                  >  2:1400846/1‑150 (MQ=255)
                                                                                     ccTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAAcc                                                                                                  <  1:1411684/118‑1 (MQ=255)
                                                                                     ccTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAAcc                                                                                                  >  2:1411684/1‑118 (MQ=255)
                                                                                     ccTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAg                                                                 <  2:340566/150‑1 (MQ=255)
                                                                                      cTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAgg                                                                >  2:472011/1‑150 (MQ=255)
                                                                                                aTTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGc                                                      >  2:107090/1‑150 (MQ=255)
                                                                                                 ttCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCg                                                                                        <  2:360624/116‑1 (MQ=255)
                                                                                                 ttCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCg                                                                                        >  1:360624/1‑116 (MQ=255)
                                                                                                 ttCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGcc                                                     >  1:785147/1‑150 (MQ=255)
                                                                                                 tacAGATTCAATAAATTAAAAAGTCATATAGCCAGGATATCTCGTTTAACCCCCAAACATAATATAAAACATATTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGcc                                                     >  2:1550009/3‑150 (MQ=255)
                                                                                                  tCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCa                                                    <  1:1167206/150‑1 (MQ=255)
                                                                                                  tCAGATTCAATAAATTGCAAAGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCTACTGAAATCATGGTGACAAAAGGATAGCGTTGCCa                                                    >  2:370027/1‑150 (MQ=255)
                                                                                                    aGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACg                                                  >  1:162783/1‑150 (MQ=255)
                                                                                                     gATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACgg                                                 <  2:288389/150‑1 (MQ=255)
                                                                                                          aaTAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTg                                            <  2:1063876/150‑1 (MQ=255)
                                                                                                          aaTAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCGGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTg                                            >  1:401310/1‑150 (MQ=255)
                                                                                                           aTAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGc                                           >  1:919101/1‑150 (MQ=255)
                                                                                                                  gCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCt                                    >  2:1113480/1‑150 (MQ=255)
                                                                                                                   cAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCtt                                   <  1:1608106/150‑1 (MQ=255)
                                                                                                                     aCGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTc                                 >  1:1232651/1‑150 (MQ=255)
                                                                                                                     aCGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTc                                 <  2:61192/150‑1 (MQ=255)
                                                                                                                      cGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCa                                <  1:457434/150‑1 (MQ=255)
                                                                                                                       gTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCat                               <  2:925964/150‑1 (MQ=255)
                                                                                                                          gcagatgCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCAtata                            >  2:264887/7‑150 (MQ=38)
                                                                                                                          aTGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCAtata                            >  2:712929/1‑150 (MQ=255)
                                                                                                                          aTGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCAtata                            <  1:1598263/150‑1 (MQ=255)
                                                                                                                           tGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATAc                           >  1:183416/1‑150 (MQ=255)
                                                                                                                           tGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATAc                           <  2:903581/150‑1 (MQ=255)
                                                                                                                             gatgCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACat                         >  1:606501/4‑150 (MQ=255)
                                                                                                                                gCCGGATGCGGCGTAAACGCCCTATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCGCCGCCTTTCATATACATAtg                      >  2:1139045/1‑150 (MQ=255)
                                                                                                                                  cGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATAtgtg                    <  1:445879/150‑1 (MQ=255)
                                                                                                                                  cGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATAtgtg                    <  1:1239398/150‑1 (MQ=255)
                                                                                                                                   ggATAGGGCGTTTACGCCGCATCTGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGa                   >  1:739421/1‑150 (MQ=255)
                                                                                                                                        gggCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCg              >  1:807406/1‑150 (MQ=255)
                                                                                                                                           cGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGc           <  2:1499712/150‑1 (MQ=255)
                                                                                                                                               tACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATc       >  2:491897/1‑150 (MQ=255)
                                                                                                                                               tACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATc       >  2:627144/1‑150 (MQ=255)
                                                                                                                                                aCGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCa      >  2:585849/1‑150 (MQ=255)
                                                                                                                                                    cgcATCCGGCATCGGCGCCATCTTCAGTATCTGCCACAAACCTATCGGTTTAACCTTTCCGCGCGACGGAAATCATGGGGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAACTCGGGCCCTCAAGAt  >  1:1210462/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGAT  >  NC_000913/4000095‑4000392

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: