Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,099,661 (A)5→6 intergenic (‑19/+164) yggN ← / ← b2959 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,099,6561.A99.6% 1044.6 / ‑4.2 275intergenic (‑14/+169)yggN/b2959hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (134/140);  total (135/140)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 305 total aligned reads displayed.

TGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGC  >  NC_000913/3099564‑3099743
                                                                                             |                                                                                       
 gACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCCCTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGa                                                                                        <  2:1658277/93‑1 (MQ=255)
 gACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGa                                                                                        <  1:5814682/93‑1 (MQ=255)
       gCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaat                                                                                  <  2:1194555/93‑3 (MQ=255)
        cACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaatc                                                                                 <  1:2223024/93‑4 (MQ=255)
        cACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaatc                                                                                 <  1:4604221/93‑4 (MQ=255)
         aCTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaatca                                                                                >  2:6366438/1‑89 (MQ=255)
         aCTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaatca                                                                                <  2:5182315/93‑5 (MQ=255)
         aCTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaatca                                                                                >  1:2110009/1‑89 (MQ=255)
         aCTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaatca                                                                                >  1:5982233/1‑89 (MQ=255)
           tGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAgg                                                                               >  1:4645034/1‑93 (MQ=255)
            ggTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGc                                                                              <  1:3037949/93‑1 (MQ=255)
             gTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGcc                                                                             <  1:4091697/93‑1 (MQ=255)
             gTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGcc                                                                             <  2:5722339/93‑1 (MQ=255)
              tAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCa                                                                            <  1:2403190/93‑1 (MQ=255)
               aGTCGGCGTGAGCGGTCATTGCCGTTACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAg                                                                           >  1:72211/1‑93 (MQ=255)
                 tCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCa                                                                         >  1:4374490/1‑93 (MQ=255)
                  cGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAt                                                                        <  1:4522964/93‑1 (MQ=255)
                   ggCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAtt                                                                       <  1:2919727/93‑1 (MQ=255)
                   ggCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTCTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAtt                                                                       <  2:5132770/93‑1 (MQ=255)
                    gCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAttt                                                                      >  2:6207752/1‑93 (MQ=255)
                      gTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGc                                                                    >  1:3510515/1‑93 (MQ=255)
                       tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                                   <  1:4809594/93‑1 (MQ=255)
                       tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                                   >  1:3673898/1‑93 (MQ=255)
                        gAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTg                                                                  >  2:6066539/1‑93 (MQ=255)
                          gCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGc                                                                >  2:3140767/1‑93 (MQ=255)
                            ggTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGccc                                                              <  1:786298/93‑1 (MQ=255)
                            ggTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGccc                                                              >  1:4431470/1‑93 (MQ=255)
                             gTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGcccc                                                             >  2:842784/1‑93 (MQ=255)
                                 ttGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAt                                                         >  1:999492/1‑93 (MQ=255)
                                 ttGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAt                                                         >  2:3197158/1‑93 (MQ=255)
                                     cGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGAt                                                     >  2:3758687/1‑93 (MQ=255)
                                     cGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGAt                                                     <  2:2813973/93‑1 (MQ=255)
                                      gTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAATAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATa                                                    <  1:2476229/93‑1 (MQ=255)
                                         aCTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACa                                                 >  1:6640340/1‑93 (MQ=255)
                                         aCTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACa                                                 >  2:6118541/1‑93 (MQ=255)
                                          cTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaa                                                >  2:4861859/1‑93 (MQ=255)
                                          cTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaa                                                >  2:3418058/1‑93 (MQ=255)
                                           tGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa                                               <  2:4682476/93‑1 (MQ=255)
                                            gAAAGTGCTGCCGCCAGCAGCGTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                                              >  1:3132398/1‑93 (MQ=255)
                                            gAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                                              <  1:1210329/93‑1 (MQ=255)
                                             aaaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGt                                             <  2:753854/93‑1 (MQ=255)
                                              aaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa                                            <  2:1877386/93‑1 (MQ=255)
                                              aaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa                                            <  2:5551087/93‑1 (MQ=255)
                                              aaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa                                            <  2:6298065/93‑1 (MQ=255)
                                              aaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa                                            <  1:3158926/93‑1 (MQ=255)
                                               aGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAc                                           <  2:594587/93‑1 (MQ=255)
                                               aGTGCTGCCGCCAGCAGCAGTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAc                                           <  1:904407/93‑1 (MQ=255)
                                                  gctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgcg                                        <  2:2546309/93‑1 (MQ=255)
                                                   ctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGt                                       <  2:1441716/93‑1 (MQ=255)
                                                   ctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGt                                       <  2:898068/93‑1 (MQ=255)
                                                   ctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGt                                       <  1:6411221/93‑1 (MQ=255)
                                                    tgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGtt                                      <  2:4311802/93‑1 (MQ=255)
                                                    tgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGtt                                      <  1:4750009/93‑1 (MQ=255)
                                                     gccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGTCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGttt                                     <  1:1777918/93‑1 (MQ=255)
                                                      ccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGtttt                                    <  2:302596/93‑1 (MQ=255)
                                                        gccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTcc                                  <  2:4221649/93‑1 (MQ=255)
                                                        gccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTcc                                  >  2:3627664/1‑93 (MQ=255)
                                                        gccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTcc                                  >  1:5128157/1‑93 (MQ=255)
                                                         ccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCa                                 >  1:2305895/1‑93 (MQ=255)
                                                          cagcagCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCAt                                >  2:4111984/1‑93 (MQ=255)
                                                            gcagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGc                              >  1:2557782/1‑93 (MQ=255)
                                                            gcagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGc                              >  1:2703180/1‑93 (MQ=255)
                                                            gcagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGc                              >  1:4556975/1‑93 (MQ=255)
                                                            gcagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTCTTCCATGc                              >  2:6599137/1‑93 (MQ=255)
                                                             cagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGcc                             <  2:5437651/93‑1 (MQ=255)
                                                               gcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCgg                           <  1:1851319/93‑1 (MQ=255)
                                                                caTTTTGCGTATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGa                          <  2:1516359/93‑1 (MQ=255)
                                                                 aTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGAt                         >  2:2916680/1‑93 (MQ=255)
                                                                  ttttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATg                        >  1:4154073/1‑93 (MQ=255)
                                                                   tttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgc                       <  1:554142/93‑1 (MQ=255)
                                                                      gcgcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcggc                    >  1:111182/1‑93 (MQ=255)
                                                                       cgcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcggcg                   >  2:5279812/1‑93 (MQ=255)
                                                                         catcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTa                 <  2:4475928/93‑1 (MQ=255)
                                                                          atcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTaa                >  1:6035719/1‑93 (MQ=255)
                                                                          atcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTaa                >  1:3859400/1‑93 (MQ=255)
                                                                           tcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTaaa               >  1:6496314/1‑93 (MQ=255)
                                                                            catAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAac              >  1:2685956/1‑93 (MQ=255)
                                                                             atAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAaca             >  1:2782180/1‑93 (MQ=255)
                                                                             atAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAaca             >  1:1478666/1‑93 (MQ=255)
                                                                              tAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAacac            >  2:4234991/1‑93 (MQ=255)
                                                                               aGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACAcc           <  1:1011902/93‑1 (MQ=255)
                                                                               aGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACAcc           >  2:229237/1‑93 (MQ=255)
                                                                               aGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACAcc           >  1:4218118/1‑93 (MQ=255)
                                                                               aGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACAcc           <  2:3345055/93‑1 (MQ=255)
                                                                               aGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACAcc           <  2:260060/93‑1 (MQ=255)
                                                                                gTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCt          >  2:5481865/1‑93 (MQ=255)
                                                                                gTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCt          <  2:2423107/93‑1 (MQ=255)
                                                                                 tCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCtt         >  1:3779927/1‑93 (MQ=255)
                                                                                 tCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCtt         <  2:3882256/93‑1 (MQ=255)
                                                                                 tCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCtt         >  1:5459530/1‑93 (MQ=255)
                                                                                  cTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTa        <  2:2110009/93‑1 (MQ=255)
                                                                                  cTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTa        <  2:5903197/93‑1 (MQ=255)
                                                                                   ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTAt       <  1:3903266/93‑1 (MQ=255)
                                                                                   ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTAt       <  2:1317664/93‑1 (MQ=255)
                                                                                    tCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATc      <  1:5915956/93‑1 (MQ=255)
                                                                                     cccTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATcc     >  2:3218851/1‑93 (MQ=255)
                                                                                     cccTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATcc     >  2:502374/1‑93 (MQ=255)
                                                                                      ccTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATccc    >  1:6159412/1‑93 (MQ=255)
                                                                                      ccTCAAGAAAAAATCAGGACAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATcac    >  2:2020766/1‑91 (MQ=37)
                                                                                       cTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCg   >  1:1941861/1‑93 (MQ=255)
                                                                                             |                                                                                       
TGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGC  >  NC_000913/3099564‑3099743

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: