Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 769,691 C→T P572S (CCG→TCG)  mngB → alpha‑mannosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913769,6910CT98.5% 246.4 / ‑5.2 67P572S (CCG→TCG) mngBalpha‑mannosidase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (33/33);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

TGAAGAAAGCTCTGATGATGGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGG  >  NC_000913/769543‑769838
                                                                                                                                                    |                                                                                                                                                   
tGAAGAAAGCTCTGATGATGGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCg                                                                                                                                                   >  2:396387/1‑151 (MQ=255)
 gaagaaAGCTCTGATGATGGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGc                                                                                                                                                  >  2:924632/1‑151 (MQ=255)
     aaaGCTCTGATGATGGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCaa                                                                                                                                              <  1:788586/151‑1 (MQ=255)
        gCTCTGATGATGGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCGGTGTCGCTCAAtt                                                                                                                                            >  1:696093/1‑150 (MQ=255)
            tgatgaTGGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTc                                                                                                                                        >  1:349501/1‑150 (MQ=255)
              atgatgGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGa                                                                                                                                     <  2:850864/151‑1 (MQ=255)
                 atgGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACg                                                                                                                                  <  1:841725/151‑1 (MQ=255)
                 atgGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACg                                                                                                                                  >  1:144772/1‑151 (MQ=255)
                   ggTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCa                                                                                                                                <  1:777718/151‑1 (MQ=255)
                      gatgaATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAgcg                                                                                                                             <  1:161750/151‑1 (MQ=255)
                      gatgaATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAgcg                                                                                                                             <  1:4744/151‑1 (MQ=255)
                             atGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGc                                                                                                                       >  1:288142/1‑150 (MQ=255)
                                   aTTCACCCGCAAATGAAGAGTGGGCAATTACCGCAGCGCACGCGAAACCGCAATGCGATATTATTCGTGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACACGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCAGGCAATCCAc                                                                                                                >  2:726508/1‑151 (MQ=255)
                                   aTTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATcc                                                                                                                  <  1:629955/149‑1 (MQ=255)
                                     tCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACt                                                                                                               <  2:308936/150‑1 (MQ=255)
                                       aCCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGc                                                                                                            >  1:1017685/1‑151 (MQ=255)
                                             aaaaGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCag                                                                                                          >  2:419687/1‑147 (MQ=255)
                                             aaaaGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTa                                                                                                      <  1:48800/151‑1 (MQ=255)
                                              aaaGCATAGTGGGTAATTACCGACGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGACCAGGGCTGTTATCCGCCATGACATCTCAGTGTCGCTCAATTTGTCAGAACGCACCGCCCGGCAATCCACTgcccgcg                                                                                                        >  2:764711/1‑142 (MQ=255)
                                               aagaagAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAgg                                                                                                    >  1:647263/1‑151 (MQ=255)
                                                  aagaGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAgggg                                                                                                  >  1:829900/1‑150 (MQ=255)
                                                         ggTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTgttg                                                                                          >  1:446875/1‑151 (MQ=255)
                                                            aaTTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTgttggt                                                                                        <  1:1063036/150‑1 (MQ=255)
                                                             aTTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTgttggttg                                                                                      <  1:232913/151‑1 (MQ=255)
                                                               tACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTc                                                                                    >  2:177684/1‑151 (MQ=255)
                                                                aCCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTc                                                                                    >  1:632342/1‑150 (MQ=255)
                                                                  cGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCAct                                                                                 <  1:481410/151‑1 (MQ=255)
                                                                  cGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCAct                                                                                 >  1:1017741/1‑151 (MQ=255)
                                                                       cGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCtt                                                                              >  2:358830/1‑149 (MQ=255)
                                                                       cGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAg                                                                            <  2:981066/151‑1 (MQ=255)
                                                                        gAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAg                                                                            >  1:765709/1‑150 (MQ=255)
                                                                        gAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGt                                                                           >  1:556168/1‑151 (MQ=255)
                                                                        gAACGCGAAACCGCAATGCGAAATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGt                                                                           <  1:755853/151‑1 (MQ=255)
                                                                         aaCGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGt                                                                           <  2:365244/150‑1 (MQ=255)
                                                                          aagcgAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTATCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGt                                                                           >  1:333794/3‑149 (MQ=255)
                                                                           cgcgAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCAt                                                                        <  1:238197/151‑1 (MQ=255)
                                                                           cgcgAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCAt                                                                        <  1:376655/151‑1 (MQ=255)
                                                                            gcgAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATa                                                                       <  2:589050/151‑1 (MQ=255)
                                                                               aaaCCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAAc                                                                     >  2:576667/1‑150 (MQ=255)
                                                                                   cGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAAcagcag                                                                <  2:427660/151‑1 (MQ=255)
                                                                                   cGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAAcagc                                                                  >  2:756110/1‑149 (MQ=255)
                                                                                    gCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAAcagcag                                                                <  2:954178/150‑1 (MQ=255)
                                                                                      aaTGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTg                                                                                                              <  1:966938/102‑1 (MQ=255)
                                                                                      aaTGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTg                                                                                                              >  2:966889/1‑102 (MQ=255)
                                                                                      aaTGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCg                                                             >  1:51169/1‑151 (MQ=255)
                                                                                      aaTGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCg                                                             <  1:945963/151‑1 (MQ=255)
                                                                                      aaTGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCg                                                             <  1:27191/151‑1 (MQ=255)
                                                                                         gCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTa                                                           <  2:275798/150‑1 (MQ=255)
                                                                                         gCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTAt                                                          >  1:100233/1‑151 (MQ=255)
                                                                                         gCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTAt                                                          <  1:609011/151‑1 (MQ=255)
                                                                                            atatTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGGTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTg                                                        >  2:465212/1‑150 (MQ=255)
                                                                                            atatTATTCATGAAGCCTGGCAGAGCAGGGCTGTTAGCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTa                                                           <  2:22465/147‑1 (MQ=255)
                                                                                                   tCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGa                                                 <  2:502848/150‑1 (MQ=255)
                                                                                                   tCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGa                                                 >  2:659541/1‑150 (MQ=255)
                                                                                                         aGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAAt                                          >  1:154811/1‑151 (MQ=255)
                                                                                                          gCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATc                                         >  2:387314/1‑151 (MQ=255)
                                                                                                                cAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGAt                                    >  1:928059/1‑150 (MQ=255)
                                                                                                                     cAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAAcc                                >  1:1048882/1‑149 (MQ=255)
                                                                                                                         gCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTg                          <  1:662107/151‑1 (MQ=255)
                                                                                                                               aTCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGatc                    >  2:582994/1‑151 (MQ=255)
                                                                                                                                 ccGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGatca                   <  1:262406/150‑1 (MQ=255)
                                                                                                                                      aTGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCt             <  1:1017417/151‑1 (MQ=255)
                                                                                                                                       tGACATGGCAGTGTCGCGCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGGTGACGATCATCGcc              <  2:972399/149‑1 (MQ=255)
                                                                                                                                          cATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCgt         <  1:516044/151‑1 (MQ=255)
                                                                                                                                           aTGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCTCTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCgt         <  2:542775/150‑1 (MQ=255)
                                                                                                                                           aTGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCgtg        <  2:136886/151‑1 (MQ=255)
                                                                                                                                                aGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTg   >  2:662280/1‑151 (MQ=255)
                                                                                                                                                 gtgtCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTgg  <  2:749710/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGAAGAAAGCTCTGATGATGGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGG  >  NC_000913/769543‑769838

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