Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,410,995 C→T 8.0% intergenic (+111/+18) 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,9950CT 8.0% 229.8 / 6.7 75intergenic (+111/+18)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 C (38/31);  new base T (4/2);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

ATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACCACTTCAGA  >  NC_000913/1410846‑1411144
                                                                                                                                                     |                                                                                                                                                     
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                                                                   >  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:479503/150‑1 (MQ=255)
                                                                                                                                             cGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAAcc          >  1:733427/1‑150 (MQ=255)
                                                                                                                                             cGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAAcc          >  1:991974/1‑150 (MQ=255)
                                                                                                                                                  gcgcATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACCACTTc     <  1:801933/150‑1 (MQ=255)
                                                                                                                                                     cATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACCACTTCAGa  >  2:739422/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTAAAGGAAAAACGTGCCGGGTGCGGTTATTAAAATAATGAAATGTTGAATTGCCGGGTGCAAGAGTAAACATCTTATTCGGGATTGCCGGATGCGACGCTGGCCGCGTCTTATCCGGCCTCCATAAGAGTAGCCCGATACGCTTGCGCATCGGGCGCTATCCTGGTTATTTCACTTCAACCACATTCAGCCGTAACTCATCCAACTGATTTTCATCTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACCACTTCAGA  >  NC_000913/1410846‑1411144

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