Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,517,481 C→G 7.1% R31R (CGC→CGG)  ydcX → DUF2566 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,517,4810CG 7.1% 290.2 / 13.0 98R31R (CGC→CGG) ydcXDUF2566 family protein
Reads supporting (aligned to +/- strand):  ref base C (39/52);  new base G (2/5);  total (41/57)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

GTCAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCAC  >  NC_000913/1517333‑1517628
                                                                                                                                                    |                                                                                                                                                   
gTCAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCACTCTATCAACACATGCTTGTTTTTTATGCGGTTATGCCAGCAATCGCATTTCTTATCACCAGGTTTCTTTCTCACGATAAGAAACGCa                                                                                                                                                    >  1:1121035/1‑150 (MQ=255)
  cAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGCCGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATc                                                                                                                                                  >  2:1541698/1‑150 (MQ=255)
  cAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATc                                                                                                                                                  >  2:13690/1‑150 (MQ=255)
   aaaaaaTGACGCATATTTGCGCGCGTTTTATTCATCTGTCTGTACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATcc                                                                                                                                                 <  2:1635113/150‑1 (MQ=255)
    aaaaaTGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCg                                                                                                                                                >  2:1460909/1‑150 (MQ=255)
    aaaaaTGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCACATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCg                                                                                                                                                >  2:1419317/1‑150 (MQ=255)
    aaaaaTGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAAAACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCg                                                                                                                                                >  2:1460911/1‑150 (MQ=255)
        aTGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTc                                                                                                                                            <  2:1152010/150‑1 (MQ=255)
        aTGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCCTTCTCACGATAAGAAACGCATCCGTTTc                                                                                                                                            <  2:743023/150‑1 (MQ=255)
            cGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTaa                                                                                                                                        >  1:634697/1‑150 (MQ=255)
             gCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatc                                                                                                                                       <  2:427867/150‑3 (MQ=255)
             gCATATTTGAGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAg                                                                                                                                       <  2:1231147/150‑1 (MQ=255)
              cATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                                                                      <  2:706246/150‑4 (MQ=255)
              cATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTCTTatcg                                                                                                                                      <  2:1383155/150‑4 (MQ=255)
              cATATTTGCGAGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTATCTATCAACAGATGCTTGTTTTTTATGAGGTTATGGCATCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTCCTTatcg                                                                                                                                      <  2:103167/150‑4 (MQ=255)
                 ttttGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTAATGGCAGCAATCTCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGcc                                                                                                                                   <  2:262873/149‑1 (MQ=255)
                 atTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGcc                                                                                                                                   >  2:908456/1‑150 (MQ=255)
                  tttGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCt                                                                                                                                  <  1:1161328/150‑1 (MQ=255)
                      cgcgcgTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCt                                                                                                                              <  2:352805/150‑1 (MQ=255)
                      cgcgcgTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCt                                                                                                                              <  2:1183765/150‑1 (MQ=255)
                         gcgTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCtggt                                                                                                                           <  1:1013073/150‑1 (MQ=255)
                         gcgTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCtggt                                                                                                                           <  2:432842/150‑1 (MQ=255)
                           gTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATTCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCtggtgg                                                                                                                         >  1:319477/1‑150 (MQ=255)
                            ttttATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATAGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTggg                                                                                                                        <  2:489746/150‑1 (MQ=255)
                              ttATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCGTGTTTTTAATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTggggg                                                                                                                      <  1:113539/150‑1 (MQ=255)
                               tATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGc                                                                                                                     <  2:344173/150‑1 (MQ=255)
                                aTTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCa                                                                                                                    >  2:358133/1‑150 (MQ=255)
                                aTTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTCTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCAGCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCa                                                                                                                    >  2:1694233/1‑150 (MQ=255)
                                 ttCATCTGACTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCaa                                                                                                                   <  1:1686389/150‑1 (MQ=255)
                                   cATCTGGCTGGACTCCAGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCAGTTTCTTAAGCGACTTTCTGGTGGGGGCacta                                                                                                                 >  2:1569926/1‑147 (MQ=255)
                                   cATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACa                                                                                                                 <  2:525209/150‑1 (MQ=255)
                                   cATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGAGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATTACAAATGCATCCGTTTCTTAAGCGCCTATCTGGTGGGGGCACCa                                                                                                                 >  1:1084767/1‑150 (MQ=255)
                                        ggctggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGTGGCAACATGGcc                                                                                                            <  1:1293744/150‑1 (MQ=255)
                                         gctggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCg                                                                                                           <  1:1325909/150‑1 (MQ=255)
                                         gctggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTGATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTATCTGGTGGGGGCAACATGGCCg                                                                                                           >  1:167233/1‑150 (MQ=255)
                                            ggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATg                                                                                                        >  2:931106/1‑150 (MQ=255)
                                            ggACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATg                                                                                                        <  1:1687228/150‑1 (MQ=255)
                                             gACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGa                                                                                                       >  1:810183/1‑150 (MQ=255)
                                              aCGGCCGGACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCTTTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAg                                                                                                      <  2:1699430/150‑1 (MQ=255)
                                              aCGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTCTCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCGTTCTGGTGGGGGCAACATGGCCGATGAg                                                                                                      <  2:1692869/150‑1 (MQ=255)
                                              aCGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAg                                                                                                      <  2:83418/150‑1 (MQ=255)
                                                gCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGtt                                                                                                    <  1:513435/150‑1 (MQ=255)
                                                gCCCGTACATGTCTCTCTATCAACACAAGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACGTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGTTGAGtt                                                                                                    <  1:726953/150‑1 (MQ=255)
                                                 cccGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCCCCTATCTGGAGGGGCCAACATGGCCGATGAGttt                                                                                                   >  2:731956/1‑150 (MQ=255)
                                                    gTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTccc                                                                                                <  2:1232329/150‑1 (MQ=255)
                                                     tACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCg                                                                                               <  1:568597/150‑1 (MQ=255)
                                                       cATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACAATGCCGATGAGTTTCCCGGt                                                                                             >  2:1091777/1‑150 (MQ=255)
                                                          gTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGgc                                                                                          <  2:1320259/150‑1 (MQ=255)
                                                           tctctcTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTc                                                                                                                                            >  1:837665/1‑99 (MQ=255)
                                                           tctctcTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGGATGCGTTTc                                                                                                                                            <  2:837665/99‑1 (MQ=255)
                                                            ctctctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGgcgc                                                                                        >  2:647935/1‑150 (MQ=255)
                                                              ctctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTATGGTGTGGGCAACATGGACGATTAGTTTCCCGTTGTCGAtg                                                                                      >  1:180759/1‑150 (MQ=255)
                                                                ctATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCtgtt                                                                                    >  2:229106/1‑150 (MQ=255)
                                                                caaTCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCtgtt                                                                                    <  1:1331864/148‑1 (MQ=255)
                                                                     aaCACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGtttt                                                                               >  1:114599/1‑150 (MQ=255)
                                                                     aaCACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGAAGAGTTTCCAGGTGGCGCTgtagttat                                                                               >  2:809914/1‑148 (MQ=255)
                                                                         caTGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGACTTTCTGGTGGGGGCAACATGGCCGATGAGTTTACCGGTGGCGCTTTTTTTTTCACt                                                                           >  1:127428/1‑150 (MQ=255)
                                                                            gCTTGTTTTTTCTGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGAAAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGtt                                                                        <  1:1616796/150‑1 (MQ=255)
                                                                             cTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACATGGTTTCTTTCTCACGATAAGAAAAGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGATGTTGTTTTCACTGttt                                                                       >  1:302877/1‑150 (MQ=255)
                                                                              tttttttttATGCGGGTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAATAAACGCATCCGTTTCTTAAGCGCCTTTCTGGGGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGtttt                                                                      <  2:1252017/147‑1 (MQ=255)
                                                                              ttGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGtttt                                                                      >  1:462834/1‑150 (MQ=255)
                                                                                 ttttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTaaa                                                                   >  2:1177801/1‑150 (MQ=255)
                                                                                   ttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATc                                                                 <  1:1509378/150‑1 (MQ=255)
                                                                                       aTGCGGTTATGGCAGCAATCGCATTTATTATCACCTGGTTTCTTTCTCACTATAACAAACGCATCCGTTTCTTAAACGCCTTTCTGGTGGGGCCAACATGGCCGATGAGTTAACCGGTGGCCCTGTTGTTTAAAATTTTTTaaaaattaa                                                             >  2:1216154/1‑144 (MQ=255)
                                                                                        tGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACt                                                            <  2:1012213/150‑1 (MQ=255)
                                                                                              tATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTgg                                                      >  2:178619/1‑150 (MQ=255)
                                                                                               aTGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGAAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGc                                                     <  1:13690/150‑1 (MQ=255)
                                                                                                tGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGTTAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCt                                                    >  1:1540075/1‑150 (MQ=255)
                                                                                                tGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTATATCTTACTCACTGGCt                                                    >  2:865808/1‑150 (MQ=255)
                                                                                                   cagcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGa                                                 >  1:834691/1‑150 (MQ=255)
                                                                                                     gcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACg                                               >  2:235327/1‑150 (MQ=255)
                                                                                                     gcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACg                                               >  2:1414307/1‑150 (MQ=255)
                                                                                                        aTCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGca                                            >  1:119113/1‑150 (MQ=255)
                                                                                                               ttCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCa                                     <  1:908456/150‑1 (MQ=255)
                                                                                                               ttCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCa                                     <  2:1341708/150‑1 (MQ=255)
                                                                                                                tCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAg                                    >  1:1138815/1‑150 (MQ=255)
                                                                                                                  ttATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGAc                                  <  2:528534/150‑1 (MQ=255)
                                                                                                                     tCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTta                               <  2:1297890/150‑1 (MQ=255)
                                                                                                                       aCCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTtata                             <  1:196897/150‑1 (MQ=255)
                                                                                                                        ccTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTGGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGCCGGCACCACGCAGACGTatat                            >  2:1042652/1‑150 (MQ=255)
                                                                                                                          tGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCa                          <  2:810183/150‑1 (MQ=255)
                                                                                                                                 tttctCTCGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGAAACATGGCCGATGAGTTTCCCGGTGGCTCTGTTGTTTTCACTGTTTTAAATCTTTCTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGAt                   <  2:1121035/150‑1 (MQ=255)
                                                                                                                                 tttctCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCTTTTGGAt                   <  1:407245/150‑1 (MQ=255)
                                                                                                                                  ttctCACGATAAGAAACGCATCCGAGTCTTAAGCTCCTTTCTGGTGGGGGCAACATTGCCGATTAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATg                  <  2:180759/150‑1 (MQ=255)
                                                                                                                                  ttctCACGATAAGAAACGCAGCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATg                  <  1:1703684/150‑1 (MQ=255)
                                                                                                                                   tctcACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTtgatga                 <  2:1194370/150‑1 (MQ=255)
                                                                                                                                   tctcACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGa                 <  1:1099950/150‑1 (MQ=255)
                                                                                                                                    ctcACGATAAGAACCGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGaa                >  2:711602/1‑150 (MQ=255)
                                                                                                                                      cACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGTGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATTAATc              <  1:1460911/150‑1 (MQ=255)
                                                                                                                                      cACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATc              >  2:59413/1‑150 (MQ=255)
                                                                                                                                      cACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATc              <  1:1419317/150‑1 (MQ=255)
                                                                                                                                      cACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATc              <  2:426652/150‑1 (MQ=255)
                                                                                                                                      cACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATc              <  1:1460909/150‑1 (MQ=255)
                                                                                                                                      cACGATAAGAAACGCATCCGTTACTTAAGCGCCGTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTGTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATc              <  2:683614/150‑1 (MQ=255)
                                                                                                                                          aTAAGAAACGCATCCGTTGCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATa          <  1:229106/150‑1 (MQ=255)
                                                                                                                                            aaGAAACGGATGCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATaaa        >  2:885094/1‑150 (MQ=255)
                                                                                                                                                 aggCAACCGTTTCTTAAGCGCATTTCTGGGGGGGGCAACATGGCCTATTCGTTACCCGG‑GGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTcac  <  2:1491310/148‑1 (MQ=255)
                                                                                                                                                 aCGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTca   <  1:647935/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCAC  >  NC_000913/1517333‑1517628

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: ‑