Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,910,368 G→C 5.2% P11A (CCG→GCG)  pstC ← phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,910,3680GC 5.2% 392.2 / 7.4 125P11A (CCG→GCG) pstCphosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (58/60);  new base C (5/2);  total (63/62)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 127 total aligned reads displayed.

TCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAAT  >  NC_000913/3910224‑3910517
                                                                                                                                                |                                                                                                                                                     
tCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGt                                                                                                                                                  <  1:1582977/150‑1 (MQ=255)
 ccATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGtt                                                                                                                                                 >  2:1641217/1‑150 (MQ=255)
 ccATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGtt                                                                                                                                                 <  2:56871/150‑1 (MQ=255)
    tAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTaaa                                                                                                                                              >  1:1174577/1‑150 (MQ=255)
      ggAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAg                                                                                                                                            >  1:1451289/1‑150 (MQ=255)
       gAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGGCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTaaaatc                                                                                                                                           >  1:279742/1‑148 (MQ=255)
          aGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAgg                                                                                                                                        >  1:431374/1‑150 (MQ=255)
          aGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTAACCCACCGGGTAAAAAgggcg                                                                                                                                        <  2:1192829/150‑5 (MQ=255)
           gCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGc                                                                                                                                       >  2:758660/1‑150 (MQ=255)
           gCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGc                                                                                                                                       <  1:1573904/150‑1 (MQ=255)
             cAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCtt                                                                                                                                     >  1:1195322/1‑150 (MQ=255)
              aGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACCTCAATAGCACAATCAGAGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTg                                                                                                                                    <  2:1042411/150‑1 (MQ=255)
              aGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTg                                                                                                                                    >  1:1614321/1‑150 (MQ=255)
                aCCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATATCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCTGTGGGTTAAAAGCAGGCTTGGt                                                                                                                                  <  2:547692/150‑1 (MQ=255)
                aCCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGGTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGt                                                                                                                                  <  2:362634/150‑1 (MQ=255)
                  cAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTg                                                                                                                                <  2:716931/150‑1 (MQ=255)
                   aaaTTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGc                                                                                                                               <  2:1262368/150‑1 (MQ=255)
                   aaaTTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGGCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGc                                                                                                                               <  1:558826/150‑1 (MQ=255)
                      tttCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGc                                                                                                                            <  2:1016671/150‑1 (MQ=255)
                       ttCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGcc                                                                                                                           <  1:1550879/150‑1 (MQ=255)
                         cTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCAt                                                                                                                         >  2:570547/1‑150 (MQ=255)
                           gAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATaa                                                                                                                       <  1:1319775/150‑1 (MQ=255)
                            aaTGCACGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACCTCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAg                                                                                                                      <  2:100527/150‑1 (MQ=255)
                               gCTAGGCCAGTAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTa                                                                                                                   <  1:1576748/150‑1 (MQ=255)
                                cTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTAc                                                                                                                  >  1:859657/1‑150 (MQ=255)
                                 tCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTAct                                                                                                                 <  2:378968/150‑1 (MQ=255)
                                   ggCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGGTGCAGCCATATGTTActct                                                                                                               >  2:100939/1‑150 (MQ=255)
                                   cgCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTActct                                                                                                               >  1:223237/2‑150 (MQ=255)
                                    gCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAATTTACTCtt                                                                                                              <  2:1125188/150‑1 (MQ=255)
                                    gCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCtt                                                                                                              >  1:191124/1‑150 (MQ=255)
                                    gCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCtt                                                                                                              >  2:1660802/1‑150 (MQ=255)
                                       aggaggAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAg                                                                                                           <  2:484900/150‑1 (MQ=255)
                                       aggaggAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAg                                                                                                           <  1:1415971/150‑1 (MQ=255)
                                          aggagATGATCAGAGAGACAATAATGCCTCCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTa                                                                                                        >  2:320967/1‑150 (MQ=255)
                                          aggagATCATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTa                                                                                                        <  1:1039153/150‑1 (MQ=255)
                                             agaTGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAAc                                                                                                     >  1:1017911/1‑150 (MQ=255)
                                                tgatCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAAcgcg                                                                                                  >  1:78513/1‑150 (MQ=255)
                                                 gatCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGt                                                                                                 <  2:787421/150‑1 (MQ=255)
                                                 gatCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGt                                                                                                 >  1:1275502/1‑150 (MQ=255)
                                                  atCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTCCCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGtt                                                                                                >  1:211047/1‑150 (MQ=255)
                                                  atCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGtt                                                                                                <  2:1726198/150‑1 (MQ=255)
                                                   tCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGttt                                                                                               <  2:446049/150‑1 (MQ=255)
                                                   gcAGTGAGAAAATAATGCCACACAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCATCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGttt                                                                                               <  2:1507476/149‑1 (MQ=255)
                                                     agagagACAATAATGGCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGGTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTAc                                                                                             <  2:682519/150‑1 (MQ=255)
                                                     agagagACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTAc                                                                                             >  2:551152/1‑150 (MQ=255)
                                                     agagagACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTAAGTTAAACGCGTATAc                                                                                             >  2:1531594/1‑150 (MQ=255)
                                                          gaCAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGt                                                                                        >  2:421184/1‑150 (MQ=255)
                                                                    ccACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTaa                                                                              <  1:488267/150‑1 (MQ=255)
                                                                      aCCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATg                                                                            <  1:565055/150‑1 (MQ=255)
                                                                       cccAACATCAATAGCACAATCAGCGCCGCCTGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGc                                                                           <  2:311133/150‑1 (MQ=255)
                                                                       cccAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGc                                                                           >  1:482722/1‑150 (MQ=255)
                                                                         ccaCATCAATAGCACACTCAGCTCCGCCATTTTTAACCGCACGCTGAAAATTCTTGCTCCCTTTTTACCCGCGTGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCg                                                                         <  2:858739/148‑1 (MQ=255)
                                                                         cAACATCAATAGCACAATTAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCAGTTTGTACTGCCg                                                                         >  2:1001642/1‑150 (MQ=255)
                                                                         cAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCTGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCg                                                                         <  2:1410940/150‑1 (MQ=255)
                                                                               cAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCAAAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTTATGCCGGATgcg                                                                   >  2:131976/1‑150 (MQ=255)
                                                                                    gCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTa                                                              >  2:819348/1‑150 (MQ=255)
                                                                                        aaTCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTATAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAac                                                          <  1:1084304/150‑1 (MQ=255)
                                                                                        aaTCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAac                                                          >  1:15323/1‑150 (MQ=255)
                                                                                        aaTCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAac                                                          >  2:1268964/1‑150 (MQ=255)
                                                                                        aaTCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAac                                                          >  1:918876/1‑150 (MQ=255)
                                                                                           cAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACAcc                                                       >  2:1066258/1‑150 (MQ=255)
                                                                                             gagcctccAGTTTTAACAGCACGCTGAAAATTATGTCTCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAACCGCGTTTACAAAGTCATTTTGTAATGCCGGTTGCGGCGTAAATCACCGt                                                     <  2:839254/148‑1 (MQ=255)
                                                                                               gccgccAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTAc                                                   >  1:395008/1‑150 (MQ=255)
                                                                                                 cgccAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTAccc                                                 <  1:1453730/150‑1 (MQ=255)
                                                                                                    cAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCa                                                                                                            <  1:1554371/88‑1 (MQ=255)
                                                                                                    cAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCa                                                                                                            >  2:1554371/1‑88 (MQ=255)
                                                                                                     aGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGATACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGcc                                             >  2:1041757/1‑150 (MQ=255)
                                                                                                        tttACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAACCGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTgg                                          >  1:474558/1‑150 (MQ=255)
                                                                                                        tttACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTgg                                          <  2:431374/150‑1 (MQ=255)
                                                                                                          tACCAGCACGCTGAAAATTATTTCGCCCTTTTTACCCTGTTGTTTAATAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTAGAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAg                                        <  1:87681/150‑1 (MQ=255)
                                                                                                            ccAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGtt                                      >  2:850125/1‑150 (MQ=255)
                                                                                                              aGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGtttt                                    >  2:1690758/1‑150 (MQ=255)
                                                                                                               gCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTtta                                   >  1:649492/1‑150 (MQ=255)
                                                                                                               gCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTtta                                   <  2:1174577/150‑1 (MQ=255)
                                                                                                                cACGCTGAAAATTCTGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAACCGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTttat                                  <  1:1531594/150‑1 (MQ=255)
                                                                                                                cACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTttat                                  >  2:1595356/1‑150 (MQ=255)
                                                                                                                cACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGGCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTttat                                  <  1:320967/150‑1 (MQ=255)
                                                                                                                  cGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTttattt                                >  2:1520349/1‑149 (MQ=255)
                                                                                                                  cGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCATAAAACACCGTACCCGGCCTGGAGTTttattt                                >  2:788758/1‑149 (MQ=255)
                                                                                                                   gCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAg                               >  2:1354277/1‑150 (MQ=255)
                                                                                                                    cTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGTCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGt                              >  2:1317875/1‑150 (MQ=255)
                                                                                                                      gAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTAc                            <  2:1195322/150‑1 (MQ=255)
                                                                                                                       aaaaTTATGTCGCCCTTTTTACCCGGTGTTTTAAATGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACa                           <  1:718255/150‑1 (MQ=255)
                                                                                                                        aaaTTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTttatta                                <  2:1296868/144‑1 (MQ=255)
                                                                                                                        aaaTTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTttatta                                >  1:1296868/1‑144 (MQ=255)
                                                                                                                                tCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTAc                  <  2:1309366/150‑1 (MQ=255)
                                                                                                                                 cctgcTTTTAACCCACCGGGTAAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAACCGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTAcc                 >  2:800649/5‑150 (MQ=255)
                                                                                                                                  gCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCg                >  2:574846/1‑150 (MQ=255)
                                                                                                                                  gCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCg                >  1:1722976/1‑150 (MQ=255)
                                                                                                                                  gCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCg                <  2:78513/150‑1 (MQ=255)
                                                                                                                                  gCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTCGATGCAGACATAAGTTACTCTTCAGTAAAACGCGTTTACAAAGTCATATTGTAATGCCGGATGCGGCGTATAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCTGCTTACCg                >  1:258844/1‑150 (MQ=255)
                                                                                                                                  ctgcTTTTAACCCACCGGGTAAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTTGTACAGCGGCTTACCg                >  1:933610/4‑150 (MQ=255)
                                                                                                                                  ctgcTTTTAACCCACCGGGTAAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCg                >  1:481245/4‑150 (MQ=255)
                                                                                                                                    ccTTTTTACCCGGTGGGTTAAATGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGATAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCt              <  1:647668/150‑1 (MQ=255)
                                                                                                                                    ccTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCt              <  2:1658851/150‑1 (MQ=255)
                                                                                                                                        tttACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTg          <  2:1451289/150‑1 (MQ=255)
                                                                                                                                        tttACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTg          <  1:1652242/150‑1 (MQ=255)
                                                                                                                                            cccGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGTGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCtt      >  1:852556/1‑150 (MQ=255)
                                                                                                                                                gTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAATTTACTCTTCAGTTAAACGCGGTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAat  <  1:1001642/150‑1 (MQ=255)
                                                                                                                                                gTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAat  <  2:420403/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAAT  >  NC_000913/3910224‑3910517

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