Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,517,484 C→G 5.8% I32M (ATC→ATG ydcX → DUF2566 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,517,4840CG 5.8% 367.2 / 10.4 116I32M (ATC→ATGydcXDUF2566 family protein
Reads supporting (aligned to +/- strand):  ref base C (50/59);  new base G (2/5);  total (52/64)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

Only 100 of 120 total aligned reads displayed.

CAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTG  >  NC_000913/1517335‑1517633
                                                                                                                                                     |                                                                                                                                                     
cAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATc                                                                                                                                                       >  2:1318446/1‑150 (MQ=255)
cAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATc                                                                                                                                                       >  2:1744446/1‑150 (MQ=255)
    aaaTGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGtt                                                                                                                                                   <  2:1395607/150‑1 (MQ=255)
     aaTGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGttt                                                                                                                                                  <  2:1076712/150‑1 (MQ=255)
      aTGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACATGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTc                                                                                                                                                 <  2:467096/150‑1 (MQ=255)
       tgccGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCt                                                                                                                                                <  2:556789/147‑1 (MQ=255)
       tGACGCATATTTGCGCGCGTTTTATTCATCTGGCTTGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCt                                                                                                                                                <  1:1075818/150‑1 (MQ=255)
       tGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCt                                                                                                                                                <  2:1549467/150‑1 (MQ=255)
            cATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                                                                           <  1:1539898/150‑4 (MQ=255)
            cATATTTGCGCGCGTTTAATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCCCATTTCTAATCACCTTGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                                                                           <  2:135216/150‑4 (MQ=255)
            cATATTTGCGCGCGTGTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                                                                           <  1:1410429/150‑4 (MQ=255)
             atatTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAgcg                                                                                                                                          <  2:351454/150‑1 (MQ=255)
              tatTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAgcgc                                                                                                                                         <  2:817223/150‑1 (MQ=255)
               atTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcgtga                                                                                                                                        <  2:29339/150‑7 (MQ=255)
                        cgTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCtggtg                                                                                                                               <  1:815376/150‑1 (MQ=255)
                             tATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGc                                                                                                                          <  2:638224/150‑1 (MQ=255)
                              aTTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCa                                                                                                                         >  1:486238/1‑150 (MQ=255)
                              aTTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCa                                                                                                                         <  1:1043823/150‑1 (MQ=255)
                                 cATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATATCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCATTTCTGGTGGGGGCAACa                                                                                                                      >  2:627376/1‑150 (MQ=255)
                                     tggctggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGc                                                                                                                  >  1:567012/1‑150 (MQ=255)
                                     tggctggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGc                                                                                                                  <  1:919721/150‑1 (MQ=255)
                                        ctggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGa                                                                                                               >  2:461104/1‑150 (MQ=255)
                                        ctggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGa                                                                                                               >  1:1751049/1‑150 (MQ=255)
                                          ggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGTGCAACATGGCCGATg                                                                                                             >  2:37748/1‑150 (MQ=255)
                                          ggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATg                                                                                                             <  2:996948/150‑1 (MQ=255)
                                           gACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCTTTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGa                                                                                                            <  2:519919/150‑1 (MQ=255)
                                            aCGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAg                                                                                                           <  2:800309/150‑1 (MQ=255)
                                             cGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGt                                                                                                          >  1:1574162/1‑150 (MQ=255)
                                             cGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGt                                                                                                          >  1:848968/1‑150 (MQ=255)
                                              gcacGTACATGTCTCTCTATCAACACATGCTTGTTTTTTAAGCGGTAATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGtt                                                                                                         <  1:1269799/147‑1 (MQ=255)
                                               cccGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGttt                                                                                                        >  1:24181/1‑150 (MQ=255)
                                               cccGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGttt                                                                                                        >  1:1606961/1‑150 (MQ=255)
                                                   tACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGATTTTCCCg                                                                                                    >  2:1201207/1‑150 (MQ=255)
                                                   tACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCg                                                                                                    <  1:440166/150‑1 (MQ=255)
                                                     cATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGt                                                                                                  >  1:1215336/1‑150 (MQ=255)
                                                       tGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTgg                                                                                                <  1:1676077/150‑1 (MQ=255)
                                                        gTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGgc                                                                                               <  1:1558605/150‑1 (MQ=255)
                                                         tctctcTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGgcg                                                                                              >  2:1164951/1‑150 (MQ=255)
                                                          ctctctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGgcgc                                                                                             <  1:1517784/150‑1 (MQ=255)
                                                           tctctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCt                                                                                            >  2:1339767/1‑150 (MQ=255)
                                                            ctctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCtg                                                                                           <  2:1261313/150‑1 (MQ=255)
                                                             tctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCTCCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCtgt                                                                                          <  2:695693/150‑1 (MQ=255)
                                                             tctATCAACACATGCTTGTTTTTTATGCGGTTATCGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCCCtgt                                                                                          >  1:475720/1‑150 (MQ=255)
                                                               tATCAACCCATGCTTGTTTTTTATGCTGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCtgttg                                                                                        <  2:994918/150‑1 (MQ=255)
                                                               tATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCtgttg                                                                                        <  2:1712784/150‑1 (MQ=255)
                                                                 tCAACACATGCTTGTTTTTTATGAGGTTATGGCAGCAATCGCATTTCTTATAACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTgttgtt                                                                                      >  1:1354745/1‑150 (MQ=255)
                                                                  cAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGttt                                                                                     <  2:361796/150‑1 (MQ=255)
                                                                    acacATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGtttta                                                                                   >  2:1408238/1‑149 (MQ=255)
                                                                     cacaTGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCa                                                                                  >  2:770162/1‑150 (MQ=255)
                                                                          gCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACTCAGCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGtt                                                                             <  2:1072563/150‑1 (MQ=255)
                                                                           cTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGttt                                                                            >  1:1228933/1‑150 (MQ=255)
                                                                            ttGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGtttt                                                                           <  2:601687/150‑1 (MQ=255)
                                                                               ttttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAATGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTaaa                                                                        <  1:1175334/150‑1 (MQ=255)
                                                                               ttttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTaaa                                                                        <  2:906520/150‑1 (MQ=255)
                                                                               ttttttATGCGGTTATGGCAGCAAGCGCATTTCTTATCACCTGGTTTCTTTCTCCCGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTaaa                                                                        <  1:997804/150‑1 (MQ=255)
                                                                                   ttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCtt                                                                    >  1:935530/1‑150 (MQ=255)
                                                                                     aTGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTAc                                                                  <  1:1059184/150‑1 (MQ=255)
                                                                                     aTGCGGATATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTAc                                                                  <  1:1263771/150‑1 (MQ=255)
                                                                                            tATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTgg                                                           <  2:1617672/150‑1 (MQ=255)
                                                                                               ggCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGCTAAGAAACGCATCCGTTTCTTCAGCGCCTTTAGGGTGGGGGCAACATGGCCGATGCGTTTCCCCGTGGAGATGTTGTTTTCACTGTTGTAAATATTACTCACTGTCtt                                                        >  1:1094163/1‑150 (MQ=255)
                                                                                                 cagcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGa                                                      >  2:973221/1‑150 (MQ=255)
                                                                                                   gcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACg                                                    >  2:1088203/1‑150 (MQ=255)
                                                                                                    caATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCTACATTGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACgg                                                   <  2:843727/150‑1 (MQ=255)
                                                                                                    caATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACgg                                                   >  1:283682/1‑150 (MQ=255)
                                                                                                    caATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACgg                                                   <  1:1362161/150‑1 (MQ=255)
                                                                                                      aTCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGcc                                                 >  2:833432/1‑149 (MQ=255)
                                                                                                      aTCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGca                                                 >  1:498124/1‑150 (MQ=255)
                                                                                                        cGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCGTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGGTTTCACTGTTTTAAATCTTACTCACTGGCTTGACgccccc                                               >  2:291409/1‑145 (MQ=255)
                                                                                                        cGCATTTCTTATCACCTGGTTTCATTCTCACGATAAGAAACGCATCCTTTTCTTAATCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTATCACTGTTTTAAATCTTACTCACTGGCTTGACGGcacc                                               <  2:73721/150‑1 (MQ=255)
                                                                                                         tcATTTCTTATCACCTGGTTTCTTTCTCACGATAAGTAACGCATCCGTTTATTAAGCGCCTTTCTGGGGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGcacca                                              <  1:697625/149‑1 (MQ=255)
                                                                                                          cATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGcaccac                                             >  2:1265202/1‑150 (MQ=255)
                                                                                                               cTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGa                                        >  2:1341369/1‑150 (MQ=255)
                                                                                                                ttATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGAc                                       >  2:486782/1‑150 (MQ=255)
                                                                                                                   tCACCTGGTTGCGATCTCACGATAAGAAACGCATCCGTATCTTAAGCGCCTTGCTGGGGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTta                                    <  2:971881/150‑1 (MQ=255)
                                                                                                                    cACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTtat                                   >  2:432436/1‑150 (MQ=255)
                                                                                                                       cTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATc                                <  2:1796380/150‑1 (MQ=255)
                                                                                                                        tGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCa                               >  2:883156/1‑150 (MQ=255)
                                                                                                                           tttctttcTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCAttt                            <  1:1635954/150‑1 (MQ=255)
                                                                                                                           tttctttcTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCAttt                            <  1:27644/150‑1 (MQ=255)
                                                                                                                             tctttctCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTgg                          <  2:1536825/150‑1 (MQ=255)
                                                                                                                               tttctCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGAt                        >  1:1751444/1‑150 (MQ=255)
                                                                                                                                ttctCACGATAAGAAACGCATGCGGTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATg                       <  1:627376/150‑1 (MQ=255)
                                                                                                                                ttctCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATg                       >  1:210475/1‑150 (MQ=255)
                                                                                                                                  ctcACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGaa                     >  2:1742020/1‑150 (MQ=255)
                                                                                                                                  ctcACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTCTTTTAAATCTGACTCACTGGCTTGCCGGCACAACGCAGACTTATATCATTTGGATGaa                     >  2:704363/1‑150 (MQ=255)
                                                                                                                                   ggcgcttAAGAAACGGATGCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAAt                    >  1:422031/7‑150 (MQ=255)
                                                                                                                                    cACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATc                   >  2:1589925/1‑150 (MQ=255)
                                                                                                                                    cACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGCGGCAACATGGCAGATGAGTTTCCCGGTTGCGCTGTAGTAATCACTGTTTTAAATCATACTCACTGGCTTGACTGCACCACGCAGACTTAAATCATTTGGATGAATc                   >  1:986110/1‑150 (MQ=255)
                                                                                                                                       cttAAGAAACGGATGCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGAt                >  1:426308/3‑150 (MQ=255)
                                                                                                                                        aTAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATa               >  1:1790818/1‑150 (MQ=255)
                                                                                                                                         tAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATaa              >  1:1135791/1‑150 (MQ=255)
                                                                                                                                         tAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATaa              >  2:936641/1‑150 (MQ=255)
                                                                                                                                         tAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGAGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATaa              >  1:31833/1‑150 (MQ=255)
                                                                                                                                            taaaCGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAAtt           <  1:334366/149‑1 (MQ=255)
                                                                                                                                            gAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATAATTTGGATGAATCGATAAAtt           >  2:1704748/1‑150 (MQ=255)
                                                                                                                                              aaCGCATCCGTTTCTTATGCGCCTTTCTGGGGGGGCCAACATGGCCGATGAGTTTCCCGGTGGCGCTTTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTc         <  2:486238/150‑1 (MQ=255)
                                                                                                                                              aaCGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATTGATAAATTTc         <  2:1606961/150‑1 (MQ=255)
                                                                                                                                                cGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTcac       <  1:1318446/150‑1 (MQ=255)
                                                                                                                                                 gCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTcaca      >  1:1477092/1‑150 (MQ=255)
                                                                                                                                                     ccGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGTCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTg  >  1:241944/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTG  >  NC_000913/1517335‑1517633

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