Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 51,909 C→A 10.0% L174L (CTG→CTT rsmA ← 16S rRNA m(6)2A1518,m(6)2A1519 dimethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091351,9090CA10.0% 226.1 / 11.3 70L174L (CTG→CTTrsmA16S rRNA m(6)2A1518,m(6)2A1519 dimethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (28/35);  new base A (4/3);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

TTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGATG  >  NC_000913/51762‑52058
                                                                                                                                                   |                                                                                                                                                     
tttACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAg                                                                                                                                                     <  2:1143502/150‑1 (MQ=255)
    cgacgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTAcc                                                                                                                                                 <  2:936181/150‑1 (MQ=255)
       cgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCggg                                                                                                                                              >  1:95362/1‑150 (MQ=255)
       cgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCggg                                                                                                                                              >  1:838019/1‑150 (MQ=255)
       cgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTGGAAGTACCgcc                                                                                                                                              <  1:517808/150‑3 (MQ=255)
        gCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTGGAAGTACCgccg                                                                                                                                             <  2:518409/150‑4 (MQ=255)
         cTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGAt                                                                                                                                            <  1:1056714/150‑1 (MQ=255)
           ggTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATca                                                                                                                                          <  2:687870/150‑1 (MQ=255)
            gTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATcac                                                                                                                                         >  1:986936/1‑150 (MQ=255)
              tAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACAt                                                                                                                                       >  2:784242/1‑150 (MQ=255)
               aaaGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACAtt                                                                                                                                      >  2:239230/1‑150 (MQ=255)
                aaGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTg                                                                                                                                     <  2:92558/150‑1 (MQ=255)
                  ggCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCa                                                                                                                                   >  2:246333/1‑150 (MQ=255)
                    cTTCGGTGGTGATGCGGCTCAACACACGAACATATTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAAt                                                                                                                                 <  2:871380/150‑1 (MQ=255)
                     ttCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATa                                                                                                                                <  2:1135009/150‑1 (MQ=255)
                        ggtggtGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTa                                                                                                                             <  1:27637/150‑1 (MQ=255)
                           ggtgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTg                                                                                                                          >  2:528258/1‑150 (MQ=255)
                             tgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATcac                                                                                                                                         >  1:860890/1‑133 (MQ=255)
                             tgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATcac                                                                                                                                         <  2:860890/133‑1 (MQ=255)
                               aTGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGcc                                                                                                                      <  2:508263/150‑1 (MQ=255)
                                tGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCa                                                                                                                     >  1:1342693/1‑150 (MQ=255)
                                 gCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCAt                                                                                                                    <  2:146615/150‑1 (MQ=255)
                                   ggCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGa                                                                                                                  >  1:901900/1‑150 (MQ=255)
                                          cacacGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTGCCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTaa                                                                                                           <  1:1090623/150‑1 (MQ=255)
                                             acGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCg                                                                                                        <  2:321635/150‑1 (MQ=255)
                                                aaCATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGAcc                                                                                                     <  1:875136/150‑1 (MQ=255)
                                                   aTCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATa                                                                                                  <  1:185468/150‑1 (MQ=255)
                                                   aTCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATa                                                                                                  <  1:246333/150‑1 (MQ=255)
                                                    tCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGa                                                                                                                  >  1:1397857/1‑133 (MQ=255)
                                                    tCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGa                                                                                                                  <  2:1397857/133‑1 (MQ=255)
                                                    tCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATAc                                                                                                 >  2:1349404/1‑150 (MQ=255)
                                                    tCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATAc                                                                                                 >  2:1000670/1‑150 (MQ=255)
                                                       ttAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCt                                                                                              <  1:1028989/150‑1 (MQ=255)
                                                          aCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTg                                                                                           >  2:1125697/1‑150 (MQ=255)
                                                           ccGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGc                                                                                          <  2:528183/150‑1 (MQ=255)
                                                             gggTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTg                                                                                        <  1:412028/150‑1 (MQ=255)
                                                                 gAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCg                                                                                    >  1:1200313/1‑150 (MQ=255)
                                                                        cGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGc                                                                                          <  2:75491/137‑1 (MQ=255)
                                                                        cGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGc                                                                                          >  1:75491/1‑137 (MQ=255)
                                                                          ttGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCaa                                                                           <  1:57535/150‑1 (MQ=255)
                                                                           tGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAAc                                                                          >  2:1001247/1‑150 (MQ=255)
                                                                            gCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCCTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAAcc                                                                         <  2:95362/150‑1 (MQ=255)
                                                                            gCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATCGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAAcc                                                                         >  1:585280/1‑150 (MQ=255)
                                                                             cATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCa                                                                        <  2:838019/150‑1 (MQ=255)
                                                                               tGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGa                                                                      >  1:1159778/1‑150 (MQ=255)
                                                                                     aCCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTc                                                                <  2:64537/150‑1 (MQ=255)
                                                                                       cAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcac                                                              >  1:951143/1‑150 (MQ=255)
                                                                                           cgcACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCaccacc                                                          <  2:1074954/150‑1 (MQ=255)
                                                                                                gacgGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCttt                                                     <  2:1200388/150‑1 (MQ=255)
                                                                                                  cggcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCttttt                                                   >  2:1082022/1‑150 (MQ=255)
                                                                                                     cggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCa                                                <  1:120107/150‑1 (MQ=255)
                                                                                                      ggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGACTCACCACCTCTTTTTGCaa                                               <  1:528258/150‑1 (MQ=255)
                                                                                                            cACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATaaa                                         >  2:997251/1‑150 (MQ=255)
                                                                                                                      ggtggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCg                               <  1:1001247/150‑1 (MQ=255)
                                                                                                                        ttgtgtAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc                             <  1:846722/148‑1 (MQ=255)
                                                                                                                            gtAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAAACGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATg                         >  1:284234/1‑150 (MQ=255)
                                                                                                                              aaaGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGc                       <  2:986936/150‑1 (MQ=255)
                                                                                                                               aaGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCa                      >  2:105954/1‑150 (MQ=255)
                                                                                                                                 ggCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATc                    >  1:1216185/1‑150 (MQ=255)
                                                                                                                                  gCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAAACGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCa                   <  2:284234/150‑1 (MQ=255)
                                                                                                                                    tGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATCGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGt                 <  2:585280/150‑1 (MQ=255)
                                                                                                                                       tcccggTACTGGAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGtata              >  2:1047436/4‑150 (MQ=255)
                                                                                                                                        cccggTACTGGAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTca                                                               >  2:1343109/3‑100 (MQ=255)
                                                                                                                                        cccggTACTGGAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTca                                                               <  1:1343109/98‑1 (MQ=255)
                                                                                                                                        cccggTACTGGAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAg             >  2:1354753/3‑150 (MQ=255)
                                                                                                                                         gcggTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGc            >  2:630533/1‑150 (MQ=255)
                                                                                                                                         ccggTACTGGAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGc            >  1:731854/2‑150 (MQ=255)
                                                                                                                                            gTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTaa         >  2:1242056/1‑150 (MQ=255)
                                                                                                                                             tACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTaaa        <  1:239230/150‑1 (MQ=255)
                                                                                                                                                ttCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCAGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAg     >  1:1060573/1‑150 (MQ=255)
                                                                                                                                                  ccAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGAt   >  1:180499/1‑150 (MQ=255)
                                                                                                                                                   cAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGATg  >  2:1138468/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGATG  >  NC_000913/51762‑52058

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: