Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,430,144 G→C 6.5% A32G (GCC→GGC)  purF ← amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,430,1440GC 6.5% 274.0 / 8.6 77A32G (GCC→GGC) purFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (37/35);  new base C (3/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

CCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCT  >  NC_000913/2430003‑2430283
                                                                                                                                             |                                                                                                                                           
ccAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCAt                                                                                                                                       <  2:838958/148‑1 (MQ=255)
                 cacaTGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCAt                                                                                                                                       <  2:228650/131‑1 (MQ=255)
                 cacaTGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCAt                                                                                                                                       >  1:228650/1‑131 (MQ=255)
                   caTGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAg                                                                                                                    >  3:773459/1‑148 (MQ=255)
                         cAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGt                                                                                                              <  4:773459/148‑1 (MQ=255)
                                 atatTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTTGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGa                                                                                                                      >  2:371596/1‑132 (MQ=255)
                                 atatTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCggcggc                                                                                                                                         <  3:544672/113‑1 (MQ=255)
                                 atatTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCggcggc                                                                                                                                         >  4:544672/1‑113 (MQ=255)
                                 atatTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGa                                                                                                                      <  1:371596/132‑1 (MQ=255)
                                        ccTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATc                                                                                                                                      <  2:507700/109‑1 (MQ=255)
                                        ccTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATc                                                                                                                                      >  1:507700/1‑109 (MQ=255)
                                               aCGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGa                                                                                                                                  >  3:466460/1‑106 (MQ=255)
                                               aCGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGa                                                                                                                                  <  4:466460/106‑1 (MQ=255)
                                               aCGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCAc                                                                                                                 <  1:783832/123‑1 (MQ=255)
                                               aCGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCAc                                                                                                                 >  2:783832/1‑123 (MQ=255)
                                                     cATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAAcc                                                                                  >  1:1091765/1‑148 (MQ=255)
                                                       tatGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGaa                                                                                                                     <  4:53137/111‑1 (MQ=255)
                                                       tatGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGaa                                                                                                                     >  3:53137/1‑111 (MQ=255)
                                                       tatGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTaa                                                                                    <  3:128545/144‑1 (MQ=255)
                                                       tatGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTaa                                                                                    >  4:128545/1‑144 (MQ=255)
                                                        atGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGc                                                                               <  2:1091765/148‑1 (MQ=255)
                                                          ggCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCAt                                                                             >  2:478572/1‑148 (MQ=255)
                                                            cGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGcat                                                                                                      <  1:63621/121‑1 (MQ=255)
                                                            cGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGcat                                                                                                      >  2:63621/1‑121 (MQ=255)
                                                                  tCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCgg                                                                     <  1:478572/148‑1 (MQ=255)
                                                                     aaTACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTa                                                                                                            >  4:335190/1‑106 (MQ=255)
                                                                     aaTACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTa                                                                                                            <  3:335190/106‑1 (MQ=255)
                                                                      aTACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAAcc                                                                                  <  4:74426/131‑1 (MQ=255)
                                                                      aTACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAAcc                                                                                  >  3:74426/1‑131 (MQ=255)
                                                                           tCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                            <  3:464982/116‑1 (MQ=255)
                                                                           tCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                            >  4:464982/1‑116 (MQ=255)
                                                                           tCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGa                                                            >  1:991586/1‑148 (MQ=255)
                                                                            cGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCAc                                                                                                                 >  3:520705/1‑94 (MQ=255)
                                                                            cGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCAc                                                                                                                 <  4:520705/94‑1 (MQ=255)
                                                                                  ccAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                            <  2:934229/109‑1 (MQ=255)
                                                                                  ccAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                            >  1:934229/1‑109 (MQ=255)
                                                                                   cAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                            <  1:1051307/108‑1 (MQ=255)
                                                                                   cAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                            >  2:1051307/1‑108 (MQ=255)
                                                                                   cAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCg                                                    <  2:991586/148‑1 (MQ=255)
                                                                                     gCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGca                                                  >  4:860139/1‑148 (MQ=255)
                                                                                               tttACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGAt                                                                 >  1:105091/1‑123 (MQ=255)
                                                                                               tttACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGAt                                                                 <  2:105091/123‑1 (MQ=255)
                                                                                                    gCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTcc                                   <  3:860139/148‑1 (MQ=255)
                                                                                                     cAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAac                                                                          <  3:247624/108‑1 (MQ=255)
                                                                                                     cAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAac                                                                          >  4:247624/1‑108 (MQ=255)
                                                                                                     cAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCt                                  >  3:759162/1‑148 (MQ=255)
                                                                                                        aCGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGt                               >  1:268658/1‑148 (MQ=255)
                                                                                                             aGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCgg                                                                     >  2:135412/1‑105 (MQ=255)
                                                                                                             aGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCgg                                                                     <  1:135412/105‑1 (MQ=255)
                                                                                                                 gTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCttt                                       <  4:422679/131‑1 (MQ=255)
                                                                                                                 gTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCttt                                       >  3:422679/1‑131 (MQ=255)
                                                                                                                  ttATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACAt                     >  1:269420/1‑148 (MQ=255)
                                                                                                                     ttGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATa                                                                >  3:1089061/1‑102 (MQ=255)
                                                                                                                     ttGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATa                                                                >  1:226008/1‑102 (MQ=255)
                                                                                                                     ttGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATa                                                                <  2:226008/102‑1 (MQ=255)
                                                                                                                     ttGGCATCAATGGTGATGATGCCGGCGGCATCCTGACCGCGCTGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATa                                                                <  4:1089061/102‑1 (MQ=255)
                                                                                                                           tCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGt                                                                                       >  3:218640/1‑73 (MQ=255)
                                                                                                                           tCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGt                                                                                       <  4:218640/73‑1 (MQ=255)
                                                                                                                              aTGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACtt         >  4:381026/1‑148 (MQ=255)
                                                                                                                               tGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGGCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTa        <  2:269420/148‑1 (MQ=255)
                                                                                                                                 gTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTa                             <  3:416727/125‑1 (MQ=255)
                                                                                                                                 gTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTa                             <  1:291567/125‑1 (MQ=255)
                                                                                                                                 gTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTa                             >  2:291567/1‑125 (MQ=255)
                                                                                                                                 gTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTa                             >  4:416727/1‑125 (MQ=255)
                                                                                                                                  tgatgaTGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGc                                                                                                        >  3:876181/1‑49 (MQ=255)
                                                                                                                                  tgatgaTGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGc                                                                                                        <  4:876181/49‑1 (MQ=255)
                                                                                                                                   caggatgCCGCCGGCATCATGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGct    >  2:1062821/4‑148 (MQ=255)
                                                                                                                                    atgatgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACAt                                               <  2:952205/104‑1 (MQ=255)
                                                                                                                                    atgatgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACAt                                               >  1:952205/1‑104 (MQ=255)
                                                                                                                                    atgatgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCag                >  2:926394/1‑135 (MQ=255)
                                                                                                                                    atgatgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCag                <  1:926394/135‑1 (MQ=255)
                                                                                                                                      gatgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGt                                                                                       <  3:865249/62‑1 (MQ=255)
                                                                                                                                      gatgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGt                                                                                       >  4:865249/1‑62 (MQ=255)
                                                                                                                                      gatgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCttttt                                     <  2:862698/112‑1 (MQ=255)
                                                                                                                                      gatgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCttttt                                     >  1:862698/1‑112 (MQ=255)
                                                                                                                                       atgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCt                                  <  4:437169/114‑1 (MQ=255)
                                                                                                                                       atgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCt                                  >  3:437169/1‑114 (MQ=255)
                                                                                                                                       atgCCGCCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGctct  >  4:329116/1‑146 (MQ=255)
                                                                                                                                       atgCCGCCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGctct  <  3:329116/146‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                           
CCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCT  >  NC_000913/2430003‑2430283

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: