Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 15,319 T→C 5.0% intergenic (+21/‑126) dnaJ → / → insL1 chaperone Hsp40, DnaK co‑chaperone/IS186 transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091315,3190TC 5.0% 284.2 / 4.4 100intergenic (+21/‑126)dnaJ/insL1chaperone Hsp40, DnaK co‑chaperone/IS186 transposase
Reads supporting (aligned to +/- strand):  ref base T (49/45);  new base C (3/2);  total (52/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

AAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAAC  >  NC_000913/15172‑15465
                                                                                                                                                   |                                                                                                                                                  
aaGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGCCCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTgg                                                                                                                                                  >  2:1033443/1‑150 (MQ=255)
 aGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCACGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGc                                                                                                                                                <  2:648446/150‑1 (MQ=255)
  ggCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCCCAACAGCCCGCGCTCAAAGCGCTGCTTTGAAGGTGTGAAGAAGTTTTTTGCAGACCTGACCAGATAACCTCCCCTAAAGCCTGCCCCTGGGc                                                                                                                                                >  2:339347/1‑150 (MQ=255)
      gAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGc                                                                                                                                            >  2:779463/1‑150 (MQ=255)
      gAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGATCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTTACCTCCACAAACGCCTGCCCGTGGGCAGGc                                                                                                                                            >  1:1093483/1‑150 (MQ=255)
       acaccgctgctGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCACCCGGCGAGCACATCAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAACAGCCTGCCCGTGGGCAGGcc                                                                                                                                           <  2:355807/145‑1 (MQ=255)
         accgctgctGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGGGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTg                                                                                                                                         <  2:1836285/147‑1 (MQ=255)
           aGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCACGGGCAGGcttttg                                                                                                                                       <  2:1629460/150‑6 (MQ=255)
            gctgctGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTCACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTgggg                                                                                                                                      >  1:1727975/1‑149 (MQ=255)
            gctgctGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGt                                                                                                                                      >  2:519520/1‑150 (MQ=255)
            gctgctGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCACGGGCAGGcttttgg                                                                                                                                      <  2:230071/150‑7 (MQ=255)
                 ttctagagctgcaagaAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGTTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTaaaaa                                                                                                                                 <  2:513512/146‑1 (MQ=255)
                 tgcaagagctgcaagaAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTaaaaa                                                                                                                                 >  2:443229/1‑150 (MQ=255)
                     agagctgcaagaAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAgg                                                                                                                             <  1:361719/150‑1 (MQ=255)
                     agagctgcaagaAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAgg                                                                                                                             <  1:795443/150‑1 (MQ=255)
                      gagctgcaagaAAGCTTCGGGGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAggg                                                                                                                            <  1:1595726/150‑1 (MQ=255)
                           gcaagaAAGCTTCGGTGGCCCAACCGGCGAGCACCACAGCCCGCGCTCAAAGAGCGTCGTTGATGGTGTGAAGAAGTTTTGTGCCGACCGGCCCCGCTAACCTCCCCAACCGCCTGCCCGTGGGCAGGCCTGGGTAACAATAGGGTGCGt                                                                                                                       >  1:106121/1‑150 (MQ=255)
                           gcaagaAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGCCCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGt                                                                                                                       >  2:1980772/1‑150 (MQ=255)
                           gcaagaAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGt                                                                                                                       >  1:1172737/1‑150 (MQ=255)
                               gaAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGaa                                                                                                                   <  1:1552941/150‑1 (MQ=255)
                                caaGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAg                                                                                                                  <  2:278852/149‑1 (MQ=255)
                                 aaGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAATAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGa                                                                                                                 <  1:1750767/150‑1 (MQ=255)
                                 aaGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGa                                                                                                                 <  2:1221939/150‑1 (MQ=255)
                                   gCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGCCGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCGTGGGTAAAACTAGGGTGCGTTGAAGata                                                                                                               >  1:111270/1‑150 (MQ=255)
                                    cTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGatat                                                                                                              >  1:1327975/1‑150 (MQ=255)
                                    cTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGatat                                                                                                              >  1:854074/1‑150 (MQ=255)
                                      tCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGc                                                                                                            <  1:1160115/150‑1 (MQ=255)
                                        ggTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGGGCGTTGAAGATATGCgc                                                                                                          >  2:1004104/1‑149 (MQ=255)
                                         gTGGCCCAACCGGCGAGCACAACAGCCCCCGCTCAAAGAGCTGCTTTGATGGTGTGAAGACGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAg                                                                                                         <  1:110444/150‑1 (MQ=255)
                                               cAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTg                                                                                                   >  2:1087878/1‑150 (MQ=255)
                                                aaCCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCTGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGt                                                                                                  <  2:1635442/150‑1 (MQ=255)
                                                   cGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGAAGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTaaa                                                                                               <  1:560524/150‑1 (MQ=255)
                                                      cGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTg                                                                                            >  2:1528120/1‑150 (MQ=255)
                                                         gCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCGGTAAAGTGGCg                                                                                         >  1:1013792/1‑150 (MQ=255)
                                                          cacaACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCgg                                                                                        <  1:585937/150‑1 (MQ=255)
                                                           acaacaGCCCGCTCTCAAAGAGCTTCTTTGATGGGGTGAAGAAGTGTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGGTGAAGATATGCGAGCACCTGTAAAGTGGCggg                                                                                       <  1:1289332/150‑1 (MQ=255)
                                                           acaacaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTATAAATAGGGTGCGTTGAAGATATGAGAGCACCTGTAAAGTGGCggg                                                                                       >  1:1303932/1‑150 (MQ=255)
                                                           acaacaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCggg                                                                                       <  1:645869/150‑1 (MQ=255)
                                                            caacaGCCCGCGCTCAAAGCGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCgggg                                                                                      <  1:1633328/150‑1 (MQ=255)
                                                            caacaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAACAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCgggg                                                                                      <  2:794361/150‑1 (MQ=255)
                                                            caacaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAATAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCgggg                                                                                      >  1:51556/1‑150 (MQ=255)
                                                            caacaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATATGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCgggg                                                                                      <  2:1866550/150‑1 (MQ=255)
                                                            caacaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCgggg                                                                                      >  2:1112656/1‑150 (MQ=255)
                                                            caacaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCgggg                                                                                      <  1:1624944/150‑1 (MQ=255)
                                                               caGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCACGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTTAAGTGGCGGGGATc                                                                                   >  2:312645/1‑150 (MQ=255)
                                                                 gcccgcccTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGTGCAGGCCTGGGTAAAAATAGGGGGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCAc                                                                                 <  2:159707/150‑1 (MQ=255)
                                                                   ccGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTc                                                                               >  2:64394/1‑150 (MQ=255)
                                                                   ccGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTc                                                                               >  2:1221373/1‑150 (MQ=255)
                                                                      cgcTCAAAGAGCGTCTTTGATGGTGGGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCa                                                                            <  1:818623/150‑1 (MQ=255)
                                                                        cTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATa                                                                          >  1:1160233/1‑150 (MQ=255)
                                                                         tCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATaa                                                                         <  2:843522/150‑1 (MQ=255)
                                                                           aaaGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGCGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAgc                                                                       <  1:760380/150‑1 (MQ=255)
                                                                                gCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTaa                                                                  >  2:1061121/1‑150 (MQ=255)
                                                                                gCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCAGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGACTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAACTGGCGGGGATCACTCCCATAAGCGCTaa                                                                  >  2:1757442/1‑150 (MQ=255)
                                                                                     tttGATGGTGTGAAGAAGTTTTTTGACGACCTCACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAACATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTaa                                                             <  1:453208/150‑1 (MQ=255)
                                                                                         aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCGGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGt                                                         >  2:1011868/1‑150 (MQ=255)
                                                                                           ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtg                                                       <  2:223569/150‑1 (MQ=255)
                                                                                            gtgtGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtgt                                                      >  2:649708/1‑150 (MQ=255)
                                                                                             ggtgAAGAAGTCTTATGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtgtg                                                     <  2:771683/149‑1 (MQ=255)
                                                                                               tgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCATGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGTGAACACTCCCATAAGCGATAACTTAAGGGTTGTGGt                                                   >  2:690678/1‑150 (MQ=255)
                                                                                                gaagcaGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAACAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAGGGGTTGTGGTa                                                  <  1:274647/150‑1 (MQ=255)
                                                                                                gaagaaGTTTTTTGACGACCTGCCCCGCTAACCTCCCCAAACGCCTGGCCGTGGGGAGGCCTGGGTTAGATTTGGGTGCGTGGAAGTTTTGCGGGCACCTGGAAAGTGGCGGGGTTCACTCCCTTTAGCGCTATCGTAAGGGGTGTgggg                                                  >  2:1330286/1‑148 (MQ=255)
                                                                                                gaagaaGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCACATAAGCGCTAACTTAAGGGTTGTGGTa                                                  >  2:1248681/1‑150 (MQ=255)
                                                                                                  agaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                                >  2:1110732/1‑150 (MQ=255)
                                                                                                      gTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                            <  2:2015164/150‑1 (MQ=255)
                                                                                                        tttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                          >  1:426503/1‑150 (MQ=255)
                                                                                                        tttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                          <  1:779463/150‑1 (MQ=255)
                                                                                                         ttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGTTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                         <  2:527982/150‑1 (MQ=255)
                                                                                                         ttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                         <  1:1371339/150‑1 (MQ=255)
                                                                                                          tttGACGACCGGCCCCGCTACCCTCCCCAATAGCCGGCCCGGGGGCAGGCCTGGGTAAAAATAGGGTGCGGTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                                        <  1:164767/150‑1 (MQ=38)
                                                                                                                gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                                  >  2:408355/1‑150 (MQ=255)
                                                                                                                 acCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                                 >  1:1661729/1‑150 (MQ=255)
                                                                                                                  ccTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                                >  1:1615515/1‑150 (MQ=255)
                                                                                                                     gACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc                             >  2:995158/1‑150 (MQ=255)
                                                                                                                      aCCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                            <  1:519520/150‑1 (MQ=255)
                                                                                                                       cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                           >  2:490016/1‑150 (MQ=255)
                                                                                                                        ccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                          >  1:1670631/1‑150 (MQ=255)
                                                                                                                         cGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGGAAAAATCGGGGGCGTGGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGc                         <  2:1727975/150‑1 (MQ=255)
                                                                                                                          gCTAAGCTCCCCACAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGGTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGcc                        <  2:297942/150‑1 (MQ=255)
                                                                                                                            tAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa                      <  1:443229/150‑1 (MQ=255)
                                                                                                                            tAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATAAGCGAGCACCTGTAAAGTGGCGGGGAGCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCGGa                      >  1:1446463/1‑150 (MQ=255)
                                                                                                                             aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                     >  2:1499098/1‑150 (MQ=255)
                                                                                                                              aCCTCCCCAACAGCCTGCCCGTGGGCAGGCCTGTGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAATGTGGCGTGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg                    <  2:992630/150‑1 (MQ=255)
                                                                                                                                cTCCCCAACAGCCTGCCCGTGGGCAGGCATGGGTATAAATAGGGTGCGTTGAAGATATGCGAGAACATGGAAAGTGGCGGGGATCACTCCCATAAGCGGTAACTTAAGGGTTGTGGTAGTACGCGTGATTTGATTTAACGTGCCGATGat                  >  2:370849/1‑149 (MQ=255)
                                                                                                                                 tCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                 <  2:1922330/150‑1 (MQ=255)
                                                                                                                                 tCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                 <  1:1997106/150‑1 (MQ=255)
                                                                                                                                  ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt                <  1:1033443/150‑1 (MQ=255)
                                                                                                                                      ccaggCCTGCCCACGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGATTTActc            >  1:1539687/5‑150 (MQ=255)
                                                                                                                                      aaaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc            <  1:64394/150‑1 (MQ=255)
                                                                                                                                         ggCCTGCCCACGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCCCTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGGCTGATATGATTTAACGTGCCGCTGAATTACTCTCa         >  1:357498/2‑150 (MQ=25)
                                                                                                                                         aGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa         >  1:404696/1‑150 (MQ=255)
                                                                                                                                          gCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAc        >  2:43898/1‑150 (MQ=255)
                                                                                                                                          gCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAc        >  1:1424450/1‑150 (MQ=255)
                                                                                                                                           ccTGCCCACGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCTAGCACCTGTAACGTGGCGGGGAGCACTCCCATAAGCGATAACTTAAGGGTTGTGGTATTAAGCCTGTGAGGATTTAACGGGCCGAGGAATTACTCTCACg       >  2:77985/1‑150 (MQ=21)
                                                                                                                                            cTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGa      >  2:1894505/1‑150 (MQ=255)
                                                                                                                                             tGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGAt     >  2:984180/1‑150 (MQ=255)
                                                                                                                                              gCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATa    <  2:1103282/150‑1 (MQ=255)
                                                                                                                                              gCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCAATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATa    >  1:1808341/1‑150 (MQ=255)
                                                                                                                                               cccGTGGGCAGGCCTTGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATaa   <  1:814837/150‑1 (MQ=255)
                                                                                                                                                ccGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAAc  >  2:121097/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAAC  >  NC_000913/15172‑15465

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: