Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,786 T→A 7.5% E166D (GAA→GAT) ‡ feaR ← DNA‑binding transcriptional activator FeaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,446,7860TA 7.5% 254.7 / 6.4 85E166D (GAA→GAT) ‡feaRDNA‑binding transcriptional activator FeaR
Reads supporting (aligned to +/- strand):  ref base T (43/35);  new base A (5/2);  total (48/37)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAA  >  NC_000913/1446639‑1446933
                                                                                                                                                   |                                                                                                                                                   
gCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTc                                                                                                                                                   <  1:1007912/150‑1 (MQ=255)
   tAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGTCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGa                                                                                                                                                >  1:1462796/1‑150 (MQ=255)
   tAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGa                                                                                                                                                >  2:940554/1‑150 (MQ=255)
    aaTATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGaa                                                                                                                                               >  2:540355/1‑150 (MQ=255)
      tatCTCTTCGCGAATATTATCGTCTATAAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTACCGCTGATGAAGTACCGGGCGAAGCAGACACACCATCGCCTGTAGAGCAGCTTCACTTTCTGTTTCAGAAAg                                                                                                                                             >  2:1025371/1‑150 (MQ=255)
       atCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGAAACAGAAAGt                                                                                                                                            <  2:579648/150‑1 (MQ=255)
          tctTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGAAACAGAAAGTgaa                                                                                                                                         <  1:56899/150‑3 (MQ=255)
           ctTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCg                                                                                                                                        >  1:759535/1‑150 (MQ=255)
           ctTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCg                                                                                                                                        >  1:733420/1‑150 (MQ=255)
            ttCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGTGGTTGAACAGATTCCCGCTTATTAAGTATCGGGTGCAGCAGACACACCATAGGCTTTAGCTCAGCTTCACTTTCTGATTCAGTAAGTgaagc                                                                                                                                       <  2:990956/150‑5 (MQ=255)
                      ttatCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATg                                                                                                                             >  2:474995/1‑150 (MQ=255)
                           gTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCtctc                                                                                                                        <  1:1260684/150‑1 (MQ=255)
                           gTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCtctc                                                                                                                        >  1:168561/1‑150 (MQ=255)
                           gTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCtctc                                                                                                                        >  1:730473/1‑150 (MQ=255)
                            tCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTcc                                                                                                                       <  2:432759/150‑1 (MQ=255)
                               aTCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCACCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGt                                                                                                                    <  2:1088394/150‑1 (MQ=255)
                                tCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTa                                                                                                                   <  1:325777/150‑1 (MQ=255)
                                  aaCGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAAc                                                                                                                 <  1:179483/150‑1 (MQ=255)
                                  aaCGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAAc                                                                                                                 >  2:885203/1‑150 (MQ=255)
                                   aCGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACa                                                                                                                <  1:434946/150‑1 (MQ=255)
                                      aaaCCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGc                                                                                                             <  1:1406163/149‑1 (MQ=255)
                                        aCCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTcc                                                                                                                       >  1:863974/1‑138 (MQ=255)
                                        aCCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTcc                                                                                                                       <  2:863974/138‑1 (MQ=255)
                                        aCCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGa                                                                                                           <  2:1031301/150‑1 (MQ=255)
                                         ccACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGAt                                                                                                          >  2:1160473/1‑150 (MQ=255)
                                         ccACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGAt                                                                                                          <  2:746796/150‑1 (MQ=255)
                                           aCTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGa                                                                                                        >  1:1286243/1‑150 (MQ=255)
                                            cTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGAc                                                                                                       >  2:182272/1‑150 (MQ=255)
                                              ttttGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTg                                                                                                     <  1:324667/150‑1 (MQ=255)
                                              ttttGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTg                                                                                                     >  1:30153/1‑150 (MQ=255)
                                              ttttGAAACTGACGTTCACGACGAGGTTGAACAGATTCCAGCTGATGAAGTACCCGGCTCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCAGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACACGCGAATACTg                                                                                                     >  2:1178481/1‑150 (MQ=255)
                                                 tGAAAGTGACGTTCACGACGAGGTTGAACAGATTCCCGCTTATGAAGTACCGGGCGCTGCAGACACACCATCTCCTGTAGCGCAGCTTCACTTTCTGTTTCAGATAGTGCCGGATTATTCATTCTCTCCTGTAACAGGCGATGACTTAGt                                                                                                  <  2:625139/150‑1 (MQ=255)
                                                 tGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGt                                                                                                  >  2:723038/1‑150 (MQ=255)
                                                 tGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGt                                                                                                  >  1:1410946/1‑150 (MQ=255)
                                                 tGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGt                                                                                                  <  2:942676/150‑1 (MQ=255)
                                                  gAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGtt                                                                                                 <  2:759535/150‑1 (MQ=255)
                                                          cGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATg                                                                                         >  1:589735/1‑150 (MQ=255)
                                                            ttCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATggg                                                                                       <  2:733420/150‑1 (MQ=255)
                                                              cACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTa                                                                                     <  2:628689/150‑1 (MQ=255)
                                                               acgacgAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTaa                                                                                    >  2:1113462/1‑150 (MQ=255)
                                                                cgacgaGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAAc                                                                                                                 >  2:1342547/1‑120 (MQ=255)
                                                                cgacgaGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAAc                                                                                                                 <  1:1342547/120‑1 (MQ=255)
                                                                cgacgaGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAg                                                                                   >  2:791537/1‑150 (MQ=255)
                                                                         tGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTcc                                                                          >  1:1167882/1‑150 (MQ=255)
                                                                         tGAACAGATTCCCGCTGATGAAGTACCGGCCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGGGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGCCAGAGTTGCACAATGGGTAAGTCAGCGTcc                                                                          >  1:1379456/1‑150 (MQ=255)
                                                                          gAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCa                                                                         >  2:971844/1‑150 (MQ=255)
                                                                               gATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCtt                                                                    >  1:428911/1‑150 (MQ=255)
                                                                                 ttCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTc                                                                  >  1:1045685/1‑150 (MQ=255)
                                                                                 ttCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTc                                                                  >  1:1051027/1‑150 (MQ=255)
                                                                                    ccGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTATCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTgc                                                               >  1:665412/1‑150 (MQ=255)
                                                                                     cGCTGATGAAGTACCGTGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTgcg                                                              <  2:1139753/150‑1 (MQ=255)
                                                                                           tgaAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATa                                                        >  2:538919/1‑150 (MQ=255)
                                                                                           tgaAGTACAGAGCTCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATa                                                        >  2:93828/1‑150 (MQ=255)
                                                                                                  ccGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGttttt                                                 <  1:1113462/150‑1 (MQ=255)
                                                                                                   cGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTg                                                >  2:50561/1‑150 (MQ=255)
                                                                                                       cgcAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATgg                                            <  1:841655/150‑1 (MQ=255)
                                                                                                         cagcagACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATgggg                                          >  1:743064/1‑150 (MQ=255)
                                                                                                          agcagACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGa                                         <  1:23286/150‑1 (MQ=255)
                                                                                                                cacacCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATTGGGAAAatat                                   <  2:317678/150‑1 (MQ=255)
                                                                                                                cacacCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAatat                                   >  2:1118009/1‑150 (MQ=255)
                                                                                                                  cacCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAg                                                           >  1:243792/1‑124 (MQ=255)
                                                                                                                  cacCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAg                                                           <  2:243792/124‑1 (MQ=255)
                                                                                                                  cacCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTg                                 >  2:937252/1‑150 (MQ=255)
                                                                                                                   acCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGt                                <  1:955597/150‑1 (MQ=255)
                                                                                                                    ccATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGtt                               >  2:595156/1‑150 (MQ=255)
                                                                                                                     cATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTc                              >  2:44389/1‑150 (MQ=255)
                                                                                                                       tCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCca                            <  2:59840/150‑1 (MQ=255)
                                                                                                                        cGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCcag                           <  1:13061/150‑1 (MQ=255)
                                                                                                                        cGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCcag                           <  2:168561/150‑1 (MQ=255)
                                                                                                                         gCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCcagc                          <  1:885203/150‑1 (MQ=255)
                                                                                                                           cTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCcagcag                        >  2:1285504/1‑150 (MQ=255)
                                                                                                                                     gCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATg                                             <  2:621755/119‑1 (MQ=255)
                                                                                                                                     gCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATg                                             >  1:621755/1‑119 (MQ=255)
                                                                                                                                     gCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGc              >  1:764636/1‑150 (MQ=255)
                                                                                                                                     ccggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGc              >  1:1449578/5‑150 (MQ=255)
                                                                                                                                      cggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCa             >  2:1443830/4‑150 (MQ=255)
                                                                                                                                      cTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCa             <  1:1160473/150‑1 (MQ=255)
                                                                                                                                       ggcACTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCaa            >  2:406925/3‑150 (MQ=255)
                                                                                                                                          aCTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAATGCGTGGCAAAAg         <  1:401792/150‑1 (MQ=255)
                                                                                                                                          aCTTTCTGAAACAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAg         >  2:180447/1‑150 (MQ=255)
                                                                                                                                              tCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAAt     >  1:1504690/1‑150 (MQ=255)
                                                                                                                                                 gTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGaa  >  1:1075995/1‑150 (MQ=255)
                                                                                                                                                 gTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGaa  <  2:1045685/150‑1 (MQ=255)
                                                                                                                                                 gTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGaa  <  2:1051027/150‑1 (MQ=255)
                                                                                                                                                 gTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGaa  >  1:958425/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACTGACGTTCACGACGAGGTTGAACAGATTCCCGCTGATGAAGTACCGGGCGCAGCAGACACACCATCGCCTGTAGCGCAGCTTCACTTTCTGTTTCAGAAAGTGCCGGATTATTCATGCTCTCCTGTAACAGGCGATGACTGAGTTGCACCATGGGTAAGTCAGCGTCCAGTCTTTCTGCGCAGATAGGTTTTTGATGGGGAAAATATTGTTCCAGCAGAGTGCGTGGCAAAAGTAATGAA  >  NC_000913/1446639‑1446933

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: