Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,895,392 A→T 8.9% K196N (AAA→AAT) ‡ pabB → aminodeoxychorismate synthase, subunit I

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

Only 100 of 107 total aligned reads displayed.

CCAGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTA  >  NC_000913/1895245‑1895539
                                                                                                                                                   |                                                                                                                                                   
ccaGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGaaaaa                                                                                                                                                     <  1:457526/148‑1 (MQ=255)
ccaGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAAt                                                                                                                                                    >  3:142014/1‑149 (MQ=255)
     gTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGc                                                                                                                                               >  3:502719/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)
               ttCTTTTCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAg                                                                                                                                      >  2:490968/1‑148 (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)
                                   atgTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAATTTTTCGccg                                                                                                                                             <  2:173042/121‑2 (MQ=255)
                                   atgTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAATTTTTCGccg                                                                                                                                             >  1:173042/1‑120 (MQ=255)
                                    tgTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGAt                                                                                                                 <  2:406777/148‑1 (MQ=255)
                                       cAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAg                                                                                                                                      >  3:798568/1‑124 (MQ=255)
                                       cAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAg                                                                                                                                      <  4:798568/124‑1 (MQ=255)
                                        aaTGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCt                                                                                                             >  2:169451/1‑148 (MQ=255)
                                              cGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAgg                                                                                                       >  1:607065/1‑148 (MQ=255)
                                               gTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAat                                                                                                                                  >  4:25827/1‑120 (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)
                                                                  ccAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTa                                                                                                            >  4:46611/1‑123 (MQ=255)
                                                                  ccAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTa                                                                                                            <  3:46611/123‑1 (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                                                                                                           >  3:456516/1‑123 (MQ=255)
                                                                   cAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTAt                                                                                                           <  4:456516/123‑1 (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:607065/148‑1 (MQ=255)
                                                                                      aaGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGa                                                               >  2:743668/1‑148 (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                                                       <  1:355101/148‑1 (MQ=255)
                                                                                              aCGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggc                                                       >  2:796376/1‑148 (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)
                                                                                                                     cAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAgg                               <  3:353797/149‑1 (MQ=255)
                                                                                                                      aaTATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTg                                                                   >  3:372537/1‑112 (MQ=255)
                                                                                                                      aaTATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTg                                                                   <  4:372537/112‑1 (MQ=255)
                                                                                                                      aaTATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGAt                                                              <  3:471928/117‑1 (MQ=255)
                                                                                                                      aaTATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGAt                                                              >  4:471928/1‑117 (MQ=255)
                                                                                                                       atatGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCt                                       <  2:805567/139‑1 (MQ=255)
                                                                                                                       atatGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCt                                       >  1:805567/1‑139 (MQ=255)
                                                                                                                       atatGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTCATCAGGc                              >  1:37208/1‑148 (MQ=255)
                                                                                                                          tGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAAc                          <  4:365040/149‑1 (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                                                                                                            <  2:775419/62‑1 (MQ=255)
                                                                                                                               cgcgAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTa                                                                                                            >  1:775419/1‑62 (MQ=255)
                                                                                                                               cgcgAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATg                                                         >  2:815346/1‑113 (MQ=255)
                                                                                                                               cgcgAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATg                                                         <  1:815346/113‑1 (MQ=255)
                                                                                                                               cgcgAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAg                                         >  1:730968/1‑129 (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)
                                                                                                                                                   aTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTa  <  1:743668/148‑1 (MQ=255)
                                                                                                                                                   aTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTa  >  1:588224/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCAGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCCGTCGGGCCTGGCTGGAAAGCCAGCAATTCTCGCCGCAGGAAGATTTCACGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTA  >  NC_000913/1895245‑1895539

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