Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,349,053 A→T 8.9% I352N (ATT→AAT)  yciW ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,349,0530AT 8.9% 266.1 / 13.3 90I352N (ATT→AAT) yciWputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (46/36);  new base T (3/5);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

CTGAAACATAATTGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAATTTCATTTTCG  >  NC_000913/1348904‑1349202
                                                                                                                                                     |                                                                                                                                                     
cTGAAACATAATTGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATtaa                                                                                                                                                       >  2:271122/1‑150 (MQ=255)
   aaaCATTATTGTCTATTTTTGTGGCCATTACAGCGCCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATGCAGCCACACAGCCCGCTCCATACCAGCAGATtaataa                                                                                                                                                    <  2:1557063/150‑1 (MQ=255)
   aaaCATAATTGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCGGATtaataa                                                                                                                                                    >  2:1519924/1‑150 (MQ=255)
     aCATAATTGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATc                                                                                                                                                  <  1:1229743/150‑1 (MQ=255)
        tAATTGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGc                                                                                                                                               >  1:258523/1‑150 (MQ=255)
         aaTTGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTATAAATCTGct                                                                                                                                              <  2:1265928/150‑2 (MQ=255)
           ttGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTATTAATCTGctgg                                                                                                                                            <  1:1009723/150‑4 (MQ=255)
           ttGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTATTAATCTGctgg                                                                                                                                            <  1:44933/150‑4 (MQ=255)
            tGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTtg                                                                                                                                           >  2:880846/1‑150 (MQ=255)
             gTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTtga                                                                                                                                          >  2:1059152/1‑150 (MQ=255)
             gTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTtga                                                                                                                                          >  2:153811/1‑150 (MQ=255)
              tCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTGACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTATTAATCTGctggtat                                                                                                                                         <  2:712325/150‑7 (MQ=255)
               cTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTtgatg                                                                                                                                        >  1:1031412/1‑150 (MQ=255)
               cTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTtgatg                                                                                                                                        >  2:135169/1‑150 (MQ=255)
                      ggtggCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTc                                                                                                                                 >  1:1360057/1‑150 (MQ=255)
                       gtggCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCg                                                                                                                                >  2:530421/1‑150 (MQ=255)
                           ccATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCgtg                                                                                                                              >  2:492/1‑148 (MQ=255)
                           ccATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCgtg                                                                                                                              <  1:492/148‑1 (MQ=255)
                           ccATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCACCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGtt                                                                                                                            >  2:513359/1‑150 (MQ=255)
                             aTTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCg                                                                                                                          <  2:277671/150‑1 (MQ=255)
                             aTTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCg                                                                                                                          <  1:463592/150‑1 (MQ=255)
                             aTTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCg                                                                                                                          <  2:1031412/150‑1 (MQ=255)
                             aTTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCAGCTCCATACAAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCg                                                                                                                          >  1:368406/1‑150 (MQ=255)
                               tAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAg                                                                                                                        <  2:353696/150‑1 (MQ=255)
                                     ggCTGACAGTTTTACGCGAATCTGTCTGACGAGGAAAGGTTAATATCTCTCACCCAACGCGATTTTTAATAGATTTATCCAGCCACACAGCCCGCTCAATACCAGCAGATTAATAATCAGCGTTGATGATTTTTCGTCTTCGAGTAAAgg                                                                                                                  >  2:1157782/1‑150 (MQ=255)
                                       cTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGc                                                                                                                >  2:1613509/1‑150 (MQ=255)
                                            aGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAAc                                                                                                           <  2:1100951/150‑1 (MQ=255)
                                             gTTTTACGCGAATCTGTCTGCCGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCCGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACt                                                                                                          >  2:1399292/1‑150 (MQ=255)
                                             gTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATTATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACt                                                                                                          <  2:1040263/150‑1 (MQ=255)
                                               tttACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTgg                                                                                                        <  2:970393/150‑1 (MQ=255)
                                                ttACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTggg                                                                                                       <  2:209261/150‑1 (MQ=255)
                                                 tACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGggc                                                                                                      <  1:153811/150‑1 (MQ=255)
                                                     cgAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGgcggcg                                                                                                  >  2:686025/1‑150 (MQ=255)
                                                      gAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGgcgc                                                                                                 <  2:919735/150‑1 (MQ=255)
                                                       aaTCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCt                                                                                                >  1:1620428/1‑150 (MQ=255)
                                                        aTCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTa                                                                                               <  2:446030/150‑1 (MQ=255)
                                                           tgtctgACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAg                                                                                            <  1:880846/150‑1 (MQ=255)
                                                            gtctgACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGc                                                                                           <  2:72308/150‑1 (MQ=255)
                                                               tgACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGAt                                                                                        >  1:1045214/1‑150 (MQ=255)
                                                                gACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATc                                                                                       >  1:62941/1‑150 (MQ=255)
                                                                gACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATc                                                                                       <  2:551460/150‑1 (MQ=255)
                                                                gACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATc                                                                                       >  2:1001489/1‑150 (MQ=255)
                                                                   gcgGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCggg                                                                                    <  1:1165047/150‑1 (MQ=255)
                                                                     ggCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTg                                                                                  <  2:1045214/150‑1 (MQ=255)
                                                                                 atGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGcc                                                                      <  1:127653/150‑1 (MQ=255)
                                                                                   gTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCAt                                                                    >  1:1422810/1‑150 (MQ=255)
                                                                                     ctcACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTg                                                                  >  2:706778/1‑150 (MQ=255)
                                                                                     ctcACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTg                                                                  >  2:5671/1‑150 (MQ=255)
                                                                                         cccAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATc                                                              >  2:764145/1‑150 (MQ=255)
                                                                                         cccAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATc                                                              >  2:1539995/1‑150 (MQ=255)
                                                                                           cAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGc                                                            >  1:257969/1‑150 (MQ=255)
                                                                                            aaCGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGcc                                                           <  1:271122/150‑1 (MQ=255)
                                                                                             aCGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCt                                                          >  1:1477372/1‑150 (MQ=255)
                                                                                             aCGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCt                                                          <  2:964173/150‑1 (MQ=255)
                                                                                              cgcgATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTg                                                         >  1:525714/1‑150 (MQ=255)
                                                                                              cgcgATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTg                                                         >  2:107522/1‑150 (MQ=255)
                                                                                               gcgATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGa                                                        >  1:1431926/1‑150 (MQ=255)
                                                                                               gcgATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGa                                                        <  2:1360057/150‑1 (MQ=255)
                                                                                               gcgATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGa                                                        >  1:709199/1‑150 (MQ=255)
                                                                                                   tttttAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACgg                                                    >  1:189327/1‑150 (MQ=255)
                                                                                                   tttttAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACgg                                                    >  2:606826/1‑150 (MQ=255)
                                                                                                      ttAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCGGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAg                                                 <  1:1519924/150‑1 (MQ=255)
                                                                                                      ttAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAg                                                 <  1:1059152/150‑1 (MQ=255)
                                                                                                         aGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTg                                                                                                              >  2:589743/1‑86 (MQ=255)
                                                                                                         aGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTg                                                                                                              <  1:589743/86‑1 (MQ=255)
                                                                                                           cGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAgcgcgc                                            <  1:1332921/150‑1 (MQ=255)
                                                                                                             aTTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTc                                          >  2:846335/1‑150 (MQ=255)
                                                                                                             aTTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTc                                          <  1:5671/150‑1 (MQ=255)
                                                                                                             aTTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTc                                          <  2:1477372/150‑1 (MQ=255)
                                                                                                              tttATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCa                                         >  1:598377/1‑150 (MQ=255)
                                                                                                                tATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCAGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAAc                                       >  2:1039479/1‑150 (MQ=255)
                                                                                                                  tCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAg                                     <  2:62941/150‑1 (MQ=255)
                                                                                                                       ccACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGa                                >  2:597400/1‑150 (MQ=255)
                                                                                                                        cacacaGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGAc                               <  2:1620428/150‑1 (MQ=255)
                                                                                                                          cacaGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACgg                             <  2:257969/150‑1 (MQ=255)
                                                                                                                                ccGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGcc                       >  1:373145/1‑150 (MQ=255)
                                                                                                                                ccGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGcc                       <  2:1081326/150‑1 (MQ=255)
                                                                                                                                     ccATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTg                  <  1:513359/150‑1 (MQ=255)
                                                                                                                                      catcaacgcagATTATTAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTgg                 >  1:991741/8‑150 (MQ=255)
                                                                                                                                      cATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTgg                 >  1:364046/1‑150 (MQ=255)
                                                                                                                                      cATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTgg                 >  2:801741/1‑150 (MQ=255)
                                                                                                                                        tACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGtt               <  1:530421/150‑1 (MQ=255)
                                                                                                                                        tACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGtt               >  1:1269800/1‑150 (MQ=255)
                                                                                                                                         caacgcagATTATTAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTa              >  2:708636/5‑150 (MQ=255)
                                                                                                                                         caacgcagATTATTAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTa              >  2:1252156/5‑150 (MQ=255)
                                                                                                                                         aCCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTa              >  1:992322/1‑150 (MQ=255)
                                                                                                                                          ccAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGCTCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTaa             >  2:1023353/1‑150 (MQ=255)
                                                                                                                                          aacgcagATTATTAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTaa             <  1:708636/147‑1 (MQ=255)
                                                                                                                                                gATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAATTTCAt       >  2:1593935/1‑150 (MQ=255)
                                                                                                                                                 aTTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGaa                                   <  1:524567/121‑1 (MQ=255)
                                                                                                                                                 aTTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGaa                                   >  2:524567/1‑121 (MQ=255)
                                                                                                                                                  ttAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTCAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAATTTCAttt     <  2:869419/150‑1 (MQ=255)
                                                                                                                                                  ttAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAATTTCAttt     >  2:355913/1‑150 (MQ=255)
                                                                                                                                                    aataatCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAATTTCATTTTc   >  1:1133213/1‑150 (MQ=255)
                                                                                                                                                     ataatCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAATTTCATTTTCg  <  2:1431926/150‑1 (MQ=255)
                                                                                                                                                     ataatCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAATTTCATTTTCg  <  1:686025/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGAAACATAATTGTCTATTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATCTGCGTTGATGATTTTTCGTGTTCGAGTAAAGGGCTGAACTGGGCGGCGCTAAAGCGATCGGGTGCTCGGGTCAGCCATTGTATCGCCTGAACGGTAGCGCGCTCAACAGAATGACGGTCAGCCCATTGGTTAATTTCATTTTCG  >  NC_000913/1348904‑1349202

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