Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,410,986 A→G 7.9% intergenic (+102/+27) dbpA → / ← ttcA ATP‑dependent RNA helicase, specific for 23S rRNA/tRNA s(2)C32 thioltransferase, iron sulfur cluster protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,410,9860AG 7.9% 213.2 / 8.5 76intergenic (+102/+27)dbpA/ttcAATP‑dependent RNA helicase, specific for 23S rRNA/tRNA s(2)C32 thioltransferase, iron sulfur cluster protein
Reads supporting (aligned to +/- strand):  ref base A (38/32);  new base G (4/2);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

CGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTT  >  NC_000913/1410839‑1411132
                                                                                                                                                   |                                                                                                                                                  
cGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACg                                                                                                                                                  <  1:388952/150‑1 (MQ=255)
 gggAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGc                                                                                                                                                 >  2:1003720/1‑150 (MQ=255)
      gATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCATGAGTAAACATCTTATTCGGGATTGCCGGATGCGCCGCTGGCAGCGTCTTATCCCGCCTCCATAAGAGTAGCCCCATACGCTTgcg                                                                                                                                            >  2:1112668/1‑150 (MQ=255)
       aTTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTgcgc                                                                                                                                           <  1:1247096/150‑1 (MQ=255)
        ttAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCa                                                                                                                                          >  2:762176/1‑150 (MQ=255)
        ttAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCa                                                                                                                                          <  1:717800/150‑1 (MQ=255)
          aaaGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATc                                                                                                                                        >  2:243150/1‑150 (MQ=255)
             ggAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCggg                                                                                                                                     <  2:632113/150‑1 (MQ=255)
              gAAAACCGTGCCGGGTGCGGGTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTCCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATGCGCAAGCGTATCGGgc                                                                                                                                    <  2:520841/150‑1 (MQ=255)
              gAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGgc                                                                                                                                    <  1:1160131/150‑1 (MQ=255)
              gAAAAACGTACCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGCTTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATGCGCAAGCGTATCGGgc                                                                                                                                    <  2:559667/150‑1 (MQ=255)
                 aaaCGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTAGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCt                                                                                                                                 <  2:626372/150‑1 (MQ=255)
                  aaCGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTACACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTa                                                                                                                                <  2:1301502/150‑1 (MQ=255)
                  aaCGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTa                                                                                                                                <  1:593975/150‑1 (MQ=255)
                     gTGCCGGGTGCGGTTATTAAAATAAAGAAATGTTCAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACTCTTGCGCATCGGGCGCTATcc                                                                                                                             <  1:953340/150‑1 (MQ=255)
                               cGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTAttta                                                                                                                   >  1:345084/1‑149 (MQ=255)
                                ggTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCa                                                                                                                  <  1:1276496/150‑1 (MQ=255)
                                 gTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCAc                                                                                                                 >  2:720204/1‑150 (MQ=255)
                                  ttattaAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACt                                                                                                                <  1:955389/150‑1 (MQ=255)
                                   tattaAAATAATGAAATGTGGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACtt                                                                                                               <  1:657533/150‑1 (MQ=255)
                                    attaAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTc                                                                                                              >  2:955238/1‑150 (MQ=255)
                                     ttaAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCCGATGCTACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGAAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCa                                                                                                             >  1:400187/1‑150 (MQ=255)
                                      taAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGAAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCaa                                                                                                            >  1:1067096/1‑150 (MQ=255)
                                             atGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACAtt                                                                                                     >  1:488340/1‑150 (MQ=255)
                                                 aaTGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGc                                                                                                 <  1:1071557/150‑1 (MQ=255)
                                                  aTGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGcc                                                                                                >  2:479503/1‑150 (MQ=255)
                                                    ggtGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCTGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACCTTCAGCCGt                                                                                              >  1:843153/3‑150 (MQ=255)
                                                    gTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCTGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGt                                                                                              >  1:843152/1‑150 (MQ=255)
                                                          ttGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCa                                                                                        <  1:1266721/150‑1 (MQ=255)
                                                              cGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCa                                                                                    <  1:152411/150‑1 (MQ=255)
                                                                ggTGCAAGGGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAAc                                                                                  >  1:1040788/1‑150 (MQ=255)
                                                                ggTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAAc                                                                                  <  1:847557/150‑1 (MQ=255)
                                                                   gCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGa                                                                               >  1:377984/1‑150 (MQ=255)
                                                                    cAAGAGTAAACATATTATTCGGGATTGCCGGATGCTACGCTGGCCGCGTCTTATCAGGCCTCAATAAGAGTAGCCCGATACGCTTGCGCATAGGTCGCTATCCTGGTTATTTCACTTCAACAACATTAAGCCGTAACTCATCCAACTGAt                                                                              >  1:349453/1‑150 (MQ=255)
                                                                      agagTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGACTCCATAAGAGTAGCCCGATACGCTTTCGCATCGGGCGCTATACTGGTTATTTCACTTCAACCACATTCAGCCGTAACTAATCCAACTGAttt                                                                            >  1:290262/1‑150 (MQ=255)
                                                                        agTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCAGCGGGCGCGATCCTGGTGATTTCTCTTCCACCACATTCAGCCGTCCCTCCTCCAAGTGATTTTc                                                                          >  1:1217057/1‑150 (MQ=255)
                                                                         gTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCa                                                                         <  2:770794/150‑1 (MQ=255)
                                                                         gTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCa                                                                         <  1:357126/150‑1 (MQ=255)
                                                                          tAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAt                                                                        <  1:283841/150‑1 (MQ=255)
                                                                            aaCATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAtct                                                                      >  2:932250/1‑150 (MQ=255)
                                                                            aaCATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCAGGTTATTTCACTTCAACAACATTCAGCCGTAACTCATCCAACTGATTTTCAgct                                                                      >  2:621971/1‑150 (MQ=255)
                                                                             aCATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATCTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAtctt                                                                     <  1:617398/150‑1 (MQ=255)
                                                                                tCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCAtcttctt                                                                  <  1:708276/150‑1 (MQ=255)
                                                                                   tATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTg                                                               <  2:506248/150‑1 (MQ=255)
                                                                                    aTTCGGTATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTgg                                                              <  2:1142783/150‑1 (MQ=255)
                                                                                    aTTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTGCATAAGAGTAGCCGGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCCCATTCGGCCGTAACTCATCCATCTGATGTTCATCTTCTTCTgg                                                              >  1:923287/1‑150 (MQ=255)
                                                                                    aTTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTga                                                              <  1:605945/150‑2 (MQ=255)
                                                                                       cGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTg                                                           >  2:1248650/1‑150 (MQ=255)
                                                                                        gggATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGc                                                          <  1:816360/150‑1 (MQ=255)
                                                                                         ggATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATTAGAGTAGCCCGATACGCATGCGCATCGGGCGCTATCGTGGTTTTGTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGcc                                                         >  1:1039430/1‑150 (MQ=255)
                                                                                           aTTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAg                                                       <  1:1248650/150‑1 (MQ=255)
                                                                                           aTTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAg                                                       >  2:869176/1‑150 (MQ=255)
                                                                                             tGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCa                                                     <  2:995107/150‑1 (MQ=255)
                                                                                                cGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTAGTTCACTTCAACCACATTCAGCCGTAACTCATCCCACTGCTTTTCATCTTCGTCTGGCTGCCCGCACGc                                                  >  2:1237281/1‑150 (MQ=255)
                                                                                                 ggATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGcc                                                 >  2:848530/1‑150 (MQ=255)
                                                                                                   aTGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCgg                                               <  1:1003720/150‑1 (MQ=255)
                                                                                                    ggCGACGCTGGCCGCGGCTTATCCGGCCTCCATAATAGTAGCCCGATACGCATGCGCAGCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCATCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGt                                              <  2:831309/149‑1 (MQ=255)
                                                                                                      cGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTg                                            >  2:444513/1‑150 (MQ=255)
                                                                                                      cGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTg                                            >  2:534634/1‑150 (MQ=255)
                                                                                                       gACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGt                                           <  2:1290193/150‑1 (MQ=255)
                                                                                                              gCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGa                                    >  1:556557/1‑150 (MQ=255)
                                                                                                                cgcgTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGAAtc                                  <  2:1036522/150‑1 (MQ=255)
                                                                                                                   gTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATctct                               >  1:765396/1‑150 (MQ=255)
                                                                                                                       tATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCGGGTTGTAGTGGGATCTCTTcgc                           >  2:531864/1‑150 (MQ=255)
                                                                                                                                   aTAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTACCTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCAGGTTGTAGTGGGATCTCTTCGCGATCAACCGCCa               >  1:548670/1‑150 (MQ=255)
                                                                                                                                      agagTAGCCTGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGAt            >  2:1307223/1‑150 (MQ=255)
                                                                                                                                          tAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCATCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCAcc        <  2:925303/150‑1 (MQ=255)
                                                                                                                                           cgCCCGATGCGCAAGCGTATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCcccc       >  2:837182/2‑150 (MQ=255)
                                                                                                                                           cgCCCGATGCGCAAGCGTATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCAccc       >  1:929626/2‑150 (MQ=255)
                                                                                                                                           cgCCCGATGCGCAAGCGTATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCAGCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGCTCTCTTCGCGATCACACGCGAGATCAccc       >  1:194358/2‑150 (MQ=255)
                                                                                                                                           aGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCAccc       >  2:970587/1‑150 (MQ=255)
                                                                                                                                             cccGATGCGCAAGCGTATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCAccccc     >  2:1025617/1‑150 (MQ=255)
                                                                                                                                             cccGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCAccccc     >  2:582498/1‑150 (MQ=255)
                                                                                                                                              ccGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCg    >  2:1007571/1‑150 (MQ=255)
                                                                                                                                              ccGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCg    >  2:519658/1‑150 (MQ=255)
                                                                                                                                                gATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGtt  >  1:726621/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGGGAAGATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTT  >  NC_000913/1410839‑1411132

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: