Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,911,825 A→G 5.5% intergenic (‑300/+14) pstS ← / ← glmS phosphate ABC transporter periplasmic binding protein/L‑glutamine:D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,911,8250AG 5.5% 270.2 / 5.3 90intergenic (‑300/+14)pstS/glmSphosphate ABC transporter periplasmic binding protein/L‑glutamine:D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (45/40);  new base G (3/2);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

GTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGC  >  NC_000913/3911676‑3911974
                                                                                                                                                     |                                                                                                                                                     
gTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCa                                                                                                                                                       <  1:1134569/150‑1 (MQ=255)
gTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCa                                                                                                                                                       <  1:1266764/150‑1 (MQ=255)
  ctaTTACAAACATTAATATCGAAGAGATGACAGAAAAATTTTCATTCTGTTACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAgg                                                                                                                                                     <  1:254976/148‑1 (MQ=255)
   aaTTCCAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTACCCGTAGATGACTTTTTTTCACATGGAGTTGGCAGGATGGTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAggg                                                                                                                                                    <  1:244335/150‑1 (MQ=255)
   aaTTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAggg                                                                                                                                                    <  2:1289786/150‑1 (MQ=255)
   aaTTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGAGGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAggg                                                                                                                                                    <  2:189934/150‑1 (MQ=255)
     ttACAAACATTAATAACGAAGAGATGACAGAAATATTTTCATTCTGTGACAGAGTAAAATTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAACTGCCCCGCTTACGCAGGGCa                                                                                                                                                  <  1:72653/150‑1 (MQ=255)
      aaCAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGCAGTTGGCAGGATGTTTGATTAAAACCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCAt                                                                                                                                                 <  2:1265242/149‑1 (MQ=255)
       aCAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCTGCGTAAGCGGGGCAtt                                                                                                                                                <  2:851265/150‑2 (MQ=255)
       aCAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATc                                                                                                                                                >  1:239887/1‑150 (MQ=255)
        cAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACTGAGAAAAAGTAGCCGAAGATGACGGGTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATcc                                                                                                                                               >  1:646188/1‑150 (MQ=255)
        cAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATcc                                                                                                                                               <  1:286112/150‑1 (MQ=255)
        aaaaCATTAATAACGAAGAGATGACAGAAAAATTTTCATTCAGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATcc                                                                                                                                               <  1:1262821/149‑1 (MQ=255)
            cATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCAttt                                                                                                                                           >  1:672461/1‑150 (MQ=255)
              ttAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATttat                                                                                                                                         >  2:1035082/1‑150 (MQ=255)
               taataaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATttatt                                                                                                                                        >  1:361099/1‑150 (MQ=255)
                  taaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACt                                                                                                                                     <  1:361266/150‑1 (MQ=255)
                  aaaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACt                                                                                                                                     <  2:1204830/149‑1 (MQ=255)
                   aaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTc                                                                                                                                    >  1:706267/1‑150 (MQ=255)
                   aaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTc                                                                                                                                    <  1:362803/150‑1 (MQ=255)
                     cGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCaa                                                                                                                                  >  1:507016/1‑150 (MQ=255)
                      gAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAAc                                                                                                                                 <  2:163681/150‑1 (MQ=255)
                        agagATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCg                                                                                                                               >  2:486974/1‑150 (MQ=255)
                           gaTGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTaa                                                                                                                            >  1:1234834/1‑150 (MQ=255)
                               aCAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGa                                                                                                                        >  1:709019/1‑150 (MQ=255)
                               aCAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGa                                                                                                                        >  2:1159242/1‑150 (MQ=255)
                                 aGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGAtt                                                                                                                      >  2:1261525/1‑150 (MQ=255)
                                  gAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATCTTTGATTAAAAAAATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATATATTACTCAACCGTAACAGAttt                                                                                                                     >  2:1279305/1‑150 (MQ=255)
                                  gAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGAATGTCACATGGAGTTGGCAGTATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCAATTTATTACTCAACCGTAACCGAttt                                                                                                                     >  2:1274298/1‑150 (MQ=255)
                                        ttttCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAg                                                                                                               <  1:602889/150‑1 (MQ=255)
                                        ttttCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAg                                                                                                               <  2:1086697/150‑1 (MQ=255)
                                        ttttCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAg                                                                                                               <  2:801504/150‑1 (MQ=255)
                                                cTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCAATTATTACTCAACCGTAACCGATTTTGCCAGGTTAcgcg                                                                                                       >  1:483382/1‑150 (MQ=255)
                                                   tgACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCt                                                                                                    >  1:307472/1‑150 (MQ=255)
                                                   tgACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCt                                                                                                    >  1:447260/1‑150 (MQ=255)
                                                     aCAGAGAAAAAGAAGCCGAAGATGACGGTATGTCACATGGAGTTGGCAGGATGTTTTATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTgg                                                                                                  <  2:333034/150‑1 (MQ=255)
                                                      cAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGt                                                                                                 >  1:1079247/1‑150 (MQ=255)
                                                       agagAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTc                                                                                                <  1:1079960/150‑1 (MQ=255)
                                                         agaAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCaa                                                                                              >  1:401518/1‑150 (MQ=255)
                                                            aaaaGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTTTTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGGCAACGt                                                                                           >  1:865374/1‑150 (MQ=255)
                                                                gTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGt                                                                                       >  2:422133/1‑150 (MQ=255)
                                                                 tAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTg                                                                                      <  2:875537/150‑1 (MQ=255)
                                                                 tAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTg                                                                                      <  2:1040987/150‑1 (MQ=255)
                                                                   gCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGcc                                                                                    <  2:1123369/150‑1 (MQ=255)
                                                                    ccGAACATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGAGTAAAAACATAACCGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCt                                                                                   <  1:655337/150‑1 (MQ=255)
                                                                      gAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCttt                                                                                 >  2:1148195/1‑150 (MQ=255)
                                                                          aTGACGGTTTGTCACATGGAGTGGGCAGGATGGTTGATTAAAAACATAACAGGAAGAAAACTGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGGTACGCGGCTGGTCAACGTCGGTGCCTTTGATc                                                                             <  2:1156877/150‑1 (MQ=255)
                                                                             aCGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAgc                                                                          >  1:668919/1‑150 (MQ=255)
                                                                             aCGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAgc                                                                          >  1:281403/1‑150 (MQ=255)
                                                                                   tGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACa                                                                    <  1:744210/150‑1 (MQ=255)
                                                                                    gTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCAGCCATTTATTACTCCACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACAt                                                                   >  1:753746/1‑150 (MQ=255)
                                                                                        caTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAACTGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTa                                                               <  1:227346/150‑1 (MQ=255)
                                                                                        caTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTa                                                               <  1:609780/150‑1 (MQ=255)
                                                                                         aTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTaa                                                              <  2:991269/150‑1 (MQ=255)
                                                                                          tGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATGCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAg                                                             <  2:121965/150‑1 (MQ=255)
                                                                                            gAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGcc                                                           >  1:1071850/1‑150 (MQ=255)
                                                                                            gAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGcc                                                           <  1:1301734/150‑1 (MQ=255)
                                                                                               ttGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCcagc                                                        <  1:755175/150‑1 (MQ=255)
                                                                                                tGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCcagca                                                       <  1:613923/150‑1 (MQ=255)
                                                                                                tGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACACAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCcagca                                                       >  2:39580/1‑150 (MQ=255)
                                                                                                 ggCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCACCCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCATTTGATCAGCTCGACATGGTAAGCcagcat                                                      >  2:662140/1‑149 (MQ=255)
                                                                                                  gCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCagcagc                                                     <  2:706267/150‑1 (MQ=255)
                                                                                                   cAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCt                                                    >  1:1214096/1‑150 (MQ=255)
                                                                                                    aGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCATGGCATACATTTAATACTCAAAAGTAACCGATTTTACAATGTAACTCGGCTGGTCAAAGTCGGTGCCTTTTATCAGCGCAACATGGTAAGCCAGCAGCTg                                                   >  1:221582/1‑150 (MQ=255)
                                                                                                        tGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGc                                               <  2:296217/150‑1 (MQ=255)
                                                                                                             gATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCCTCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGCTCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAAc                                          >  1:879349/1‑150 (MQ=255)
                                                                                                                tAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGgt                                       >  1:688606/1‑150 (MQ=255)
                                                                                                                 aaaaaCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGgtg                                      >  2:260513/1‑150 (MQ=255)
                                                                                                                 aaaaaCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCCCGACATGGTAACCCAGCAGCTGCAGCGGAACGgtg                                      >  1:222261/1‑150 (MQ=255)
                                                                                                                          aCAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTACCCGATGTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCGGCGCGCCAAGGTAAGCCAGCAGCTGCAGCGGACCGGTGTAGAAGCTc                             >  1:209983/1‑150 (MQ=255)
                                                                                                                          aCAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATc                             <  1:1159242/150‑1 (MQ=255)
                                                                                                                           cAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTACGCCAGCAGCTGCAGCGGAACGGTGTAGAAGTTCg                            >  2:1023020/1‑150 (MQ=255)
                                                                                                                            aGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCgg                           >  2:322599/1‑150 (MQ=255)
                                                                                                                             ggAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGt                          <  1:1261525/150‑1 (MQ=255)
                                                                                                                              gaagaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGTTCGGTg                         >  2:264974/1‑150 (MQ=255)
                                                                                                                              gaagaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTg                         <  2:765616/150‑1 (MQ=255)
                                                                                                                              gaagaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGAAAGGTg                         >  2:295118/1‑150 (MQ=255)
                                                                                                                               aagaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGc                        <  1:986837/150‑1 (MQ=255)
                                                                                                                                agaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCa                       >  1:168296/1‑150 (MQ=255)
                                                                                                                                 gaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCCACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCaa                      >  1:646709/1‑150 (MQ=255)
                                                                                                                                    aaaTGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCTGCGGAACGGTGTAGAAGATCGGTGCAATcc                   >  1:573311/1‑149 (MQ=255)
                                                                                                                                     gaTGCCCTGCGTAAGCGGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTag                                   <  1:330559/132‑1 (MQ=255)
                                                                                                                                     gaTGCCCTGCGTAAGCGGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTag                                   >  2:330559/2‑133 (MQ=255)
                                                                                                                                     gaTGCCCTGCGTAAGCGGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCAc                  >  2:784425/2‑150 (MQ=255)
                                                                                                                                       tGCCCTGCGTAAGCGGGGCATCCATTTATTACTCACCCGTAACCGATTTTGCCAGGTTACGCGGCGGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGACCGGTGTAGAAGCTCGGTGGAATcccca                >  1:164487/1‑149 (MQ=255)
                                                                                                                                         ccccGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTa                                                               <  2:390916/101‑1 (MQ=255)
                                                                                                                                         ccccGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTa                                                               >  1:390916/1‑101 (MQ=255)
                                                                                                                                            cGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTcc           <  1:1257771/150‑1 (MQ=255)
                                                                                                                                            cGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTcc           >  1:729119/1‑150 (MQ=255)
                                                                                                                                            cGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTcc           <  1:1035082/150‑1 (MQ=255)
                                                                                                                                              cTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGAAATCACCTCTTCcac         >  2:669964/1‑150 (MQ=255)
                                                                                                                                                 aCGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATg      <  1:247460/150‑1 (MQ=255)
                                                                                                                                                     aGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGc  <  2:1084816/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGC  >  NC_000913/3911676‑3911974

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