Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,895,391 A→T 8.8% K196I (AAA→ATA) ‡ pabB → aminodeoxychorismate synthase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,895,3910AT 8.8% 343.1 / 20.0 102K196I (AAA→ATA) ‡pabBaminodeoxychorismate synthase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (52/41);  new base T (7/2);  total (59/43)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

Only 100 of 106 total aligned reads displayed.

CACCAGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTT  >  NC_000913/1895243‑1895538
                                                                                                                                                    |                                                                                                                                                   
caccaGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGaaaa                                                                                                                                                     >  3:540579/1‑149 (MQ=255)
  ccaGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGaaaaa                                                                                                                                                    <  1:457526/148‑1 (MQ=255)
  ccaGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAAt                                                                                                                                                   >  3:142014/1‑149 (MQ=255)
       gTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGc                                                                                                                                              <  3:122126/149‑1 (MQ=255)
         cATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGcc                                                                                                                                             >  1:343163/1‑148 (MQ=255)
            aCAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAgg                                                                                                                                          >  1:726883/1‑148 (MQ=255)
               gTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTAc                                                                                                                                       >  1:273980/1‑148 (MQ=255)
               gTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACa                                                                                                                                      >  4:645954/1‑149 (MQ=255)
                tttCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAg                                                                                                                                     >  3:385338/1‑149 (MQ=255)
                 ttCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAgg                                                                                                                                    <  3:461052/149‑1 (MQ=255)
                              cATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGcaca                                                                                                                        >  2:116882/1‑148 (MQ=255)
                               aTAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGc                                                                                                                      >  3:91422/1‑149 (MQ=255)
                                tAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGc                                                                                                                      <  2:523345/148‑1 (MQ=255)
                                tAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAACCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCg                                                                                                                     >  4:353797/1‑149 (MQ=255)
                                   tgatgTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGt                                                                                                                   >  1:394706/1‑148 (MQ=255)
                                   tgatgTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTg                                                                                                                  >  3:129218/1‑149 (MQ=255)
                                    gatgTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGa                                                                                                                 <  3:645954/149‑1 (MQ=255)
                                     atgTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAATTTTTCGccg                                                                                                                                            >  1:173042/1‑120 (MQ=255)
                                     atgTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAATTTTTCGccg                                                                                                                                            <  2:173042/121‑2 (MQ=255)
                                      tgTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGAt                                                                                                                <  2:406777/148‑1 (MQ=255)
                                         cAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAg                                                                                                                                     >  3:798568/1‑124 (MQ=255)
                                         cAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAg                                                                                                                                     <  4:798568/124‑1 (MQ=255)
                                          aaTGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCt                                                                                                            >  1:617657/1‑148 (MQ=255)
                                          aaTGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCt                                                                                                            >  2:169451/1‑148 (MQ=255)
                                                cGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAgg                                                                                                      >  1:607065/1‑148 (MQ=255)
                                                 gTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAat                                                                                                                                 >  4:25827/1‑120 (MQ=255)
                                                 gTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAat                                                                                                                                 <  3:25827/120‑1 (MQ=255)
                                                   cGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGaa                                                                                                  <  4:463938/149‑1 (MQ=255)
                                                      gCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAAt                                                                                                                                                   >  3:876750/1‑97 (MQ=255)
                                                         tggctggAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAGTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCg                                                                                             >  1:676045/1‑148 (MQ=255)
                                                         tggctggAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCg                                                                                             >  2:31344/1‑148 (MQ=255)
                                                          ggctggAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGc                                                                                                             <  4:582820/131‑1 (MQ=255)
                                                          ggctggAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGc                                                                                                             >  3:582820/1‑131 (MQ=255)
                                                          ggctggAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGcc                                                                                           <  4:540579/149‑1 (MQ=255)
                                                              ggAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCaa                                                                                        >  2:180935/1‑148 (MQ=255)
                                                                    ccAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTa                                                                                                           <  3:46611/123‑1 (MQ=255)
                                                                    ccAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTa                                                                                                           >  4:46611/1‑123 (MQ=255)
                                                                    ccAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTc                                                                                 >  3:443156/1‑149 (MQ=255)
                                                                    ccAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTc                                                                                 <  4:145692/149‑1 (MQ=255)
                                                                     cAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGAt                                                                                                                <  3:222687/117‑1 (MQ=255)
                                                                     cAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGAt                                                                                                                >  4:222687/1‑117 (MQ=255)
                                                                     cAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTAt                                                                                                          <  4:456516/123‑1 (MQ=255)
                                                                     cAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTAt                                                                                                          >  3:456516/1‑123 (MQ=255)
                                                                      aGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCAt                                                                               >  3:810726/1‑149 (MQ=255)
                                                                      aGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCAt                                                                               >  3:607025/1‑149 (MQ=255)
                                                                       gCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATg                                                                              >  3:539606/1‑149 (MQ=255)
                                                                        cAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATg                                                                              <  2:343163/148‑1 (MQ=255)
                                                                        cAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGc                                                                             <  4:31225/149‑1 (MQ=255)
                                                                          aTTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCg                                                                            <  1:331358/148‑1 (MQ=255)
                                                                          aTTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGa                                                                           <  4:234439/149‑1 (MQ=255)
                                                                                   gcAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCt                                                                   >  2:479989/1‑148 (MQ=255)
                                                                                    cAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCAGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTg                                                                  >  2:774329/1‑148 (MQ=255)
                                                                                        aaGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGa                                                              >  2:743668/1‑148 (MQ=255)
                                                                                        aaGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGa                                                              <  2:607065/148‑1 (MQ=255)
                                                                                         aGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGcc                                                                                                                                             <  2:804741/68‑1 (MQ=255)
                                                                                         aGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGcc                                                                                                                                             >  1:804741/1‑68 (MQ=255)
                                                                                              tCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATg                                                        <  2:570124/148‑1 (MQ=255)
                                                                                               cACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggc                                                      >  3:713056/1‑149 (MQ=255)
                                                                                                aCGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggca                                                     <  4:658124/149‑1 (MQ=255)
                                                                                                aCGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggc                                                      >  2:796376/1‑148 (MQ=255)
                                                                                                aCGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggc                                                      <  1:355101/148‑1 (MQ=255)
                                                                                                      aCTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCAt                                                <  2:236297/148‑1 (MQ=255)
                                                                                                         tCCGACTGGCAATCCAATATGACCCTCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTcc                                             >  2:354375/1‑148 (MQ=255)
                                                                                                          ccGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCt                                            <  1:116882/148‑1 (MQ=255)
                                                                                                            gACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTc                                          <  1:750253/148‑1 (MQ=255)
                                                                                                              cTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAg                                        <  2:544112/148‑1 (MQ=255)
                                                                                                              cTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGc                                       <  4:502719/149‑1 (MQ=255)
                                                                                                                   aaTCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAAt                                  >  4:444532/1‑149 (MQ=255)
                                                                                                                      ccAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCa                                >  1:864065/1‑148 (MQ=255)
                                                                                                                      ccAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAg                               <  4:385338/149‑1 (MQ=255)
                                                                                                                       cAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAgg                              <  3:353797/149‑1 (MQ=255)
                                                                                                                        aaTATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTg                                                                  <  4:372537/112‑1 (MQ=255)
                                                                                                                        aaTATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTg                                                                  >  3:372537/1‑112 (MQ=255)
                                                                                                                        aaTATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGAt                                                             <  3:471928/117‑1 (MQ=255)
                                                                                                                        aaTATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGAt                                                             >  4:471928/1‑117 (MQ=255)
                                                                                                                         atatGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCt                                      >  1:805567/1‑139 (MQ=255)
                                                                                                                         atatGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCt                                      <  2:805567/139‑1 (MQ=255)
                                                                                                                         atatGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTCATCAGGc                             >  1:37208/1‑148 (MQ=255)
                                                                                                                             gACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAAc                         <  2:726883/148‑1 (MQ=255)
                                                                                                                              aCCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAAcc                        >  2:891369/1‑148 (MQ=255)
                                                                                                                               cccGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACcgc                      <  3:358170/149‑1 (MQ=255)
                                                                                                                               cccGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACcg                       >  1:173968/1‑148 (MQ=255)
                                                                                                                                ccGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACcgc                      >  2:634966/1‑148 (MQ=255)
                                                                                                                                 cgcgAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTa                                                                                                           >  1:775419/1‑62 (MQ=255)
                                                                                                                                 cgcgAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTa                                                                                                           <  2:775419/62‑1 (MQ=255)
                                                                                                                                 cgcgAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATg                                                        <  1:815346/113‑1 (MQ=255)
                                                                                                                                 cgcgAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAg                                        <  2:730968/129‑1 (MQ=255)
                                                                                                                                         gTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTa            <  4:810726/149‑1 (MQ=255)
                                                                                                                                           ctggCGAAATTTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAgc          >  3:478079/3‑149 (MQ=255)
                                                                                                                                            tggCGAAATTTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAgcg         >  3:104822/2‑149 (MQ=255)
                                                                                                                                             ggCGAAATTTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAgcg         >  2:541449/1‑148 (MQ=255)
                                                                                                                                             ggCGAAATTTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAgcg         >  1:440135/1‑148 (MQ=255)
                                                                                                                                             ggCGAAATTTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAgcg         >  1:147225/1‑148 (MQ=255)
                                                                                                                                              gCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCt       >  4:148903/1‑149 (MQ=255)
                                                                                                                                               cGAAATTTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCAttt             <  2:639106/142‑1 (MQ=255)
                                                                                                                                               cGAAATTTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCAttt             >  1:102056/1‑142 (MQ=255)
                                                                                                                                               cGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCt       >  2:142137/1‑148 (MQ=255)
                                                                                                                                                gAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCtttt    >  4:66431/1‑149 (MQ=255)
                                                                                                                                                 aaaaaTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCtttt    >  4:584868/1‑149 (MQ=255)
                                                                                                                                                    aaTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCtttttt  <  2:676045/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CACCAGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTT  >  NC_000913/1895243‑1895538

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