Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,517,484 C→G 7.4% 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 7.4% 279.5 / 12.8 95I32M (ATC→ATG) ydcXDUF2566 family protein
Reads supporting (aligned to +/- strand):  ref base C (36/52);  new base G (2/5);  total (38/57)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-02

CAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTG  >  NC_000913/1517335‑1517633
                                                                                                                                                     |                                                                                                                                                     
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)
                                                                                                                                                   aTCCGTTTCTTAAGCGCCTTTCTGGGGGGGGCAACATGGCCGATGAGTTTCCCGGGGGCGCGGTTGTTGTCACTGTTTTAAATATTACTCACTGGCTCGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAg    <  2:202145/150‑1 (MQ=255)
                                                                                                                                                     ccGGTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTg  <  2:114599/150‑1 (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 < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑