Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 993,258 G→T 10.3% intergenic (+24/+19) pepN → / ← ssuB aminopeptidase N/aliphatic sulfonate ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,2580GT10.3% 268.4 / 16.0 87intergenic (+24/+19)pepN/ssuBaminopeptidase N/aliphatic sulfonate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (36/42);  new base T (5/4);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

CTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCAAATCCA  >  NC_000913/993109‑993407
                                                                                                                                                      |                                                                                                                                                     
cTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAg                                                                                                                                                       <  1:147384/150‑1 (MQ=255)
    ttCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGcgcc                                                                                                                                                   <  2:1542113/150‑1 (MQ=255)
      cGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGcgccgc                                                                                                                                                 >  1:574667/1‑150 (MQ=255)
           gAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTc                                                                                                                                            <  2:1966532/150‑1 (MQ=255)
           gAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTc                                                                                                                                            <  1:1367251/150‑1 (MQ=255)
              aCGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                                                                         <  2:731192/150‑1 (MQ=255)
              aCGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTGGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTCggt                                                                                                                                         <  2:1790956/150‑2 (MQ=255)
                    cGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATTCCCCGAATGGCGGCGCTATTGCCGCCATTCggttattt                                                                                                                                    <  1:1858304/150‑7 (MQ=37)
                    cGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATTCCCCGAATGGCGGCGCTATTGCCGCCATTCggttattt                                                                                                                                    <  1:1481537/150‑7 (MQ=37)
                     gATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTa                                                                                                                                  <  2:1656390/150‑1 (MQ=255)
                         ccAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAcccc                                                                                                                              <  2:644298/150‑1 (MQ=255)
                          cAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCt                                                                                                                             >  1:1342266/1‑150 (MQ=255)
                             aCGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGtt                                                                                                                          <  1:76125/150‑1 (MQ=255)
                                  aGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGTGGCAATAGCGCCGCCATTCGGGGAATTTACCACTGTTTtctc                                                                                                                     >  1:1698201/1‑150 (MQ=255)
                                   ggAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCa                                                                                                                    >  1:1236469/1‑150 (MQ=255)
                                   ggAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCa                                                                                                                    >  1:949218/1‑150 (MQ=255)
                                    gagaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGTGAATTTAATCCTGTTTTCTCAg                                                                                                                   <  2:901702/150‑1 (MQ=255)
                                     agaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAgg                                                                                                                  <  1:1518604/150‑1 (MQ=255)
                                        aaaTGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGa                                                                                                               >  1:1789174/1‑150 (MQ=255)
                                          aTGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGt                                                                                                             <  2:689805/150‑1 (MQ=255)
                                             cgcgcgGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTc                                                                                                          >  2:1421127/1‑150 (MQ=255)
                                                 cggcgCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGAt                                                                                                      <  2:398725/150‑1 (MQ=255)
                                                     gcTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcac                                                                                                  >  2:1096546/1‑150 (MQ=255)
                                                      cTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGCTCTGTACCACAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGGTTTCTCAGGCGAGTTTCAGATTcacc                                                                                                 <  2:766245/150‑1 (MQ=255)
                                                      cTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcacc                                                                                                 <  1:1432316/150‑1 (MQ=255)
                                                      cTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcacc                                                                                                 >  1:3159/1‑150 (MQ=255)
                                                            cAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCAGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTg                                                                                           >  1:1435239/1‑150 (MQ=255)
                                                               ttGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCa                                                                                                         >  1:1611775/1‑133 (MQ=255)
                                                               ttGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCa                                                                                                         <  2:1611775/133‑1 (MQ=255)
                                                                tGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATc                                                                                       <  1:1415991/150‑1 (MQ=255)
                                                                  aaaGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAc                                                                                     >  2:1957065/1‑150 (MQ=255)
                                                                   aaGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAcc                                                                                    >  2:1806813/1‑150 (MQ=255)
                                                                   aaGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAcc                                                                                    >  1:1614183/1‑150 (MQ=255)
                                                                    tgggCTGGATTATCTCTCTGGCGATCTGTAGGTGAAGATAACTAAAGCACTGTCTTGATAATTAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAccc                                                                                   <  1:1833956/149‑1 (MQ=255)
                                                                    aGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAccc                                                                                   >  1:1529386/1‑150 (MQ=255)
                                                                        cTGGAAAATCTCTCTGGCGATCTGTACGATAAGATAACTAAAGCACTGTCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACGACGTGGCATCACCCGTTg                                                                               <  2:906505/150‑1 (MQ=255)
                                                                         tGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGc                                                                              <  1:964461/150‑1 (MQ=255)
                                                                         tGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGc                                                                              <  1:1571465/150‑1 (MQ=255)
                                                                          ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCa                                                                             >  1:539078/1‑150 (MQ=255)
                                                                          ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCa                                                                             <  1:1062616/150‑1 (MQ=255)
                                                                            aaaaTCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAac                                                                           >  1:1690024/1‑150 (MQ=255)
                                                                            aaaaTCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAac                                                                           >  1:1179850/1‑150 (MQ=255)
                                                                             aaaTCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAaca                                                                          <  1:497131/150‑1 (MQ=255)
                                                                                  tctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTcc                                                                     <  1:1743611/150‑1 (MQ=255)
                                                                                   ctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCg                                                                    <  2:1767165/150‑1 (MQ=255)
                                                                                     ctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCt                                                                  >  1:972748/1‑150 (MQ=255)
                                                                                      tGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCtt                                                                 <  2:1434728/150‑1 (MQ=255)
                                                                                        gCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTcc                                                               <  1:1957065/150‑1 (MQ=255)
                                                                                             cTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTc                                                          >  2:1673284/1‑150 (MQ=255)
                                                                                               gTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCg                                                        <  2:1342266/150‑1 (MQ=255)
                                                                                                 aCGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCt                                                      >  1:619061/1‑150 (MQ=255)
                                                                                                 aCGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCt                                                      <  2:1567887/150‑1 (MQ=255)
                                                                                                 aCGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCt                                                      <  2:1591166/150‑1 (MQ=255)
                                                                                                  cGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTa                                                     <  2:1624297/150‑1 (MQ=255)
                                                                                                    agaagaTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAAt                                                   <  1:497210/150‑1 (MQ=255)
                                                                                                       agaTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGc                                                >  1:311142/1‑150 (MQ=255)
                                                                                                          tAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCAcc                                             >  2:823412/1‑150 (MQ=255)
                                                                                                           aaCTAAAGCACTGGCATGATAAATAACCGAATGGCGGCAATAGCGCCGCCATACGGGGAATTTACCCCGGTATTCTCAGGCGAGTTTCAGATTCACCACGTTGCCTCACCCGTTGCAACCCCTCCTCTTCCAGCTCCGCTAATAGCACCg                                            >  2:1658716/1‑150 (MQ=255)
                                                                                                             cTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAc                                          >  2:1104497/1‑150 (MQ=255)
                                                                                                             atAAAGCACTTGCTTGATAAATAACCTAATGGTGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAc                                          <  2:1698201/149‑1 (MQ=255)
                                                                                                              tAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAcc                                         <  2:98054/150‑1 (MQ=255)
                                                                                                                  gCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTaa                                     >  1:1193412/1‑150 (MQ=255)
                                                                                                                   cACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAAt                                    >  2:1486839/1‑150 (MQ=255)
                                                                                                                     cTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCg                                  >  1:1606872/1‑150 (MQ=255)
                                                                                                                       ggCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAc                                >  1:739365/1‑150 (MQ=255)
                                                                                                                           tGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCgg                            <  1:1096546/150‑1 (MQ=255)
                                                                                                                           tGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCgg                            >  1:446102/1‑150 (MQ=255)
                                                                                                                             ataaataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGAc                          >  2:1299396/1‑150 (MQ=255)
                                                                                                                                   tccccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACggggg                     >  1:266927/4‑150 (MQ=255)
                                                                                                                                   taaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACggggg                     <  2:949218/150‑1 (MQ=255)
                                                                                                                                   taaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACggggg                     >  2:10701/1‑150 (MQ=255)
                                                                                                                                    ccccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGa                    >  2:1550370/3‑150 (MQ=255)
                                                                                                                                    aaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCAGCGACCCTAATCGACGCGGACGGGGGa                    <  2:539078/150‑1 (MQ=255)
                                                                                                                                      ccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGata                  >  2:922161/1‑150 (MQ=255)
                                                                                                                                      ccGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGata                  <  2:619061/150‑1 (MQ=255)
                                                                                                                                       cGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTcc                                                               <  2:75041/104‑1 (MQ=255)
                                                                                                                                       cGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTcc                                                               >  1:75041/1‑104 (MQ=255)
                                                                                                                                        gAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATc                >  2:1072455/1‑150 (MQ=255)
                                                                                                                                        gAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATc                <  2:1789174/150‑1 (MQ=255)
                                                                                                                                         aaTGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATcc               >  1:46999/1‑150 (MQ=255)
                                                                                                                                         aaTGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATcc               <  1:1659170/150‑1 (MQ=255)
                                                                                                                                         aaTGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATcc               <  2:1236469/150‑1 (MQ=255)
                                                                                                                                            ggcggcGATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACt            >  1:1743146/1‑150 (MQ=255)
                                                                                                                                             gcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTg           >  1:992641/1‑150 (MQ=255)
                                                                                                                                             gcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCATAATCGACGCGGACGGGGGATATCCACTg           <  2:1179850/150‑1 (MQ=255)
                                                                                                                                             gcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGATAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCGACTg           >  2:356657/1‑150 (MQ=255)
                                                                                                                                                   aTAGCGCCTCCATTCGGTGAATTTACCCCTTTTTTCTCAGGCGAGTTTCAGATTCACCACGTTTCATCACCCGTTGCAACACCTCATCTTCCAGCTCCGCTAATCGCACCTACCCTAATCGACGCTGACGGGTGATATCCACTTTCAAAt     <  2:137571/150‑1 (MQ=255)
                                                                                                                                                     aGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCAAATcc   <  2:192183/150‑1 (MQ=255)
                                                                                                                                                      gcgcCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCAAATCCa  <  1:823412/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCAAATCCA  >  NC_000913/993109‑993407

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