Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 51,908 C→G 10.0% E175Q (GAA→CAA) ‡ 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,9080CG10.0% 222.2 / 13.9 70E175Q (GAA→CAA) ‡rsmA16S rRNA m(6)2A1518,m(6)2A1519 dimethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (28/35);  new base G (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.92e-01

TTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGAT  >  NC_000913/51762‑52057
                                                                                                                                                  |                                                                                                                                                     
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                                                                                                                                        <  2:860890/133‑1 (MQ=255)
                             tgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATcac                                                                                                                                        >  1:860890/1‑133 (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)
                                                                                                                                                  |                                                                                                                                                     
TTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGAT  >  NC_000913/51762‑52057

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: