Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,430,144 G→C 7.7% 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 7.7% 397.8 / 20.2 117A32G (GCC→GGC) purFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (58/50);  new base C (3/6);  total (61/56)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 119 total aligned reads displayed.

CCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTGGGCA  >  NC_000913/2430003‑2430288
                                                                                                                                             |                                                                                                                                                
ccAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATc                                                                                                                                           >  1:313900/1‑149 (MQ=255)
ccAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATc                                                                                                                                           <  1:829115/149‑1 (MQ=255)
ccAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATc                                                                                                                                           <  1:754832/149‑1 (MQ=255)
 cAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATc                                                                                                                                           >  4:1030818/1‑148 (MQ=255)
 cAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATc                                                                                                                                           >  3:438467/1‑148 (MQ=255)
 cAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCATca                                                                                                                                          <  2:114211/149‑2 (MQ=255)
  aGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCt                                                                                                                                         <  2:801763/149‑1 (MQ=255)
   gCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCATcatc                                                                                                                                        <  2:207261/149‑4 (MQ=255)
   gCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCATcatc                                                                                                                                        <  1:430116/149‑4 (MQ=255)
    ccGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCATcatc                                                                                                                                        <  3:956733/148‑4 (MQ=255)
    ccGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCATcatc                                                                                                                                        <  4:163443/148‑4 (MQ=255)
            tAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATg                                                                                                                               <  1:633382/149‑1 (MQ=255)
                 cacaTGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTg                                                                                                                                        <  3:858480/135‑1 (MQ=255)
                 cacaTGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTg                                                                                                                                        >  4:858480/1‑135 (MQ=255)
                 cacaTGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGaa                                                                                                                          <  1:774207/149‑1 (MQ=255)
                 cacaTGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGa                                                                                                                           <  3:569648/148‑1 (MQ=255)
                   caTGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAg                                                                                                                         >  4:1014232/1‑148 (MQ=255)
                     tGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCAc                                                                                                                      <  2:682596/149‑1 (MQ=255)
                        ccAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCg                                                                                                                    >  3:851540/1‑148 (MQ=255)
                         cAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAGTACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGAAGTTATTGGCTTCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGCAGCACCGt                                                                                                                   >  4:651467/1‑148 (MQ=255)
                          aaTGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGtt                                                                                                                  >  3:272054/1‑148 (MQ=255)
                          aaTGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTa                                                                                                                 >  2:964204/1‑149 (MQ=255)
                            tGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAg                                                                                                               <  2:828291/149‑1 (MQ=255)
                             gCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAgg                                                                                                              <  1:520220/149‑1 (MQ=255)
                              cccATATTGCCCTGCAACCGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGGCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGACGCTGAAGCACCGTTAAGGc                                                                                                             <  2:235592/149‑1 (MQ=255)
                               ccATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGGTGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGc                                                                                                             >  3:1051599/1‑148 (MQ=255)
                               ccATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGca                                                                                                            >  1:427878/1‑149 (MQ=255)
                                 atatTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGcatc                                                                                                          <  2:274579/149‑1 (MQ=255)
                                   atTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGcatcat                                                                                                        >  2:909541/1‑149 (MQ=255)
                                   atTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGcatcat                                                                                                        <  2:313900/149‑1 (MQ=255)
                                      gCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATaaa                                                                                                     <  1:391471/149‑1 (MQ=255)
                                      gCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATaaa                                                                                                     >  1:224987/1‑149 (MQ=255)
                                        ccTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCa                                                                                                                       >  4:169433/1‑129 (MQ=255)
                                        ccTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCa                                                                                                                       <  3:169433/129‑1 (MQ=255)
                                        ccTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCATcatca                                                                                                                                       >  1:1005328/1‑109 (MQ=255)
                                        ccTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCATcatca                                                                                                                                       <  2:1005328/113‑5 (MQ=255)
                                          tGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCg                                                                                                  >  3:755924/1‑148 (MQ=255)
                                          tGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCg                                                                                                  >  4:844519/1‑148 (MQ=255)
                                          tGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                                 >  1:997387/1‑149 (MQ=255)
                                           gCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                                 <  3:226421/148‑1 (MQ=255)
                                           gCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGAc                                                                                                <  2:899495/149‑1 (MQ=255)
                                             aaaCGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACt                                                                                               <  4:895024/148‑1 (MQ=255)
                                             aaaCGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTg                                                                                              >  2:959700/1‑149 (MQ=255)
                                                      atatGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCgg                                                                                     <  2:883924/149‑1 (MQ=255)
                                                       tatGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCgg                                                                                     >  3:62377/1‑148 (MQ=255)
                                                           gCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATaa                                                                                >  1:825216/1‑149 (MQ=255)
                                                                cTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGATGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCg                                                                           >  1:597367/1‑149 (MQ=255)
                                                                 ttCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCgg                                                                          <  2:467108/149‑1 (MQ=255)
                                                                 ttCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCg                                                                           <  4:159072/148‑1 (MQ=255)
                                                                  tCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCgg                                                                          >  4:978157/1‑148 (MQ=255)
                                                                  tCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGc                                                                         >  1:129668/1‑149 (MQ=255)
                                                                     aaTACATCGCTCACCAGCGCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                                 <  1:822408/122‑1 (MQ=255)
                                                                     aaTACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                                 >  2:822408/1‑122 (MQ=255)
                                                                     aaTACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGa                                                                       <  3:880100/148‑1 (MQ=255)
                                                                         cATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATAcc                                                                   >  4:717393/1‑148 (MQ=255)
                                                                          aTCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGTTTAACCGGCATAACACCGGCGATACCGa                                                                 >  2:1061901/1‑149 (MQ=255)
                                                                          aTCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGa                                                                 <  2:369595/149‑1 (MQ=255)
                                                                           tCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTa                                                                                                                 <  4:562681/100‑1 (MQ=255)
                                                                           tCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTa                                                                                                                 >  3:562681/1‑100 (MQ=255)
                                                                           tCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGa                                                                 >  3:19388/1‑148 (MQ=255)
                                                                           tCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGa                                                                 >  4:1029503/1‑148 (MQ=255)
                                                                           tCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGa                                                                 <  4:349501/148‑1 (MQ=255)
                                                                           tCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGa                                                                 >  4:271511/1‑148 (MQ=255)
                                                                             gCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACa                                                               <  4:438467/148‑1 (MQ=255)
                                                                              cTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGcatc                                                                                                          <  3:348750/104‑1 (MQ=255)
                                                                               tCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATa                                                            >  2:674465/1‑149 (MQ=255)
                                                                                 aCCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATAcc                                                          >  2:867412/1‑149 (MQ=255)
                                                                                 aCCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATAc                                                           <  4:25525/148‑1 (MQ=255)
                                                                                    aGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGc                                                        <  3:896753/148‑1 (MQ=255)
                                                                                        cGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTTACCGCGATGCTGAAGCACCATTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACAt                                                    <  3:651467/148‑1 (MQ=255)
                                                                                         gTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATa                                                   <  3:38933/148‑1 (MQ=255)
                                                                                          ttCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGcatc                                                                                                          >  1:1055028/1‑92 (MQ=255)
                                                                                          ttCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGcatc                                                                                                          <  2:1055028/92‑1 (MQ=255)
                                                                                          ttCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACg                                                 >  1:622293/1‑149 (MQ=255)
                                                                                             gCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTc                                               <  4:851540/148‑1 (MQ=255)
                                                                                                   cGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCAGGCGGCAGCCTGACCGCGATGCTGAAGTACCGTTAAGGCATCATAAATCGCCTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACAGACGTCTTTTTc                                         >  4:654128/1‑148 (MQ=255)
                                                                                                    gCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCt                                       >  1:840909/1‑149 (MQ=255)
                                                                                                      aaaCGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTc                                      >  3:960921/1‑148 (MQ=255)
                                                                                                        aCGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGt                                    <  3:691618/148‑1 (MQ=255)
                                                                                                        aCGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGt                                    >  4:718622/1‑148 (MQ=255)
                                                                                                          ggAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTaa                                 <  1:730253/149‑1 (MQ=255)
                                                                                                                 gTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCcaca                           <  3:622951/148‑1 (MQ=255)
                                                                                                                     ttGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCa                      <  2:867844/149‑1 (MQ=255)
                                                                                                                       gggaaccaTGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCag                     <  4:617831/141‑1 (MQ=255)
                                                                                                                       ggCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCaga                    <  1:14876/149‑1 (MQ=255)
                                                                                                                       ggCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCag                     <  3:1029503/148‑1 (MQ=255)
                                                                                                                         cATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCagag                   >  4:565894/1‑148 (MQ=255)
                                                                                                                          aTCTATGTTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCa                 >  2:586016/1‑149 (MQ=255)
                                                                                                                          aTCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGc                  <  3:33043/148‑1 (MQ=255)
                                                                                                                           tCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCa                 >  3:736223/1‑148 (MQ=255)
                                                                                                                             tATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACtt              <  1:724522/149‑1 (MQ=255)
                                                                                                                              aTGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTa             >  2:803473/1‑149 (MQ=255)
                                                                                                                                  tgatgaTGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGct         >  2:516572/1‑149 (MQ=255)
                                                                                                                                  tcaggatgCCGCCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGct         >  2:392422/5‑149 (MQ=255)
                                                                                                                                      gatgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTg      <  4:62377/148‑1 (MQ=255)
                                                                                                                                      gatgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTg      <  3:235409/148‑1 (MQ=255)
                                                                                                                                        tgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTa             >  2:17433/1‑139 (MQ=255)
                                                                                                                                        tgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTa             <  1:17433/139‑1 (MQ=255)
                                                                                                                                             gcggcATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCag                     >  4:1055381/1‑126 (MQ=255)
                                                                                                                                             gcggcATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCag                     <  3:1055381/126‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
CCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTGGGCA  >  NC_000913/2430003‑2430288

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: