Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,822,073 T→C 7.0% I62V (ATC→GTC) ‡ osmE ← osmotically‑inducible lipoprotein OsmE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,822,0730TC 7.0% 297.0 / 9.2 100I62V (ATC→GTC) ‡osmEosmotically‑inducible lipoprotein OsmE
Reads supporting (aligned to +/- strand):  ref base T (44/49);  new base C (5/2);  total (49/51)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

GCAGCCTGTGGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGC  >  NC_000913/1821925‑1822220
                                                                                                                                                    |                                                                                                                                                   
gCAGCCTGTGGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATc                                                                                                                                                    <  2:1232725/150‑1 (MQ=255)
    ccTGTGGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGc                                                                                                                                                <  1:250220/150‑1 (MQ=255)
    ccTGTGGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGc                                                                                                                                                <  1:1596344/150‑1 (MQ=255)
     cTGTGGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCt                                                                                                                                               >  2:1446463/1‑150 (MQ=255)
        tgGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCAc                                                                                                                                            >  1:583768/1‑150 (MQ=255)
          gATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACtt                                                                                                                                          >  1:1583192/1‑150 (MQ=255)
          gATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGATCCATGCTcgcgg                                                                                                                                          <  2:687630/150‑5 (MQ=255)
          gATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTCCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACtt                                                                                                                                          >  2:1165885/1‑150 (MQ=255)
            tCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCa                                                                                                                                        <  2:69470/150‑1 (MQ=255)
             cAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAg                                                                                                                                       <  1:753176/150‑1 (MQ=255)
             cAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAg                                                                                                                                       <  2:1616654/150‑1 (MQ=255)
             cAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGATCCATGCTcgcggtac                                                                                                                                       <  1:47495/150‑8 (MQ=255)
               gtgtCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGAc                                                                                                                                     >  2:982650/1‑150 (MQ=255)
                 gtCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGa                                                                                                                                   >  2:1564959/1‑150 (MQ=255)
                   cGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAg                                                                                                                                 <  1:368957/150‑1 (MQ=255)
                    gTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAgg                                                                                                                                <  1:350477/150‑1 (MQ=255)
                        tCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTa                                                                                                                            <  1:174520/150‑1 (MQ=255)
                           gCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTAccc                                                                                                                         >  1:1105826/1‑150 (MQ=255)
                            cacaGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTAcc                                                                                                                          >  1:1148090/1‑148 (MQ=255)
                            cacaGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTAcc                                                                                                                          <  2:1148090/148‑1 (MQ=255)
                                ggTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAAt                                                                                                                    <  1:1048194/150‑1 (MQ=255)
                                ggTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAAt                                                                                                                    >  1:164984/1‑150 (MQ=255)
                                 gTCTGATAACCGTAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAAtt                                                                                                                   <  1:1165885/150‑1 (MQ=255)
                                     gATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTtgtg                                                                                                               <  2:1583192/150‑1 (MQ=255)
                                     gATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTtgtg                                                                                                               <  2:944289/150‑1 (MQ=255)
                                     gATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTtgtg                                                                                                               <  1:1608858/150‑1 (MQ=255)
                                     gATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTtgtg                                                                                                               >  2:876082/1‑150 (MQ=255)
                                      aTAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGc                                                                                                              <  1:544255/150‑1 (MQ=255)
                                        aaCCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCaa                                                                                                            >  1:296568/1‑150 (MQ=255)
                                         aCCGGAGTTGATGACATGTCCGGTATCATCTACCGCGACAAAGAAGGTTTCTGCTTTACCATCACGTTTACCCAGGATGTAGGTCTGGCAAGTACAGCGAGCATGGATCATGCTCACATCAGACGAAGGTTTACCCGCAATTTGTGCAAc                                                                                                           >  1:1311048/1‑150 (MQ=255)
                                          ccGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTAGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCAAGCTCACATCAGACGAAGGTTTACCCGCAATTTGTGCAAcc                                                                                                          >  1:1532616/1‑150 (MQ=255)
                                           cGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTTACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCt                                                                                                         >  1:1231427/1‑150 (MQ=255)
                                            ggAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATTGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTg                                                                                                        <  2:673181/150‑1 (MQ=255)
                                             gAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTgc                                                                                                       >  1:971343/1‑150 (MQ=255)
                                               gTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTgcgc                                                                                                     >  2:117112/1‑150 (MQ=255)
                                                 tgatgaCATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGccc                                                                                                   >  2:1348038/1‑150 (MQ=255)
                                                     gaCATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCt                                                                                               <  2:997766/150‑1 (MQ=255)
                                                      aCATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTc                                                                                              >  2:1180592/1‑150 (MQ=255)
                                                      aCATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTc                                                                                              >  2:747786/1‑150 (MQ=255)
                                                       cATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTtg                                                                                                                 <  2:206946/130‑1 (MQ=255)
                                                       cATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTtg                                                                                                                 >  1:206946/1‑130 (MQ=255)
                                                         tGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATg                                                                                           >  2:578504/1‑150 (MQ=255)
                                                          gTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGc                                                                                          <  1:1555898/150‑1 (MQ=255)
                                                          gTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGc                                                                                          >  2:929105/1‑150 (MQ=255)
                                                           tCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGcc                                                                                         >  1:10329/1‑150 (MQ=255)
                                                           tCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGAAGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGcc                                                                                         >  1:777477/1‑150 (MQ=255)
                                                            ccGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCt                                                                                        >  1:1084373/1‑150 (MQ=255)
                                                            ccGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCt                                                                                        >  1:851990/1‑150 (MQ=255)
                                                               gTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCtttt                                                                                     <  2:971343/150‑1 (MQ=255)
                                                                tATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCttttt                                                                                    <  1:812061/150‑1 (MQ=255)
                                                                tATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCttttt                                                                                    <  1:299894/150‑1 (MQ=255)
                                                                   catcTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGa                                                                                 <  2:290684/150‑1 (MQ=255)
                                                                    atcTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGAc                                                                                <  1:1484603/150‑1 (MQ=255)
                                                                    atcTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCGCATGCCTTTTTTGAc                                                                                <  2:49603/150‑1 (MQ=255)
                                                                       tAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCttttt                                                                                    >  1:123557/1‑143 (MQ=255)
                                                                       tAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCttttt                                                                                    <  2:123557/143‑1 (MQ=255)
                                                                          cgcgACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCttt                                                                          >  1:1384492/1‑150 (MQ=255)
                                                                            cgACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTca                                                                        <  2:1526786/150‑1 (MQ=255)
                                                                             gACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTcac                                                                       >  2:1407676/1‑150 (MQ=255)
                                                                              aCAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTcacc                                                                      >  2:231928/1‑150 (MQ=255)
                                                                                  aGTAGGTTTCTGCTTAACCATCACGTTGACCCAGGATGTAGGTCTGGCAATTACCGCGAGGATTGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAg                                                                  <  2:1160422/150‑1 (MQ=255)
                                                                                   gTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAgg                                                                 <  1:1194448/150‑1 (MQ=255)
                                                                                    tAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGc                                                                <  1:658375/150‑1 (MQ=255)
                                                                                        tttCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTa                                                            <  1:892304/150‑1 (MQ=255)
                                                                                          tCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTTACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACa                                                          <  2:1231427/150‑1 (MQ=255)
                                                                                                ttACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTgg                                                    <  1:1564959/150‑1 (MQ=255)
                                                                                                ttACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTgg                                                    <  1:1130699/150‑1 (MQ=255)
                                                                                                ttACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTgg                                                    <  2:628609/150‑1 (MQ=255)
                                                                                                 tACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGt                                                   >  1:1094538/1‑150 (MQ=255)
                                                                                                  aCCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTc                                                  <  1:876082/150‑1 (MQ=255)
                                                                                                    cATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCtt                                                <  2:183467/150‑1 (MQ=255)
                                                                                                     aTCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCttt                                               >  1:1196346/1‑150 (MQ=255)
                                                                                                     aTCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCttt                                               <  1:801368/150‑1 (MQ=255)
                                                                                                       cACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTgg                                             <  1:756970/150‑1 (MQ=255)
                                                                                                          gTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTAc                                          <  1:929105/150‑1 (MQ=255)
                                                                                                               cccAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCa                                     <  1:231928/150‑1 (MQ=255)
                                                                                                               cccAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCa                                     >  1:1368276/1‑150 (MQ=255)
                                                                                                                ccAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCAt                                    <  1:982650/150‑1 (MQ=255)
                                                                                                                  aGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATaa                                  <  1:1446463/150‑1 (MQ=255)
                                                                                                                    gATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGc                                >  2:1315300/1‑150 (MQ=255)
                                                                                                                           gTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACaa                         >  1:901768/1‑150 (MQ=255)
                                                                                                                            tCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGACTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAAc                        >  2:1126403/1‑150 (MQ=255)
                                                                                                                             cTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAAcc                       <  2:1084373/150‑1 (MQ=255)
                                                                                                                               ggCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCg                     >  1:1608291/1‑150 (MQ=255)
                                                                                                                               ggCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCg                     >  2:1098841/1‑150 (MQ=255)
                                                                                                                               ggCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGACTTTTTTGACGTCTTTCACCACATGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCg                     <  2:1311048/150‑1 (MQ=255)
                                                                                                                                 cAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGcc                   <  2:1158851/150‑1 (MQ=255)
                                                                                                                                 cAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGcc                   >  1:476065/1‑150 (MQ=255)
                                                                                                                                    ctgaagtgAGCATGATCCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGc                >  1:1225757/8‑150 (MQ=255)
                                                                                                                                    ctgaagtgAGCATGATCCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGc                >  1:712853/8‑150 (MQ=255)
                                                                                                                                        agtgAGCATGATCCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCGTGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATgg            >  1:550087/4‑150 (MQ=255)
                                                                                                                                        agtgAGCATGATCCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATgg            >  1:69904/4‑150 (MQ=255)
                                                                                                                                         gtgAGCATGATCCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGt           >  2:897575/3‑150 (MQ=255)
                                                                                                                                         gcgAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGt           >  2:587914/1‑150 (MQ=255)
                                                                                                                                          cgAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGtt          <  1:1098841/150‑1 (MQ=255)
                                                                                                                                             gCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAAt       >  1:506269/1‑150 (MQ=255)
                                                                                                                                              cATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATa      <  1:1407676/150‑1 (MQ=255)
                                                                                                                                                 ggATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCACCCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATAccg   >  2:1569408/1‑150 (MQ=255)
                                                                                                                                                 ggATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATAccg   >  1:745262/1‑150 (MQ=255)
                                                                                                                                                  gATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATAccgc  >  1:797048/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCAGCCTGTGGATCAGTGTCGTATTCAGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGTACAACCCGCCAGCATGGTTAATACCGC  >  NC_000913/1821925‑1822220

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