Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 608,037 A→G 5.1% intergenic (+49/‑28) hokE → / → insL1 toxic polypeptide, small/IS186 transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913608,0370AG 5.1% 224.1 / 3.1 80intergenic (+49/‑28)hokE/insL1toxic polypeptide, small/IS186 transposase
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base G (2/2);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

CGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTC  >  NC_000913/607897‑608155
                                                                                                                                            |                                                                                                                      
cGCTTCTGGTCGGAGCCTCGGTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAAGCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                                                                                               >  2:790740/1‑150 (MQ=255)
cGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                                                                                               <  2:1097745/150‑1 (MQ=255)
cGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                                                                                               >  2:1304916/1‑150 (MQ=255)
cGCTTCTGGTCGGAGACGCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                                                                                               <  1:41771/150‑1 (MQ=255)
  cTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                                                                                             <  1:1262059/150‑1 (MQ=255)
          cGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTGTTACGCCTGATATGATTTa                                                                                                     >  2:799775/1‑150 (MQ=255)
          cGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTGTTACGCCTGATATGATTTa                                                                                                     >  2:791056/1‑150 (MQ=255)
          cGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                                                                                                     >  1:462287/1‑150 (MQ=255)
            gagaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTTAGGCTTTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc                                                                                                   <  1:498391/150‑1 (MQ=255)
             cgaCTCGCTGGGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                                                                                                  <  1:942056/149‑1 (MQ=255)
             agaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                                                                                                  <  1:761732/150‑1 (MQ=255)
             agaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCAGAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                                                                                                  <  1:363476/150‑1 (MQ=255)
              gaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCTAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTAATACGCCTGATATGATTTAACGt                                                                                                 <  2:1002904/150‑1 (MQ=255)
                cTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGc                                                                                               <  1:602300/150‑1 (MQ=255)
                 tCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGcc                                                                                              <  1:865221/150‑1 (MQ=255)
                  cGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTCGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGAATTACGCCTGATATGATTTAACGTGCCg                                                                                             >  2:768926/1‑150 (MQ=255)
                  cGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                                                             >  2:1093252/1‑150 (MQ=255)
                  cGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                                                             >  2:748806/1‑150 (MQ=255)
                  cGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGATTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTANCGTGCCg                                                                                             >  2:1087939/1‑150 (MQ=255)
                    cTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTTTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                                                           >  1:1263561/1‑150 (MQ=255)
                      gtgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa                                                                                         <  2:1304433/150‑1 (MQ=255)
                      gtgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa                                                                                         >  2:1196467/1‑150 (MQ=255)
                       tgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa                                                                                        >  2:423063/1‑150 (MQ=255)
                        gtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                                                                                       >  2:1083138/1‑150 (MQ=255)
                          gAGTTTACGGTGAAGGAACGTAATATTTAGTTTAAGGCTGTTCTCGCTTACTAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAATCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa                                                                                     <  2:61408/150‑1 (MQ=255)
                           cgattACGGTGAAGGAACGTAATATAGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAAACCCATAAGCGCTAACTTAAGGGTAGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc                                                                                    <  2:455574/147‑1 (MQ=255)
                             tttACGGTGAAGGAAGGTAATATTGAGTTTATGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc                                                                                  <  1:216494/150‑1 (MQ=255)
                               tACGGGTAATGAACTTAATATTGAGTTTAAGTCTGTTCTCGCTTACGAACCGAATATGTAGCCGTTGAGCGGTGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc                                                                                <  1:9304/150‑1 (MQ=21)
                                 cGGCGACGGAACGTAATATTGCGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGGTGTGCTGGGAGTAATCCCATAATCGCTAACTTAATGGTTGTGGTTTTACGCCTGATATGATTTAACGTGCCGTTGAATTACTCTCAc                                                                              <  2:813659/150‑1 (MQ=25)
                                   gTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATACCATAAGCGCTAACTTAAGGGTTGTGGGATTACGCCGGATATGATTTACCGTGCCGATGAATTACTCTCACGa                                                                            >  1:212224/1‑150 (MQ=255)
                                        ggAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACt                                                                       >  1:575091/1‑150 (MQ=255)
                                        ggAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACt                                                                       <  1:447509/150‑1 (MQ=255)
                                        ggAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAAATTAAGGGTTGTGGTATGACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACt                                                                       >  2:1155243/1‑150 (MQ=255)
                                        ggAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAAGCCCATAAGCGCTAACATAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACt                                                                       <  2:210906/150‑1 (MQ=255)
                                          aaCGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTgg                                                                     >  1:524179/1‑150 (MQ=255)
                                          aaCGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTgg                                                                     >  1:1129929/1‑150 (MQ=255)
                                          aaCGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTgg                                                                     <  1:562128/150‑1 (MQ=255)
                                             gTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTATCTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGACGATGAATTACTCTCACGATAACTGGTCa                                                                  >  1:135062/1‑150 (MQ=255)
                                             gTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCAATAATCGCTAACTTAAGGGTTGTTGTATTACGCATGTTATGATTTAACGTGCCGATGAATTACTCTAACGATAACTGGTct                                                                  >  1:823632/1‑149 (MQ=25)
                                              tAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                                                                 <  1:1196467/150‑1 (MQ=255)
                                                  atTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAAt                                                             >  2:272302/1‑150 (MQ=255)
                                                     gAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTcg                                                          >  1:1241838/1‑149 (MQ=255)
                                                      tgTTTAAGGCTGTTCTCGCTTCCGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAGGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATtgtg                                                         <  2:607881/149‑1 (MQ=255)
                                                       gTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGATTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg                                                        <  1:765716/150‑1 (MQ=255)
                                                         ttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGGTGTGGTATTACGCCTGATAGGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGcc                                                      <  2:328779/150‑1 (MQ=255)
                                                         ttAAGGCTGTTCTCGCTTAAGAACCGAAGAAGTAGCAGTTGTGCGGGTAGTAATCCCATAAGCGCTAACTTAAGGGTTTTGGTATTACGCCTGATATTATTTAACGTGCCGATTAATTACTCTCACGATAACTGGTCAGCAATTCTGGcc                                                      <  2:391390/150‑1 (MQ=18)
                                                          tAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTATGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGccc                                                     >  1:730732/1‑150 (MQ=255)
                                                            aGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCat                                                   >  2:1256415/1‑150 (MQ=255)
                                                             ggCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGAAATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCata                                                  >  2:988499/1‑150 (MQ=255)
                                                              gcagTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCatat                                                 >  1:488947/4‑150 (MQ=255)
                                                                tGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTg                                               >  1:1158289/1‑150 (MQ=255)
                                                                ggatctcGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTg                                               <  1:799862/147‑1 (MQ=255)
                                                                   tctcGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTGTTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTa                                            <  1:791056/150‑1 (MQ=255)
                                                                   tctcGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTGTTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTa                                            <  1:799775/150‑1 (MQ=255)
                                                                   tctcGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTa                                            >  2:839792/1‑150 (MQ=255)
                                                                   tctcGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTa                                            >  2:1049223/1‑150 (MQ=255)
                                                                      cGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTAAGCATTTCTGGCCCATATTGGTAAGc                                         >  1:841911/1‑150 (MQ=255)
                                                                        cTTACGATCCTAAGAAGTAGCGGTAGTGCGGGGAGTATTCCCATAAGCGCTTACTTTAGGGTTGTGGTATTATGCCTGATATTATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGccc                                       <  1:624745/150‑1 (MQ=14)
                                                                        cTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGccc                                       >  2:416364/1‑150 (MQ=255)
                                                                             gAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaaga                                  >  2:21653/1‑150 (MQ=255)
                                                                              aaCCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAAGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaagaa                                 >  1:996309/1‑150 (MQ=255)
                                                                                ccGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGTCCCATAATGGTAAGCCCGAAGAACt                               >  2:809180/1‑150 (MQ=255)
                                                                                ccGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACt                               <  1:460523/150‑1 (MQ=255)
                                                                                ccGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACt                               >  2:1215735/1‑150 (MQ=255)
                                                                                  gaagaaGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTgg                             >  2:450451/1‑150 (MQ=255)
                                                                                  gaagaaGTAGACGGGGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGAGGTGGTATTACGCCTGATATGATTTAACGTGCAGATGAAAAAGACTCAGGATAACTTGTCACCAATTCGAGCCCATATAGCTAAACCCGAATACCTgg                             >  2:207347/1‑150 (MQ=9)
                                                                                   aagaagTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGa                            >  2:773256/1‑150 (MQ=255)
                                                                                    agaagTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTATCGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGAt                           <  1:818729/150‑1 (MQ=255)
                                                                                    agaagTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGAt                           <  2:45740/150‑1 (MQ=255)
                                                                                       agaagacGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATCTGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACt                        <  2:1139193/144‑1 (MQ=255)
                                                                                         ttgCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGCTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATAACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATCGGTAAGCCCGAAGAACTGGATACTTc                      <  2:975535/148‑1 (MQ=255)
                                                                                            ccGTTGTGCGGGGAGTAATCGCATAAGCGCTTACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGTATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGc                   <  1:333369/150‑1 (MQ=255)
                                                                                            ccGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCCCGATAACTGGTCAGCAATTCTGGCCAATATTGGTAAGCCCGAAGAACTGGATACTTCGGc                   >  1:1127299/1‑150 (MQ=255)
                                                                                            ccGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGccc                                       >  2:697605/1‑130 (MQ=255)
                                                                                            ccGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGccc                                       <  1:697605/130‑1 (MQ=255)
                                                                                                       ggAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAAATCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCg        <  1:90801/150‑1 (MQ=25)
                                                                                                        gAGTAATCCCATAAGCGCTAACTTAAGGGTAGTGGTATTACGCCTGATATGATTTAAAGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTTGCCAATATTGGTAAGCCCGATGAACTGGATACTTCGGCACGTAATGCCgg       >  2:704281/1‑150 (MQ=11)
                                                                                                        gAGTAATCCAATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATTACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGATATGGATACTTCGGCACGTAATGCCgg       >  2:822282/1‑150 (MQ=255)
                                                                                                            aaTCCCATAAGCGCTAACTTAAGGGTTTTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGACCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGc  >  2:1107738/1‑149 (MQ=14)
                                                                                                             aTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGctc  <  1:1287461/150‑1 (MQ=11)
                                                                                                                                            |                                                                                                                      
CGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTC  >  NC_000913/607897‑608155

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