Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,923,204 G→A G159G (GGC→GGT rsmG ← 16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,923,2040GA94.1% 278.1 / 6.6 85G159G (GGC→GGTrsmG16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (39/41);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

AAATTTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGC  >  NC_000913/3923059‑3923348
                                                                                                                                                 |                                                                                                                                                
aaattTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCtt                                                                                                                                              <  1:354289/150‑1 (MQ=255)
aaattTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTg                                                                                                                                             <  2:454175/151‑1 (MQ=255)
      tATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCAc                                                                                                                                        <  1:776945/150‑1 (MQ=255)
       aTTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCa                                                                                                                                      <  1:1189564/151‑1 (MQ=255)
        tttGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCAcc                                                                                                                                       >  1:91916/1‑149 (MQ=255)
            cTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAGGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGa                                                                                                                                   >  2:334970/1‑149 (MQ=255)
            cTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGa                                                                                                                                   >  2:1039052/1‑149 (MQ=255)
              ttAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGa                                                                                                                                <  2:409866/150‑1 (MQ=255)
               tAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGAt                                                                                                                               <  2:506298/150‑1 (MQ=255)
               tAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGGTTAACCACTGATTCGACCTGATATTCTTCGGCCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGAtg                                                                                                                              >  2:12430/1‑151 (MQ=255)
                 aTCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGAtggt                                                                                                                            <  2:495019/151‑1 (MQ=255)
                     ccaccaGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGc                                                                                                                         >  2:803560/1‑150 (MQ=255)
                     ccaccaGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGc                                                                                                                         >  1:193582/1‑150 (MQ=255)
                        ccaGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCa                                                                                                                     <  2:155549/151‑1 (MQ=255)
                           gATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCt                                                                                                                  <  1:748936/151‑1 (MQ=255)
                              gACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCAc                                                                                                               >  1:978167/1‑151 (MQ=255)
                              gACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCAc                                                                                                               >  1:863286/1‑151 (MQ=255)
                                       cATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGtt                                                                                                      <  2:1117132/151‑1 (MQ=255)
                                        aTCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTc                                                                                                     >  2:1005083/1‑151 (MQ=255)
                                         tCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTc                                                                                                     >  2:129150/1‑150 (MQ=255)
                                            aGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTGCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCag                                                                                                   >  1:354134/1‑149 (MQ=255)
                                            aGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCagag                                                                                                 <  1:743454/151‑1 (MQ=255)
                                            aGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCagag                                                                                                 <  2:107776/151‑1 (MQ=255)
                                            aGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCag                                                                                                   >  2:491120/1‑149 (MQ=255)
                                              ggCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCagagag                                                                                               <  1:206097/151‑1 (MQ=255)
                                               gCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTCCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAgg                                                                                              <  2:888097/151‑1 (MQ=255)
                                                  ggAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGc                                                                                             <  2:1152451/149‑1 (MQ=255)
                                                    aaCCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCaaaa                                                                                         <  1:695037/151‑1 (MQ=255)
                                                       cTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAgcg                                                                                      >  2:853266/1‑151 (MQ=255)
                                                         gAAGTTTAACCACTGATTCGATCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGcgcg                                                                                    >  1:73634/1‑151 (MQ=255)
                                                         gAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAgcg                                                                                      <  1:1131224/149‑1 (MQ=255)
                                                         gAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGcgcg                                                                                    <  1:436908/151‑1 (MQ=255)
                                                            gTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCt                                                                                 >  1:84560/1‑151 (MQ=255)
                                                                aaCCGCTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCAGCATCCGGCATTTGACCTTTCAGCGCGTAGAAACGGCCTTGCTCTCCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCATAAGGCCCGCTaa                                                                               <  1:816359/149‑1 (MQ=255)
                                                                aaCCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAAt                                                                              <  2:1074707/150‑1 (MQ=255)
                                                                 aCCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAAtt                                                                             >  2:1020046/1‑150 (MQ=255)
                                                                   cACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACg                                                                          <  2:1034441/151‑1 (MQ=255)
                                                                        aTTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATc                                                                     <  2:1148019/151‑1 (MQ=255)
                                                                        aTTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATc                                                                     <  1:299124/151‑1 (MQ=255)
                                                                         ttCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATc                                                                     >  1:1215814/1‑150 (MQ=255)
                                                                            gACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAAt                                                                 >  1:260171/1‑151 (MQ=255)
                                                                            gACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAAt                                                                 >  2:118/1‑151 (MQ=255)
                                                                            gACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAAt                                                                 <  1:225194/151‑1 (MQ=255)
                                                                              ccTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATgg                                                               >  2:154013/1‑151 (MQ=255)
                                                                              ccTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATgg                                                               <  1:477290/151‑1 (MQ=255)
                                                                                tGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATggcg                                                             >  1:128820/1‑151 (MQ=255)
                                                                                 gATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATggcgg                                                            <  1:167476/151‑1 (MQ=255)
                                                                                  atatTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATggcggc                                                           >  2:706655/1‑151 (MQ=255)
                                                                                      tcttcGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTg                                                       >  1:882584/1‑151 (MQ=255)
                                                                                         tcGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGa                                                      <  2:105540/149‑1 (MQ=255)
                                                                                          cGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAgg                                                   <  1:299790/151‑1 (MQ=255)
                                                                                          cGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAg                                                    <  1:432102/150‑1 (MQ=255)
                                                                                           gggCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGa                                                  <  1:912397/151‑1 (MQ=255)
                                                                                              caacaaAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAAc                                               >  2:108700/1‑151 (MQ=255)
                                                                                               aacaaAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAAct                                              <  1:129144/151‑1 (MQ=255)
                                                                                                  aaaGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACtctt                                           <  1:456302/151‑1 (MQ=255)
                                                                                                    aGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACtcttct                                         >  2:363198/1‑151 (MQ=255)
                                                                                                        aTTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTAccc                                     <  1:106828/151‑1 (MQ=255)
                                                                                                         tttCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCt                                    >  2:797834/1‑151 (MQ=255)
                                                                                                          ttCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCt                                    <  1:833880/150‑1 (MQ=255)
                                                                                                           tCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGc                                  >  2:937975/1‑151 (MQ=255)
                                                                                                           tCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGc                                  >  1:245309/1‑151 (MQ=255)
                                                                                                             aTCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGctc                                >  1:422162/1‑151 (MQ=255)
                                                                                                                 tCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGctc                                <  1:89959/147‑1 (MQ=255)
                                                                                                                 tCCGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGt                             >  1:918643/1‑150 (MQ=255)
                                                                                                                  ccGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGGAACTCTTCTACCCTGCTCTGTAc                           >  2:193569/1‑151 (MQ=255)
                                                                                                                  ccGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTAc                           >  1:546461/1‑151 (MQ=255)
                                                                                                                  ccGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTAc                           >  2:316269/1‑151 (MQ=255)
                                                                                                                   cGGCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACt                          <  2:306128/151‑1 (MQ=255)
                                                                                                                     gCATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTgg                        <  1:714270/151‑1 (MQ=255)
                                                                                                                      cATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGt                       >  1:528129/1‑151 (MQ=255)
                                                                                                                      cATTTGCCCTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGt                       <  1:837155/151‑1 (MQ=255)
                                                                                                                       aTTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGtt                      <  2:675734/151‑1 (MQ=255)
                                                                                                                            cccTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAat                  >  2:690157/1‑150 (MQ=255)
                                                                                                                             ccTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAatat                <  1:1075363/151‑1 (MQ=255)
                                                                                                                             ccTTTCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAatat                >  2:1206935/1‑151 (MQ=255)
                                                                                                                               tttCAGCGCGTCGAAACGACCTTGCTCACCAGGAAGATGGTGGCACTAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGTCACTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATtc              >  1:735328/1‑151 (MQ=255)
                                                                                                                               tttCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATtc              >  2:398693/1‑151 (MQ=255)
                                                                                                                                ttCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATtc              <  1:742443/150‑1 (MQ=255)
                                                                                                                                 tCAGCGCGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATtctc            >  2:1046135/1‑151 (MQ=255)
                                                                                                                                    gcgcGTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGATACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCa          <  1:973368/150‑1 (MQ=255)
                                                                                                                                       cgTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGttt      >  2:439554/1‑151 (MQ=255)
                                                                                                                                        gTAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGttt      >  1:884580/1‑150 (MQ=255)
                                                                                                                                         tAGAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTa     >  2:947412/1‑150 (MQ=255)
                                                                                                                                           gAAACGACCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGc  <  2:404711/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AAATTTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCAGCGCGTAGAAACGGCCTTGCTCACCAGGAAGATGGTGGCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGC  >  NC_000913/3923059‑3923348

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: