Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,909,566 G→A 6.0% A278V (GCG→GTG)  pstC ← phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,909,5660GA 6.0% 348.6 / 8.9 101A278V (GCG→GTG) pstCphosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (37/57);  new base A (3/3);  total (41/60)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 106 total aligned reads displayed.

ATAGCCATTAGCGTGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCAT  >  NC_000913/3909432‑3909705
                                                                                                                                      |                                                                                                                                           
aTAGCCATTAGCGTGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGAttccgcttcc                                                                                                                                        >  2:680585/1‑140 (MQ=255)
aTAGCCATTAGCGTGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGAttccgcttcc                                                                                                                                        >  2:1556353/1‑140 (MQ=255)
aTAGCCATTAGCGTGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGAAGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGAttccgcttcc                                                                                                                                        >  2:1054401/1‑140 (MQ=255)
      aTTAGCGTGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAAt                                                                                                                                  >  1:951612/1‑140 (MQ=255)
        tAGCGTGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAAttc                                                                                                                                >  1:1858171/1‑140 (MQ=255)
        tAGCGTGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAAttc                                                                                                                                >  1:1343013/1‑140 (MQ=255)
         aGCGTGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAAttcg                                                                                                                               >  2:1037804/1‑140 (MQ=255)
         aGCGTGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAAttcg                                                                                                                               <  2:1836649/140‑1 (MQ=255)
         aGCGTGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAAttcg                                                                                                                               <  2:1837261/140‑1 (MQ=255)
            gTGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAAttcgttc                                                                                                                            >  1:842206/1‑140 (MQ=255)
             tGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCACTTCCGCAAAttcgttcg                                                                                                                           >  1:868607/1‑140 (MQ=255)
             tGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCACTTCCGCAAAttcgttcg                                                                                                                           >  1:868574/1‑140 (MQ=255)
             tGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAAGCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCACTTCCGCAAAttcgttcg                                                                                                                           >  1:917140/1‑140 (MQ=255)
              gCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGc                                                                                                                          <  2:1142028/140‑1 (MQ=255)
              gCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGc                                                                                                                          <  1:922527/140‑1 (MQ=255)
              gCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGc                                                                                                                          <  1:1185849/140‑1 (MQ=255)
               cccccTCATTCTTAGCCAGGCGCATAATCCTAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCACCGTGCAGACCGGATTCCTCGACCGCAAATTCGTTCGcc                                                                                                                         >  2:340989/1‑140 (MQ=255)
                ccccTCTTTCTTAGCCAGGCGCTTAATCTTAAACTTCGATGCTGCGAGGACGATGAAGGTAATCAGATACAGGATCAGGCCCAGTTCCATCAGTGCGGCATCGTGCTGACCGGATTCCGCTTCCGCAAATTCGTTCGCCa                                                                                                                        <  2:232110/140‑1 (MQ=255)
                    tCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAgcgc                                                                                                                    <  1:1759175/140‑1 (MQ=255)
                      aTTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCag                                                                                                                  <  1:1820108/140‑1 (MQ=255)
                      aTTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCag                                                                                                                  <  1:1666973/140‑1 (MQ=255)
                      aTTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCag                                                                                                                  <  1:1820604/140‑1 (MQ=255)
                       ttCTTAGCCAGGCGCATAATCATATACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCaga                                                                                                                 <  2:1104403/140‑1 (MQ=255)
                       ttCTTAGCCAGGCGAATAATCATCAACTTCGATGCCGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAATTTCGTTCGCCAGCGCaga                                                                                                                 <  2:781198/140‑1 (MQ=255)
                          ttAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGt                                                                                                              >  1:207871/1‑140 (MQ=255)
                          tgaGCCAGGCGCATAATCATAAACCTCGCTGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGt                                                                                                              <  1:1114511/138‑1 (MQ=255)
                           tAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTg                                                                                                             <  1:1163634/140‑1 (MQ=255)
                           tAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTg                                                                                                             <  1:1163878/140‑1 (MQ=255)
                           tAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTg                                                                                                             <  1:383497/140‑1 (MQ=255)
                            aGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGa                                                                                                            >  1:1722589/1‑140 (MQ=255)
                            aGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGa                                                                                                            <  1:1248884/140‑1 (MQ=255)
                                 ggCGCATCATCAAAAACTTCGCTGCGGCGAGGTCGATGAAGGTAATCACAAAGAGGATCAGGCCCTGTTCCATCAGTGCGGCTACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTg                                                                                                       <  2:773440/140‑1 (MQ=255)
                                 ggCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCTGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTg                                                                                                       <  2:772235/140‑1 (MQ=255)
                                  gcgcATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGt                                                                                                      >  1:1517397/1‑140 (MQ=255)
                                  gcgcATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGt                                                                                                      >  1:865894/1‑140 (MQ=255)
                                  gcgcATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGt                                                                                                      >  2:191221/1‑140 (MQ=255)
                                     cATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTATTCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGc                                                                                                   <  1:984381/140‑1 (MQ=255)
                                     cATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGc                                                                                                   <  1:1650006/140‑1 (MQ=255)
                                     cATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGc                                                                                                   <  1:1126284/140‑1 (MQ=255)
                                     cATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGc                                                                                                   <  1:1507258/140‑1 (MQ=255)
                                     cATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGc                                                                                                   <  1:250301/140‑1 (MQ=255)
                                     cATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGc                                                                                                   <  2:1625260/140‑1 (MQ=255)
                                     cATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGc                                                                                                   <  1:957482/140‑1 (MQ=255)
                                     cATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGc                                                                                                   <  1:1317513/140‑1 (MQ=255)
                                        aaTCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCg                                                                                                <  2:1383758/140‑1 (MQ=255)
                                            aTAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCat                                                                                            <  1:1462586/140‑1 (MQ=255)
                                                aCTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATAc                                                                                        >  2:38887/1‑140 (MQ=255)
                                                  ttCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCTCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAg                                                                                      >  1:1163321/1‑140 (MQ=255)
                                                  ttCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAg                                                                                      >  1:125001/1‑140 (MQ=255)
                                                  ttCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAg                                                                                      >  1:778784/1‑140 (MQ=255)
                                                    cGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCg                                                                                    >  2:950640/1‑140 (MQ=255)
                                                    cGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCg                                                                                    >  1:505752/1‑140 (MQ=255)
                                                           gcgAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCt                                                                             >  1:1744448/1‑140 (MQ=255)
                                                                 aCGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAg                                                                       <  2:951877/140‑1 (MQ=255)
                                                                 aCGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAg                                                                       <  2:1702432/140‑1 (MQ=255)
                                                                 aCGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAg                                                                       <  2:1916697/140‑1 (MQ=255)
                                                                 aCGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAg                                                                       <  2:611037/140‑1 (MQ=255)
                                                                 aCGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAg                                                                       <  2:210441/140‑1 (MQ=255)
                                                                 aCGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAg                                                                       <  2:211191/140‑1 (MQ=255)
                                                                 aCGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTCCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAg                                                                       <  1:411187/140‑1 (MQ=255)
                                                                     tGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTgg                                                                   >  2:1338785/1‑140 (MQ=255)
                                                                      gcaGGTAATCACACACAGTATCAGGCCCAGTTCCATCAGTGCGGAAACGTGCAGACCGGATTCCGCTTCCGCAAATTAGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGt                                                                  <  2:574629/138‑1 (MQ=255)
                                                                      gAAGGTAATCACCAACAGGACCAGGCACAGTTCCATCAGAGCCGCAACGTTCAGACCTGATTCAGCTTCCGCTAATTCGTTCGCCTGATCAGATGTGATACTGTTGCCCGGCATATACAGCGAGGCGATGTCGAGCTTGt                                                                  <  2:574672/140‑1 (MQ=255)
                                                                      gAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGt                                                                  <  2:1181558/140‑1 (MQ=255)
                                                                      gAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGt                                                                  >  1:1834252/1‑140 (MQ=255)
                                                                                acaAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTAcc                                                        <  1:1054401/140‑1 (MQ=255)
                                                                                acaAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTAcc                                                        <  2:194005/140‑1 (MQ=255)
                                                                                    aCAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCgatg                                                    >  1:411223/1‑140 (MQ=255)
                                                                                        gATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATaa                                                <  2:1458267/140‑1 (MQ=255)
                                                                                        gATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATaa                                                >  2:581015/1‑140 (MQ=255)
                                                                                        gATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATaa                                                >  2:143437/1‑140 (MQ=255)
                                                                                         aTCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTACCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCGGGGAGGGGTTACCGAGGATaca                                               >  2:1207696/1‑138 (MQ=255)
                                                                                            aGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGt                                            <  1:879300/140‑1 (MQ=255)
                                                                                                  aGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGc                                      <  1:619692/140‑1 (MQ=255)
                                                                                                  aGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGc                                      >  1:1631911/1‑140 (MQ=255)
                                                                                                    ttCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCTAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCa                                    <  2:1643718/140‑1 (MQ=255)
                                                                                                    ttCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCa                                    <  2:1286715/140‑1 (MQ=255)
                                                                                                    ttCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCa                                    <  2:90246/140‑1 (MQ=255)
                                                                                                    ttCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCa                                    <  1:1166027/140‑1 (MQ=255)
                                                                                                    ttCCATCAGTGCGGCAAAGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCa                                    <  2:1598237/140‑1 (MQ=255)
                                                                                                        aTCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGt                                >  1:1579590/1‑140 (MQ=255)
                                                                                                         tCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGtt                               >  1:1824224/1‑140 (MQ=255)
                                                                                                             tGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCAc                           >  1:1091205/1‑140 (MQ=255)
                                                                                                                     cGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAgcgcg                   <  2:730103/140‑1 (MQ=255)
                                                                                                                      gTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAgcgcgc                  <  1:626592/140‑1 (MQ=255)
                                                                                                                      gTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAgcgcgc                  <  1:627046/140‑1 (MQ=255)
                                                                                                                      gTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAgcgcgc                  <  1:1461037/140‑1 (MQ=255)
                                                                                                                       tGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGcgcgcg                 >  1:9275/1‑140 (MQ=255)
                                                                                                                       tGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGcgcgcg                 >  1:75903/1‑140 (MQ=255)
                                                                                                                       tGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGcgcgcg                 >  2:481484/1‑140 (MQ=255)
                                                                                                                          cgACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcc              <  1:110008/139‑1 (MQ=255)
                                                                                                                          cgACCGGATTCCGCTACCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcc              <  1:1896670/139‑1 (MQ=255)
                                                                                                                          aGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcc              <  1:1319315/140‑1 (MQ=255)
                                                                                                                          aGACCGGATTCCACTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcc              <  2:868574/140‑1 (MQ=255)
                                                                                                                          aGACCGGATTCCACTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcc              <  2:868607/140‑1 (MQ=255)
                                                                                                                          aGACCGGATTCCACTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcc              <  2:917140/140‑1 (MQ=255)
                                                                                                                                gATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcccagccc        >  2:178218/1‑140 (MQ=255)
                                                                                                                                gATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGAGATACTGTTGCCCGGCATATACAGCGAGCCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACAGAGCGCGCGGcccagccc        >  2:178197/1‑140 (MQ=255)
                                                                                                                                      gcttccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCat  >  1:1303026/1‑140 (MQ=255)
                                                                                                                                      gcttccgcAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCat  >  1:1302970/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                           
ATAGCCATTAGCGTGCCCCCTCATTCTTAGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCAT  >  NC_000913/3909432‑3909705

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: