Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,296,382 +GC intergenic (+587/+55) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,296,3811.C99.0% 605.9 / ‑2.5 201intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (80/118);  total (80/121)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_bme_pgi4,296,3812.G99.0% 604.1 / NA 201intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  major base G (80/117);  minor base C (0/2);  total (80/121)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTT  >  NC_000913_3_bme_pgi/4296294‑4296520
                                                                                        ||                                                                                                                                            
agCGAGGCTTGCCGCGTCTTTTCAGGCCTAAGCCAGAAACAGCAATAGCCTGACTTCGCGTGATTTTGTAGGGCGTATAAGGCGTTTACCCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTaa                                                                                  <  2:404483/148‑1 (MQ=11)
 gctgcGCTCATCGCGTCTTATCCGGCCTACGCCAGACAGCGCACTAGCCTGATTGAGCGTCATTTTGTAGGTCGGACAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAAc                                                                                 <  2:422176/145‑1 (MQ=1)
 gCGACGCTTATCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTACGTCGTATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAAc                                                                                 <  1:386837/149‑1 (MQ=1)
    aCGCTTAACGCCTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGt                                                                              <  1:95352/149‑1 (MQ=255)
     cGCTTAACGCGTCTTATCAGGCCTACGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGtt                                                                             <  1:741974/149‑1 (MQ=255)
     cGCTTAACGCGTCTTATCAGGCCTACCGCAGGCAGCGGTATTGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGtt                                                                             <  2:13697/149‑1 (MQ=11)
      ggttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTg                                                                            <  1:121779/147‑1 (MQ=255)
      gCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTg                                                                            <  1:122849/149‑1 (MQ=255)
      gCTTAACGCGTCTTATCAGGCCTACGCCGGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTg                                                                            <  1:484250/149‑1 (MQ=255)
      gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTg                                                                            <  1:673223/149‑1 (MQ=255)
      gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTg                                                                            <  2:640284/149‑1 (MQ=255)
        ttAACGCGTCTTATCAGCCCTACGCCAGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGATAAGTCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGt                                                                          <  2:729718/149‑1 (MQ=11)
           aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTCTCTTAACCGTTGGTTaa                                                                       >  1:554626/1‑149 (MQ=255)
           aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGCTaa                                                                       >  2:657296/1‑149 (MQ=255)
             gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTCACAGGTGGTTAAtt                                                                     >  1:524889/1‑149 (MQ=11)
             gcgTCTGTTCAGGCGTACGGCATACGGCGAATTAGTCTGATTGAGCGTGATTTTGTGAGGCGGGCAAGGCGTTTACGCCGTGTCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAAtt                                                                     <  2:270883/149‑1 (MQ=9)
              cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCGGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAAttt                                                                    <  1:580394/149‑1 (MQ=255)
                tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAAttttt                                                                  >  2:733351/1‑149 (MQ=255)
                 cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTc                                                                 <  1:514355/149‑1 (MQ=255)
                 cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTc                                                                 <  2:335758/149‑1 (MQ=255)
                    aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTGGGACGTATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTTCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCaaa                                                              >  2:197213/1‑149 (MQ=14)
                    aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCaaa                                                              <  1:755121/149‑1 (MQ=255)
                     tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                             >  1:497165/1‑149 (MQ=255)
                     tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                             >  2:79524/1‑149 (MQ=255)
                      cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACTGTTGGTTAATTTTTCAAAct                                                            <  2:477260/149‑1 (MQ=255)
                       aGGCCTACGCCAGACGGCCCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                           <  1:421345/149‑1 (MQ=17)
                       aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                           >  1:290264/1‑149 (MQ=255)
                         gCCTACGCCAGACCGCGCAATAGCCCGCTTTCTGGCAGTTGTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                         >  1:185716/1‑149 (MQ=17)
                         gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                         <  1:295887/149‑1 (MQ=255)
                         gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                         <  1:208518/149‑1 (MQ=255)
                          ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                        >  2:376663/1‑149 (MQ=255)
                            tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                      >  2:636277/1‑149 (MQ=255)
                             aCGGCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca                                                     <  2:489667/149‑1 (MQ=255)
                              ccccAGACAGCGCAATACCCTAACTTAGCGTGATTTTGTAGGTCGGATAATGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac                                                    <  1:442639/147‑1 (MQ=17)
                               gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGCTGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAACCTCTTCCCAcc                                                   >  1:656714/1‑149 (MQ=17)
                                ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTGGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAAATCTTCACAccc                                                  >  1:377133/1‑149 (MQ=17)
                                 cAGACAGGGCAATAGCCTGATTTACCGTGATTGTGAAGCTCGGACAAGGCGTTTAGGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                 <  2:152208/149‑1 (MQ=18)
                                 cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                 >  1:765707/1‑149 (MQ=255)
                                 cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                 >  1:235778/1‑149 (MQ=255)
                                 cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                 <  1:594003/149‑1 (MQ=255)
                                  aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGt                                                                          >  2:562985/1‑123 (MQ=255)
                                   gACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGGTAAGGCGTTTACCCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                               <  2:559722/149‑1 (MQ=18)
                                    aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                              <  1:532789/149‑1 (MQ=255)
                                    aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                              >  1:426662/1‑149 (MQ=255)
                                     cAGCGCAAGAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt                                             <  1:14643/149‑1 (MQ=255)
                                         gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                         <  1:707294/149‑1 (MQ=255)
                                         gcAATAGCCTGATTTAGCGTGATTTTGGAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                         <  1:218407/149‑1 (MQ=255)
                                         gcAATAGCCTGATTTAGCGTGATGTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAACTTTTCAAACTCTTCACACCCGGTATCaaa                                         <  1:650875/149‑1 (MQ=18)
                                          cAATAGCCTGATTTAGCTGGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                        <  1:47989/149‑1 (MQ=18)
                                          cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAAc                                        >  1:675434/1‑149 (MQ=255)
                                           aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                       >  1:340790/1‑149 (MQ=255)
                                               gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                         <  1:679252/127‑1 (MQ=255)
                                               gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                         >  2:679252/1‑127 (MQ=255)
                                               gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAACCCCTTc                                   >  2:464921/1‑149 (MQ=255)
                                               gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                   <  1:635825/149‑1 (MQ=255)
                                                ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                  <  2:593442/149‑1 (MQ=255)
                                                ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAAGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                  >  1:346257/1‑149 (MQ=255)
                                                ccTGATTTAGCGTGATTTTGTACGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                  >  2:666545/1‑149 (MQ=255)
                                                   ggtctAGCGTGATTTTGTAGGTCGGATAAGGCCTTTACGCCACATCCGACAGGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                               <  1:56704/145‑1 (MQ=21)
                                                   gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                               <  1:438258/149‑1 (MQ=255)
                                                   gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                               <  1:292555/149‑1 (MQ=255)
                                                    aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                              >  1:519589/1‑149 (MQ=255)
                                                     tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTTAACCGCTGCACACCCGGTATCAAACCCTTCCATACa                             <  1:185395/149‑1 (MQ=18)
                                                     tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGCATCAAACCCTTCCATACa                             >  1:768209/1‑149 (MQ=255)
                                                        aGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                          <  1:112896/149‑1 (MQ=255)
                                                         gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                         <  1:256103/149‑1 (MQ=255)
                                                           gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTACTTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                       >  1:401358/1‑149 (MQ=255)
                                                           gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                       <  1:641227/149‑1 (MQ=255)
                                                           gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                       <  1:225438/149‑1 (MQ=255)
                                                           gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTCGCCGCGTCTTATCAGG‑CCTAACTTAACCGTAGGTAAATTTTTCAAACTCTTCACACACGGTAACAAACCCTTCCATACAGCTCAg                       >  1:285106/1‑149 (MQ=18)
                                                            tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                               <  1:209041/140‑1 (MQ=255)
                                                            tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                      <  2:274141/149‑1 (MQ=255)
                                                            tGATTTTGGAGGCCGGATAAGGCGTTTACGCCGCACCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCGTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                      <  1:137312/149‑1 (MQ=25)
                                                            tGATTTGGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                      <  1:458642/149‑1 (MQ=255)
                                                               ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTc                   <  1:89230/149‑1 (MQ=255)
                                                                  tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                <  2:599755/149‑1 (MQ=255)
                                                                  tGTACGTCGGATAAGGCGTTTACGCCGCATCCCACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                <  1:589492/149‑1 (MQ=25)
                                                                   gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                      <  1:312848/142‑1 (MQ=255)
                                                                    tAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa              <  2:169206/149‑1 (MQ=255)
                                                                     tggTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAg             <  1:113957/148‑1 (MQ=255)
                                                                      ggTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt            <  1:714792/149‑1 (MQ=255)
                                                                      ggTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTAGCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACCGCTCCGGCTCAACCAgc            >  2:298145/1‑148 (MQ=21)
                                                                      ggTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACACTTGCCGCGTCTTATCAGG‑CATCTCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACACTTCCATACAGCTAAGGTTCACCCAGt            >  2:52990/1‑149 (MQ=25)
                                                                       gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAGTTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAAAGCTCAGGTTCAACCAGTg           >  2:213609/1‑149 (MQ=25)
                                                                       gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                  <  2:738852/110‑1 (MQ=255)
                                                                       gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg           <  2:254647/149‑1 (MQ=255)
                                                                       gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg           >  1:607795/1‑149 (MQ=255)
                                                                        tCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCAGCGTCTTA‑CAGGACCTATCTTCACAGTGCGGTAATTTTTCCAACTCTTCACTCCGGGGATCAAACCCGTCCGTACGGCTCAGGTTCAACCAGTGc          >  1:374242/1‑149 (MQ=12)
                                                                           gATAAGGCGTTTACGCCGCCGCCGACATTAATGTCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGCTCAACCAGTGCAgc       <  2:129734/149‑1 (MQ=37)
                                                                           gATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCATTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc       >  2:278245/1‑149 (MQ=255)
                                                                            aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAAGCGTTGGTTAATTTTTCCAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCACTGCAccg      >  2:77860/1‑149 (MQ=25)
                                                                            aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGACGCGACGCTTGCCGCGTCTTATCAGG‑CATATCTTAACCGCTGGTTCATTTTTCAAACTCTTCACACCCGGTATCACACCCGTCCATACAGCGCAGGTTCAACCAGTGCAgcg      >  1:486567/1‑149 (MQ=25)
                                                                              aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                         >  2:105481/1‑112 (MQ=255)
                                                                              aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                         <  1:105481/112‑1 (MQ=255)
                                                                              aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAGACCCTTCCATACAGCTAAGGTTCAACCACTCCAGCGcc    >  2:230451/1‑149 (MQ=25)
                                                                              aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGcc    <  1:552921/149‑1 (MQ=255)
                                                                              aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGcc    <  1:724378/149‑1 (MQ=255)
                                                                              aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGcc    >  2:31057/1‑149 (MQ=255)
                                                                               aGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCt   <  1:349421/149‑1 (MQ=255)
                                                                                ggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCtt  <  2:320959/149‑1 (MQ=255)
                                                                                        ||                                                                                                                                            
TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGG‑CCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTT  >  NC_000913_3_bme_pgi/4296294‑4296520

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