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% 313.4 / 4.7 91intergenic (‑297/+17)pstS/glmSphosphate ABC transporter periplasmic binding protein/L‑glutamine:D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (29/56);  new base A (2/3);  total (32/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACAT  >  NC_000913/3911674‑3911969
                                                                                                                                                    |                                                                                                                                                   
cAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTAc                                                                                                                                                     >  4:320631/1‑149 (MQ=255)
   tCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCa                                                                                                                                                  >  3:166076/1‑149 (MQ=255)
   tCAATTACAAACATAAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACCAGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTAGGCa                                                                                                                                                  >  4:331872/1‑149 (MQ=255)
    cAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAg                                                                                                                                                 <  1:483259/149‑1 (MQ=255)
     aaTTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAgg                                                                                                                                                <  2:111261/149‑1 (MQ=255)
     aaTTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAgg                                                                                                                                                <  3:95015/149‑1 (MQ=255)
     aaTTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAgg                                                                                                                                                <  1:230094/149‑1 (MQ=255)
     aaTTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAgg                                                                                                                                                <  4:166076/149‑1 (MQ=255)
             aCATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCAt                                                                                                                                        <  1:662729/149‑1 (MQ=255)
              cATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCTGCGTAAGCGGGGCAtttttc                                                                                                                                       <  1:249461/149‑6 (MQ=255)
                 taataaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATttat                                                                                                                                    <  3:128099/149‑1 (MQ=255)
                   ataaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGATGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATttatta                                                                                                                                  <  1:518546/149‑1 (MQ=255)
                      aCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTc                                                                                                                               >  1:446934/1‑149 (MQ=255)
                        gAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCaa                                                                                                                             >  3:14593/1‑149 (MQ=255)
                         aaGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAAc                                                                                                                            <  4:262436/149‑1 (MQ=255)
                           gagaTGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCg                                                                                                                          <  1:561598/149‑1 (MQ=255)
                           gagaTGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCg                                                                                                                          >  1:14536/1‑149 (MQ=255)
                                    gAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGAtt                                                                                                                 >  1:166598/1‑149 (MQ=255)
                                    gAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGAtt                                                                                                                 >  2:448164/1‑149 (MQ=255)
                                     aaaaaTTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGAttt                                                                                                                >  1:159982/1‑149 (MQ=255)
                                     aaaaaTTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGAttt                                                                                                                >  2:366544/1‑149 (MQ=255)
                                      aaaaTTTTTATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGAtttt                                                                                                               <  2:551514/149‑1 (MQ=255)
                                          ttttCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCa                                                                                                           <  4:543963/149‑1 (MQ=255)
                                          ttttCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCa                                                                                                           <  2:661977/149‑1 (MQ=255)
                                          ttttCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCa                                                                                                           <  1:646285/149‑1 (MQ=255)
                                             tCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGt                                                                                                        <  1:448164/149‑1 (MQ=255)
                                              cATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGtt                                                                                                       <  3:345123/149‑1 (MQ=255)
                                               aTTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTa                                                                                                      >  1:228436/1‑149 (MQ=255)
                                                  cTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTAcgc                                                                                                   <  4:359222/149‑1 (MQ=255)
                                                     tgACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGCAACATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGc                                                                                                <  2:325528/149‑1 (MQ=255)
                                                     tgACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGAttt                                                                                                                >  1:229725/1‑133 (MQ=255)
                                                     tgACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGAttt                                                                                                                <  2:229725/133‑1 (MQ=255)
                                                     tgACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGc                                                                                                <  2:485550/149‑1 (MQ=255)
                                                      gACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCt                                                                                               <  1:386698/149‑1 (MQ=255)
                                                       aCAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTg                                                                                              <  2:603763/149‑1 (MQ=255)
                                                         agagAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGt                                                                                            >  2:516693/1‑149 (MQ=255)
                                                            gaAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCaa                                                                                         <  2:166598/149‑1 (MQ=255)
                                                                 aGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTaa                                                                                                                       <  3:214845/114‑1 (MQ=255)
                                                                 aGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTaa                                                                                                                       >  4:214845/1‑114 (MQ=255)
                                                                 aGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCg                                                                                    <  2:228436/149‑1 (MQ=255)
                                                                 aGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCg                                                                                    <  2:61871/149‑1 (MQ=255)
                                                                 aGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCg                                                                                    <  3:673153/149‑1 (MQ=255)
                                                                  gTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCgg                                                                                   <  2:195819/149‑1 (MQ=255)
                                                                      ccGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGcc                                                                               >  3:45067/1‑149 (MQ=255)
                                                                       cGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAATTGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCt                                                                              >  4:392373/1‑149 (MQ=255)
                                                                         aaGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCttt                                                                            <  3:299225/149‑1 (MQ=255)
                                                                          aGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTg                                                                           <  4:88122/149‑1 (MQ=255)
                                                                           gATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGa                                                                          <  4:428088/149‑1 (MQ=255)
                                                                              gACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCa                                                                       <  4:26344/149‑1 (MQ=255)
                                                                              gACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCa                                                                       <  4:673522/149‑1 (MQ=255)
                                                                                  gTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAgcgc                                                                   <  3:285535/149‑1 (MQ=255)
                                                                                   tttGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGcgcg                                                                  >  1:627424/1‑149 (MQ=255)
                                                                                      gTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCAACa                                                               <  1:27498/149‑1 (MQ=255)
                                                                                        cacaTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATg                                                             <  2:270385/149‑1 (MQ=255)
                                                                                         acaTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATgg                                                            <  3:127020/149‑1 (MQ=255)
                                                                                         acaTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATgg                                                            >  1:267815/1‑149 (MQ=255)
                                                                                          caTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGt                                                           <  1:492935/149‑1 (MQ=255)
                                                                                           aTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTa                                                          >  4:290419/1‑149 (MQ=255)
                                                                                           aTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTa                                                          <  3:552121/149‑1 (MQ=255)
                                                                                           aTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTa                                                          <  4:487419/149‑1 (MQ=255)
                                                                                           aTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTa                                                          >  1:131863/1‑149 (MQ=255)
                                                                                            tGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTaa                                                         <  4:90274/149‑1 (MQ=255)
                                                                                               aGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGcc                                                      <  4:147053/149‑1 (MQ=255)
                                                                                                       ggATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTg                                                                           <  2:508281/120‑1 (MQ=255)
                                                                                                       ggATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTg                                                                           >  1:508281/1‑120 (MQ=255)
                                                                                                            tttGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCg                                         >  3:146617/1‑149 (MQ=255)
                                                                                                                aTTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAAc                                     <  2:106285/149‑1 (MQ=255)
                                                                                                                    aaaaCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTaa                                                         <  2:513044/125‑1 (MQ=255)
                                                                                                                    aaaaCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTaa                                                         >  1:513044/1‑125 (MQ=255)
                                                                                                                     aaaCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGgtgt                                <  1:344955/149‑1 (MQ=255)
                                                                                                                      aaCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGAt                                                                         <  4:365274/107‑1 (MQ=255)
                                                                                                                      aaCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGAt                                                                         >  3:365274/1‑107 (MQ=255)
                                                                                                                         aTAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTagaa                            <  1:366544/149‑1 (MQ=255)
                                                                                                                            aCAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGAt                         <  3:290419/149‑1 (MQ=255)
                                                                                                                             cAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATc                        >  2:606001/1‑149 (MQ=255)
                                                                                                                              aGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTg                                                                                              >  3:565292/1‑78 (MQ=255)
                                                                                                                              aGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTg                                                                                              <  4:565292/78‑1 (MQ=255)
                                                                                                                                gaagaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGt                     <  4:74196/149‑1 (MQ=255)
                                                                                                                                gaagaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGt                     <  2:159982/149‑1 (MQ=255)
                                                                                                                                gaagaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGt                     >  3:46274/1‑149 (MQ=255)
                                                                                                                                gaagaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGt                     >  3:318488/1‑149 (MQ=255)
                                                                                                                                  agaaAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGc                   >  2:649608/1‑149 (MQ=255)
                                                                                                                                   aatggaTGCCCTGCGTAAGCGGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTaa                                                         <  1:75545/105‑1 (MQ=255)
                                                                                                                                   aatggaTGCCCTGCGTAAGAGGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTaa                                                         >  2:75545/6‑110 (MQ=255)
                                                                                                                                    atggaTGCCCTGCGTAAGCGGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGc                                             <  4:316854/117‑1 (MQ=255)
                                                                                                                                    atggaTGCCCTGCGTAAGCGGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGc                                             >  3:316854/5‑121 (MQ=255)
                                                                                                                                     aaaaTGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAAt                <  1:516693/149‑1 (MQ=255)
                                                                                                                                       aaTGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGATAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCa              <  3:318618/149‑1 (MQ=255)
                                                                                                                                        aTGCCCCGCTTACACAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCAc             <  4:580879/149‑1 (MQ=255)
                                                                                                                                            cccGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACctct         >  3:244077/1‑149 (MQ=255)
                                                                                                                                               gCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTcc      >  4:509062/1‑149 (MQ=255)
                                                                                                                                                   aCGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGAGCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACAt  <  4:617994/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGTCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACAT  >  NC_000913/3911674‑3911969

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