Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,195,078 (G)6→7 intergenic (+253/‑63) ibsD → / → ibsE toxic membrane protein/toxic membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,195,0721.G94.3% 239.4 / 7.0 70intergenic (+247/‑69)ibsD/ibsEtoxic membrane protein/toxic membrane protein
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base G (31/35);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACTCATCAGAGGGGG  >  NC_000913/3194922‑3195217
                                                                                                                                                       |                                                                                                                                                    
cccTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATc                                                                                                                                                       >  2:50091/1‑151 (MQ=255)
    ggCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCgggg                                                                                                                                                  >  1:621059/1‑151 (MQ=255)
    ggCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCggg                                                                                                                                                   <  2:830469/150‑1 (MQ=255)
    ggCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATCCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCgg                                                                                                                                                    <  1:986489/149‑1 (MQ=255)
          aaaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCgggggggc                                                                                                                                             >  2:630660/1‑147 (MQ=255)
            aaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGgcagca                                                                                                                                         <  1:593613/151‑1 (MQ=255)
               ccGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAgg                                                                                                                                      <  1:969221/151‑1 (MQ=255)
                 ggggATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGGTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGTCAGCAAggag                                                                                                                                    <  1:1016798/151‑3 (MQ=255)
                  gggATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGc                                                                                                                                   >  1:1213523/1‑151 (MQ=255)
                    gATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTg                                                                                                                                 >  2:231114/1‑151 (MQ=255)
                       tttttATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGaa                                                                                                                               <  2:1191996/150‑1 (MQ=255)
                            atatCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACggg                                                                                                                          >  1:963184/1‑150 (MQ=255)
                             tatCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGa                                                                                                                         >  2:507928/1‑150 (MQ=255)
                              atCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGaa                                                                                                                        >  2:980466/1‑150 (MQ=255)
                                 tGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGACAGcc                                                                                                                    >  1:186920/1‑151 (MQ=255)
                                  gCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGc                                                                                                                     <  1:263549/149‑1 (MQ=255)
                                   cGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGcccc                                                                                                                  >  2:831938/1‑151 (MQ=255)
                                    gTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCt                                                                                                                 >  2:384625/1‑151 (MQ=255)
                                        cGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCg                                                                                                             <  1:1196204/151‑1 (MQ=255)
                                         gCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCga                                                                                                            <  2:713237/151‑1 (MQ=255)
                                         gCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCg                                                                                                             >  2:1176227/1‑150 (MQ=255)
                                           tAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCga                                                                                                            >  2:1139075/1‑149 (MQ=255)
                                            aaaaGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaa                                                                                                           <  1:484010/149‑1 (MQ=255)
                                            aaaaGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaa                                                                                                           <  1:1128195/149‑1 (MQ=255)
                                            aaaaGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCga                                                                                                            >  2:430473/1‑148 (MQ=255)
                                             aaaGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga                                                                                                         >  2:21243/1‑150 (MQ=255)
                                              aaGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGaagaag                                                                                                      <  2:548655/151‑1 (MQ=255)
                                               aGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAgg                                                                                                     <  2:303248/151‑1 (MQ=255)
                                                ggTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAggg                                                                                                    <  1:1089612/151‑1 (MQ=255)
                                                 gTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAgggg                                                                                                   <  2:635100/151‑1 (MQ=255)
                                                    cAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCtt                                                                                               <  1:437482/151‑1 (MQ=255)
                                                     aaaTGCTCAGGACGTTGCAGCGTTTTGCGTGCCCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCAGCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCt                                                                                                <  1:391299/150‑1 (MQ=255)
                                                       aTGCTCAGGACGTTGCAGCGTTTTGCGTGCCCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCATCCCGAAGAAGGGGCCTTg                                                                                              >  2:165355/1‑150 (MQ=255)
                                                       aTGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACACTCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCt                                                                                                >  2:74155/1‑148 (MQ=255)
                                                       aTGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGt                                                                                             <  1:836443/151‑1 (MQ=255)
                                                         gCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGtat                                                                                           >  2:571746/1‑151 (MQ=255)
                                                         gCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGta                                                                                            <  1:45850/150‑1 (MQ=255)
                                                             aGGACGGTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAgg                                                                                       >  1:846552/1‑151 (MQ=255)
                                                               gACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGa                                                                                      >  1:237584/1‑150 (MQ=255)
                                                                     gCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGtt                                                                               <  2:1114053/151‑1 (MQ=255)
                                                                     gCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGt                                                                                <  2:108922/150‑1 (MQ=255)
                                                                      cAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGtt                                                                               >  2:1024345/1‑150 (MQ=255)
                                                                      cAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGtt                                                                               >  1:940281/1‑150 (MQ=255)
                                                                       aGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTat                                                                             >  2:1290587/1‑151 (MQ=255)
                                                                                          cGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcatca                                                            <  1:345254/149‑1 (MQ=255)
                                                                                          cGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTcatcat                                                           <  2:289814/150‑1 (MQ=255)
                                                                                          cGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa                                                          <  2:349847/151‑1 (MQ=255)
                                                                                                  gCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgt                                                  >  1:675180/1‑151 (MQ=255)
                                                                                                    aaaaTTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgt                                                >  1:1145486/1‑151 (MQ=255)
                                                                                                     aaaTTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgtt                                               >  1:58764/1‑151 (MQ=255)
                                                                                                        ttGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTgttgtt                                            >  1:1106190/1‑151 (MQ=255)
                                                                                                         tgcatcagagAAAGCATTGGGTAGGGTCTGGCGCTCATCAACAGTCGGGGAGTCAGCAAGGGGCTGAAACGGGAAAGCCACTAACGAAGAAGGGGGCGTTACAGAAGGAAAGGGGTATGATGAGGCTCGTAATACTGCTAATTGTgttgtt                                            >  1:961893/10‑151 (MQ=255)
                                                                                                         tGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCTCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTgttgtt                                            <  2:167253/150‑1 (MQ=255)
                                                                                                         tGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTgttgtt                                            >  1:327939/1‑150 (MQ=255)
                                                                                                          gCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTa                                           <  2:754697/150‑1 (MQ=255)
                                                                                                             tctGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGt                                       >  2:1211585/1‑151 (MQ=255)
                                                                                                              ctGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTg                                        <  2:1058842/149‑1 (MQ=255)
                                                                                                                gggAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt                                      <  2:512991/149‑1 (MQ=255)
                                                                                                                gggAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTaa                                    <  1:409446/151‑1 (MQ=255)
                                                                                                                  gAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGt                                  <  1:300194/151‑1 (MQ=255)
                                                                                                                  gAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGt                                  <  2:820469/151‑1 (MQ=255)
                                                                                                                        aTTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCg                            <  1:68180/151‑1 (MQ=255)
                                                                                                                           gcgcGGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCTTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGAc                          <  1:392071/150‑1 (MQ=255)
                                                                                                                             gcgGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACtt                        >  1:119000/1‑150 (MQ=255)
                                                                                                                              cgGGGTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGGCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTAc                      >  2:44553/1‑151 (MQ=255)
                                                                                                                                gggTCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACt                     >  2:1287323/1‑150 (MQ=255)
                                                                                                                                 ggTCCGGCGCTCATCAACAATCGGGGTGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTaa                   <  1:82189/151‑1 (MQ=255)
                                                                                                                                   tCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACt                                                        <  2:1292718/112‑1 (MQ=255)
                                                                                                                                   tCCGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACt                                                        >  1:1292679/1‑112 (MQ=255)
                                                                                                                                    ccGGCGCTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTaacaa                <  2:758295/151‑1 (MQ=255)
                                                                                                                                        cgcTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACtca            <  2:408478/151‑1 (MQ=255)
                                                                                                                                         gcTCATCAACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACtcat           <  2:542787/151‑1 (MQ=255)
                                                                                                                                             atcaACAATCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACTCATCaga       <  1:278517/151‑1 (MQ=255)
                                                                                                                                               caacaaTCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACTCATCagag      <  2:56865/150‑1 (MQ=255)
                                                                                                                                                  caaTCGGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACTCATCAGAggggg  >  2:970480/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
CCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAGGCAAAATTGCCTCTGGGAAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACTCATCAGAGGGGG  >  NC_000913/3194922‑3195217

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