Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,911,822 C→A 5.3% intergenic (‑297/+17) 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,8220CA 5.3% 405.2 / 9.7 114intergenic (‑297/+17)pstS/glmSphosphate ABC transporter periplasmic binding protein/L‑glutamine:D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (56/52);  new base A (2/4);  total (58/56)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

Only 100 of 115 total aligned reads displayed.

CAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTT  >  NC_000913/3911674‑3911963
                                                                                                                                                    |                                                                                                                                             
cAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAGCAGGAAGAAAAATGCCCCGCTTAc                                                                                                                                               >  2:584277/1‑149 (MQ=255)
 aGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACg                                                                                                                                              >  2:202272/1‑149 (MQ=255)
   tCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCa                                                                                                                                            >  4:364397/1‑149 (MQ=255)
    cAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAg                                                                                                                                           <  1:81433/149‑1 (MQ=255)
     aaTTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAgg                                                                                                                                          <  1:820782/149‑1 (MQ=255)
       ttACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGc                                                                                                                                        <  4:29000/149‑1 (MQ=255)
          cAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTGTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGATGGAAGGATGTTTGATTAAAAACATAACAGGCAGAAAAATGCCCCGCTTACGCAGGGCATc                                                                                                                                     >  2:98607/1‑149 (MQ=255)
           aaaCATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCTGCGTAAGCGGGGCAttt                                                                                                                                    <  1:404000/149‑3 (MQ=255)
           aaaCATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATcc                                                                                                                                    <  3:864124/149‑1 (MQ=255)
           aaaCATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATcc                                                                                                                                    >  4:396355/1‑149 (MQ=255)
           aaaCATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATcc                                                                                                                                    >  1:225485/1‑149 (MQ=255)
               aTTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCAttt                                                                                                                                >  1:19777/1‑149 (MQ=255)
                ttAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATtta                                                                                                                               <  3:124999/149‑1 (MQ=255)
                 taataaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATttat                                                                                                                              <  1:191394/149‑1 (MQ=255)
                 taataaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATttat                                                                                                                              <  4:272102/149‑1 (MQ=255)
                   ataaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATttatta                                                                                                                            <  3:81946/149‑1 (MQ=255)
                   ataaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATttatta                                                                                                                            <  2:300416/149‑1 (MQ=255)
                     aaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACt                                                                                                                          <  2:97168/149‑1 (MQ=255)
                     aaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACt                                                                                                                          >  3:689364/1‑149 (MQ=255)
                        gAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCaa                                                                                                                       >  1:101361/1‑149 (MQ=255)
                         aaGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAAc                                                                                                                      >  4:808482/1‑149 (MQ=255)
                              aTGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTaa                                                                                                                 >  2:59265/1‑149 (MQ=255)
                              aTGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTaa                                                                                                                 >  4:431978/1‑149 (MQ=255)
                                     aaaaaTTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGAttt                                                                                                          >  1:96689/1‑149 (MQ=255)
                                      aaaaTTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGAtttt                                                                                                         <  3:640504/149‑1 (MQ=255)
                                        aaTTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGc                                                                                                       <  2:572571/149‑1 (MQ=255)
                                        aaTTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGc                                                                                                       <  3:634571/149‑1 (MQ=255)
                                         aTTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGcc                                                                                                      >  2:521830/1‑149 (MQ=255)
                                          ttttCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCa                                                                                                     <  1:1578/149‑1 (MQ=255)
                                           tttCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAg                                                                                                    <  1:411848/149‑1 (MQ=255)
                                              cATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGtt                                                                                                 <  2:327824/149‑1 (MQ=255)
                                                ttCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGAGTAACAACATAATAGGCAGAAAAATGCCCCCCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTAc                                                                                               >  4:184678/1‑149 (MQ=255)
                                                       aCAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCTGCGTAAGCGGGGCAtttttc                                                                                                                                 <  2:564011/108‑6 (MQ=38)
                                                       aCAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCTGCGTAAGCGGGGCAtttttc                                                                                                                                 >  1:564011/1‑103 (MQ=38)
                                                       aCAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGCAAGAAAAATGCCCCGCTTACGCAGGGCGTCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTg                                                                                        >  3:554105/1‑149 (MQ=255)
                                                        cAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTgg                                                                                       <  1:602216/149‑1 (MQ=255)
                                                         agagAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGt                                                                                      <  1:637592/149‑1 (MQ=255)
                                                             aaaaaGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCGAGGTTACGCGGCTGGTCAAc                                                                                  >  4:170533/1‑149 (MQ=255)
                                                             aaaaaGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAAc                                                                                  >  4:665739/1‑149 (MQ=255)
                                                                aaGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTc                                                                               <  2:19777/149‑1 (MQ=255)
                                                                   tAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGt                                                                            <  2:562710/149‑1 (MQ=255)
                                                                        gAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCtt                                                                       >  4:556031/1‑149 (MQ=255)
                                                                        gAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCtt                                                                       >  2:765530/1‑149 (MQ=255)
                                                                         aaGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCttt                                                                      <  1:607948/149‑1 (MQ=255)
                                                                            aTGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGAt                                                                   <  3:788800/149‑1 (MQ=255)
                                                                            aTGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGAt                                                                   <  1:102651/149‑1 (MQ=255)
                                                                              gACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCa                                                                 <  1:266683/149‑1 (MQ=255)
                                                                               aCGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCAGCTGGTCAACGTCGGTACCTTTGTTCAg                                                                >  4:873314/1‑149 (MQ=255)
                                                                                       tCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACAt                                                        <  1:175087/149‑1 (MQ=255)
                                                                                        cacaTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATg                                                       <  4:3122/149‑1 (MQ=255)
                                                                                         acaTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATgg                                                      >  2:396004/1‑149 (MQ=255)
                                                                                         acaTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATgg                                                      >  1:561088/1‑149 (MQ=255)
                                                                                           aTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTa                                                    <  2:449681/149‑1 (MQ=255)
                                                                                            tGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTaa                                                   >  2:42354/1‑149 (MQ=255)
                                                                                            tGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTaa                                                   <  1:584277/149‑1 (MQ=255)
                                                                                            tGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGAGTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTaa                                                   <  2:170463/149‑1 (MQ=255)
                                                                                                 ttGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCcag                                              >  3:813626/1‑149 (MQ=255)
                                                                                                 ttGGCAGGATGTTTGAATAAAAACATTACAGGAAGAAATAAGCCCCGCTTACGCAGAGCATCCATTTATTACTCAACCGTATCCGATTTTGCTAGGTTACGCGGCTGGTCAACTTCGCTGCCTTTGATCAGCGCGCCATGGTAAGCcag                                              <  2:299286/149‑1 (MQ=255)
                                                                                                    gCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTaa                                                                                                                 >  3:805267/1‑79 (MQ=255)
                                                                                                    gCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTaa                                                                                                                 <  4:805267/79‑1 (MQ=255)
                                                                                                    gCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCcagcag                                           <  1:668969/149‑1 (MQ=255)
                                                                                                    gCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCcagcag                                           >  1:753189/1‑149 (MQ=255)
                                                                                                    gCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCcagcag                                           >  2:591822/1‑149 (MQ=255)
                                                                                                           gTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGc                                    <  3:314366/149‑1 (MQ=255)
                                                                                                            tttGATTAAAAACATAACAGGAAGAAAGATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCg                                   >  3:590918/1‑149 (MQ=255)
                                                                                                              tGACTACAACCATAACAGCAAGAAAAATGCCCCGCTTACGAAGGGCGTTAATTTATTTCTCAACCGTAACCGATTTTGCCTTGTTACGCGGCTGGTAAACGTCGGAGCCTTTGATCACCGCGACATGGTAAGCCAGAAGCTGCAGCGGa                                 >  2:257558/1‑149 (MQ=255)
                                                                                                                aTTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAAc                               >  1:697399/1‑149 (MQ=255)
                                                                                                                 ttAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACg                              <  4:68598/149‑1 (MQ=255)
                                                                                                                   aaaaaCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGgt                            <  2:711710/149‑1 (MQ=255)
                                                                                                                   aaaaaCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGgt                            <  1:257558/149‑1 (MQ=255)
                                                                                                                     aaaCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGgtgt                          >  4:351071/1‑149 (MQ=255)
                                                                                                                     aaaCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGgtgt                          >  3:698731/1‑149 (MQ=255)
                                                                                                                         aTAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTagaa                      >  3:668260/1‑149 (MQ=255)
                                                                                                                         aTAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTagaa                      >  2:109630/1‑149 (MQ=255)
                                                                                                                          tAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGTCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTagaag                     >  4:18832/1‑149 (MQ=255)
                                                                                                                          tAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTagaag                     >  1:586593/1‑149 (MQ=255)
                                                                                                                           aaCAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTagaaga                    >  3:267339/1‑149 (MQ=255)
                                                                                                                            aCAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGAt                   >  2:509321/1‑149 (MQ=255)
                                                                                                                             cAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATc                  <  1:765530/149‑1 (MQ=255)
                                                                                                                             cAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATc                  <  3:556031/149‑1 (MQ=255)
                                                                                                                              aGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCg                 <  3:808482/149‑1 (MQ=255)
                                                                                                                              aGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCg                 >  3:33020/1‑149 (MQ=255)
                                                                                                                              aGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCg                 >  2:259820/1‑149 (MQ=255)
                                                                                                                                gaagaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACAt                                                        >  2:774029/1‑108 (MQ=255)
                                                                                                                                gaagaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGt               <  3:130648/149‑1 (MQ=255)
                                                                                                                                gaagaaAAATGCCCCGCTTACGCAGGGCATCCATTTACTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGt               <  3:574478/149‑1 (MQ=255)
                                                                                                                                 aagaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATg                                                       <  3:91120/108‑1 (MQ=255)
                                                                                                                                 aagaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATg                                                       >  4:91120/1‑108 (MQ=255)
                                                                                                                                  agaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGc             >  3:399262/1‑149 (MQ=255)
                                                                                                                                  agaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGc             >  1:497147/1‑149 (MQ=255)
                                                                                                                                   gaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAgcgc                                                             <  3:520750/100‑1 (MQ=255)
                                                                                                                                   gaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAgcgc                                                             >  4:520750/1‑100 (MQ=255)
                                                                                                                                    atggaTGCCCTGCGTAAGCGGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCaa           >  3:522019/5‑149 (MQ=255)
                                                                                                                                     aaaaTGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAAt          >  1:722719/1‑149 (MQ=255)
                                                                                                                                       aaTGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCa        <  1:111921/149‑1 (MQ=255)
                                                                                                                                          gCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACct     <  1:400729/149‑1 (MQ=255)
                                                                                                                                            cccGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACctct   >  1:773641/1‑149 (MQ=255)
                                                                                                                                             ccGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCtt  >  4:371567/1‑149 (MQ=255)
                                                                                                                                             ccGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCtt  <  4:388982/149‑1 (MQ=255)
                                                                                                                                             ccGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCtt  >  3:707148/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
CAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTT  >  NC_000913/3911674‑3911963

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: