Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,400,986 A→C N34K (AAT→AAG)  mreB ← cell wall structural complex MreBCD, actin‑like component MreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,400,9860AC98.7% 304.6 / NA 78N34K (AAT→AAG) mreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (44/33);  minor base G (0/1);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCGCTT  >  NC_000913/3400838‑3401133
                                                                                                                                                    |                                                                                                                                                   
agTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCtt                                                                                                                                                   >  1:439667/1‑151 (MQ=255)
   cGGCGATAACGACGACATTCATTGGGCGAATGGCAGAAATAATGCACGGCGTACGACCCAGCATCTGCTTCGCGTCATGACATACTGCAGCTTCGCTTTTCGGTGAAACGGCAAGATCCTGACGAATGGCCACCACGGAAGGCTCGTTct                                                                                                                                                 <  1:525948/150‑2 (MQ=255)
    ggCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAAt                                                                                                                                               >  2:889797/1‑151 (MQ=255)
     gCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATa                                                                                                                                              <  1:778040/151‑1 (MQ=255)
      cGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATAc                                                                                                                                             >  2:407039/1‑151 (MQ=255)
       gATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACg                                                                                                                                            <  1:547343/151‑1 (MQ=255)
        aTAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACg                                                                                                                                            >  1:70882/1‑150 (MQ=255)
         tAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGAt                                                                                                                                          >  1:1040411/1‑151 (MQ=255)
               cGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTg                                                                                                                                    >  2:329869/1‑151 (MQ=255)
                gTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCtt                                                                                                                                     <  1:402330/149‑1 (MQ=255)
                   tttCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCt                                                                                                                                >  2:652395/1‑151 (MQ=255)
                    ttCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCtt                                                                                                                               >  2:241062/1‑151 (MQ=255)
                      cATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCtttt                                                                                                                             <  2:845911/151‑1 (MQ=255)
                      cATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCtttt                                                                                                                             <  1:318365/151‑1 (MQ=255)
                       aTTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCtttt                                                                                                                             >  1:275977/1‑150 (MQ=255)
                          gggCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACa                                                                                                                          <  2:423618/150‑1 (MQ=255)
                            gCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATa                                                                                                                        >  2:86048/1‑150 (MQ=255)
                             cGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATaa                                                                                                                       >  2:765807/1‑150 (MQ=255)
                              gAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAAt                                                                                                                     >  1:471503/1‑151 (MQ=255)
                              gAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAAt                                                                                                                     <  1:870795/151‑1 (MQ=255)
                               aaTGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTCCGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTACTTCAATACGATGCCTTGTCCTTTTACATAAAt                                                                                                                     >  2:22830/1‑150 (MQ=255)
                                    cagcaATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGt                                                                                                               <  2:269010/151‑1 (MQ=255)
                                     agcaATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGt                                                                                                               >  2:498632/1‑150 (MQ=255)
                                      gcaATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTAt                                                                                                             >  2:1030637/1‑151 (MQ=255)
                                      gcaATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTAt                                                                                                             >  2:582769/1‑151 (MQ=255)
                                       caATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTAtt                                                                                                            <  1:636087/151‑1 (MQ=255)
                                          tatTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGc                                                                                                         >  1:877880/1‑151 (MQ=255)
                                                ccGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTAcc                                                                                                   >  2:1091/1‑151 (MQ=255)
                                                  ggCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTAccc                                                                                                  >  1:124448/1‑150 (MQ=255)
                                                    cGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAgg                                                                                               <  1:580888/151‑1 (MQ=255)
                                                    cGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAgg                                                                                               <  2:709530/151‑1 (MQ=255)
                                                      tACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGt                                                                                              <  2:8903/150‑1 (MQ=255)
                                                      tACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTc                                                                                             >  2:188372/1‑151 (MQ=255)
                                                        cGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCaa                                                                                           <  1:790739/151‑1 (MQ=255)
                                                          gCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATg                                                                                         >  1:127590/1‑151 (MQ=255)
                                                          gCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATg                                                                                         <  2:530164/151‑1 (MQ=255)
                                                            ccAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATgg                                                                                        >  2:1139846/1‑150 (MQ=255)
                                                             cAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGAc                                                                                      >  2:1099848/1‑151 (MQ=255)
                                                              aGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACa                                                                                     <  1:1009946/151‑1 (MQ=255)
                                                                cATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAg                                                                                   <  2:309525/151‑1 (MQ=255)
                                                                   cTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTc                                                                                 >  2:1029219/1‑150 (MQ=255)
                                                                   cTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTc                                                                                 >  2:521105/1‑150 (MQ=255)
                                                                           cgTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGaaa                                                                         >  2:1051930/1‑150 (MQ=255)
                                                                             acATGACCGAATGCAGCTACGCTTTTCGGTAAAGCTGCACGAACCTGACGAATGTCCACCACGGAATACTCCTTCAGTACGATGCCTTGTCCTTTTACATAAATGAAGTGATTCGCAGTAGCAAAGTCAATGGACAATCCATTGGaaaa                                                                        >  1:509563/2‑149 (MQ=255)
                                                                              cATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACAt                                                                     >  2:312911/1‑151 (MQ=255)
                                                                               aTGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATg                                                                    >  1:101949/1‑151 (MQ=255)
                                                                                   ccTGCTGCAGCTTCGCTTTTCGGTGAACCGGCACGATCCTGGCGAAAGGGCACCCCGGAAGGCTCCTTCAATGCGATGGCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCa                                                                 <  2:161326/150‑1 (MQ=255)
                                                                                    cTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACg                                                               >  1:139375/1‑151 (MQ=255)
                                                                                         gCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAAt                                                           <  1:899697/150‑1 (MQ=255)
                                                                                          cAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAAttt                                                         >  1:399835/1‑151 (MQ=255)
                                                                                          cAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAAttt                                                         >  1:116773/1‑151 (MQ=255)
                                                                                              tACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTc                                                     >  2:673948/1‑151 (MQ=255)
                                                                                                   ttttCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACAt                                                >  1:866488/1‑151 (MQ=255)
                                                                                                    tttCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACa                                                 <  2:768907/149‑1 (MQ=255)
                                                                                                     ttCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATa                                               <  1:149781/150‑1 (MQ=255)
                                                                                                             aCCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGAt                                      >  2:941180/1‑151 (MQ=255)
                                                                                                                   aCGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCt                                <  2:621324/151‑1 (MQ=255)
                                                                                                                   aCGATCCTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCt                                <  2:270734/151‑1 (MQ=255)
                                                                                                                        ccTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaaa                            <  1:807805/150‑1 (MQ=255)
                                                                                                                        ccTGACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaa                             >  2:1144299/1‑149 (MQ=255)
                                                                                                                           gACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCt                         <  2:447175/150‑1 (MQ=255)
                                                                                                                           gACGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCt                         <  2:747178/150‑1 (MQ=255)
                                                                                                                             cGAATGGCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTggg                      <  2:933699/151‑1 (MQ=255)
                                                                                                                                  ggCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGga                 <  1:765782/151‑1 (MQ=255)
                                                                                                                                   gCCACCACGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGga                 >  2:768844/1‑150 (MQ=255)
                                                                                                                                    ccaccaCGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaa               <  1:900422/151‑1 (MQ=255)
                                                                                                                                      accacGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaa               >  2:192424/1‑149 (MQ=255)
                                                                                                                                       ccacGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaa               >  2:932008/1‑148 (MQ=255)
                                                                                                                                       ccacGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaa               >  2:1014637/1‑148 (MQ=255)
                                                                                                                                        cacGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaaga            >  2:1130650/1‑150 (MQ=255)
                                                                                                                                         acGGAAGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaaa          <  1:1130607/151‑1 (MQ=255)
                                                                                                                                             aaGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATcc      <  1:309800/151‑1 (MQ=255)
                                                                                                                                              aGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCg     <  2:477722/151‑1 (MQ=255)
                                                                                                                                              aGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCg     <  1:764926/151‑1 (MQ=255)
                                                                                                                                              aGGCTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCg     >  2:975110/1‑151 (MQ=255)
                                                                                                                                                 cTCCTTCAATACGATGCCTTGTCCTTTTCCATAAATGAGGGTATGCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAACCTGGGGCGGAAAAGAAAATCCGCt   <  2:128671/150‑1 (MQ=255)
                                                                                                                                                 cTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCGc    >  1:548292/1‑149 (MQ=255)
                                                                                                                                                 cTCCTTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCGCtt  >  2:206537/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCGCTT  >  NC_000913/3400838‑3401133

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