Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,410,991 T→A 7.3% intergenic (+107/+22) 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,9910TA 7.3% 212.3 / 4.1 77intergenic (+107/+22)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 T (39/32);  new base A (4/2);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACC  >  NC_000913/1410845‑1411136
                                                                                                                                                  |                                                                                                                                                 
gATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCATGAGTAAACATCTTATTCGGGATTGCCGGATGCGCCGCTGGCAGCGTCTTATCCCGCCTCCATAAGAGTAGCCCCATACGCTTgcg                                                                                                                                                >  2:1112668/1‑150 (MQ=255)
 aTTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTgcgc                                                                                                                                               <  1:1247096/150‑1 (MQ=255)
  ttAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCa                                                                                                                                              <  1:717800/150‑1 (MQ=255)
  ttAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCa                                                                                                                                              >  2:762176/1‑150 (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                                                           >  2:869176/1‑150 (MQ=255)
                                                                                     aTTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAg                                                           <  1:1248650/150‑1 (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:519658/1‑150 (MQ=255)
                                                                                                                                        ccGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCg        >  2:1007571/1‑150 (MQ=255)
                                                                                                                                          gATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGtt      >  1:726621/1‑150 (MQ=255)
                                                                                                                                              cGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAAcc  >  1:733427/1‑150 (MQ=255)
                                                                                                                                              cGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAAcc  >  1:991974/1‑150 (MQ=255)
                                                                                                                                              cGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAAcc  <  1:479503/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACC  >  NC_000913/1410845‑1411136

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: