Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,895,489 G→C 6.7% A229P (GCA→CCA)  pabB → aminodeoxychorismate synthase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,895,4890GC 6.7% 196.5 / 5.5 60A229P (GCA→CCA) pabBaminodeoxychorismate synthase, subunit I
Reads supporting (aligned to +/- strand):  ref base G (23/32);  new base C (2/2);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTAAAGGCACGCTAC  >  NC_000913/1895341‑1895636
                                                                                                                                                    |                                                                                                                                                   
gCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggcagg                                                                                                                                                     >  2:170065/1‑149 (MQ=255)
 cTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggcagg                                                                                                                                                     >  8:28970/1‑148 (MQ=255)
  tCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggcaggca                                                                                                                                                   >  2:53202/1‑149 (MQ=255)
     cTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTc                                                                                                                                                <  8:44342/149‑1 (MQ=255)
     cTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTc                                                                                                                                                >  2:309716/1‑149 (MQ=255)
     cTTCCGACTGGCAATACAATATGACCCGGGAGCAGTACGGTGAAAATTTTCGCCAGGTACAGGAATATCTGCAGAGCGGTTATTTCTATCAGGATTATCTCGCCCATCGTTTTCCTTAGACCTAATCTGCCGATGAATGGCAGGCATTc                                                                                                                                                <  6:122403/149‑1 (MQ=255)
         cGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCtt                                                                                                                                             <  2:321200/148‑1 (MQ=255)
          gACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCa                                                                                                                                           <  5:98608/149‑1 (MQ=255)
             tGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCt                                                                                                                                        >  1:277951/1‑149 (MQ=255)
               gCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCtt                                                                                                                                       >  5:38339/1‑148 (MQ=255)
               gCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCtt                                                                                                                                       >  2:231723/1‑148 (MQ=255)
               gCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCtt                                                                                                                                       >  6:82898/1‑148 (MQ=255)
               gCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCtt                                                                                                                                       >  4:158164/1‑148 (MQ=255)
                 aaTCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAAt                                                                                                                                    >  3:283227/1‑149 (MQ=255)
                  aTCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAAt                                                                                                                                    >  4:80322/1‑148 (MQ=255)
                   tCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGGGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCa                                                                                                                                  <  5:241178/149‑1 (MQ=255)
                     cAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAg                                                                                                                                 <  8:12563/148‑1 (MQ=255)
                       atatGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGc                                                                                                                               >  6:69333/1‑148 (MQ=255)
                          tGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCaa                                                                                                                            <  7:173850/148‑1 (MQ=255)
                           gACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAAcc                                                                                                                          >  5:55885/1‑149 (MQ=255)
                               cgcgAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACcgcgc                                                                                                                      <  1:53202/149‑1 (MQ=255)
                                          cGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAgc                                                                                                            <  1:309716/148‑1 (MQ=255)
                                           ggCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAgcgc                                                                                                          <  6:254464/149‑1 (MQ=255)
                                                   aTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTAc                                                                                                  >  7:281805/1‑149 (MQ=255)
                                                       cGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTc                                                                                               <  1:55093/148‑1 (MQ=255)
                                                             gTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGaa                                                                                          >  1:11923/1‑147 (MQ=255)
                                                             gTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACa                                                                                        <  7:28970/149‑1 (MQ=255)
                                                                     atatCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCaa                                                                                <  3:80322/149‑1 (MQ=255)
                                                                        tCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCtt                                                                                             <  3:93816/133‑1 (MQ=255)
                                                                        tCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCtt                                                                                             >  4:93816/1‑133 (MQ=255)
                                                                        tCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAAttt                                                                             >  2:136848/1‑149 (MQ=255)
                                                                            cacaGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCTACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTaa                                                                          <  4:101430/148‑1 (MQ=255)
                                                                            cacaGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTaa                                                                          <  3:158164/148‑1 (MQ=255)
                                                                            cacaGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTa                                                                           <  5:82898/147‑1 (MQ=255)
                                                                                      aTTTCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATGCTGGCGATTAATGTCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGcc                                                                <  7:228807/148‑1 (MQ=255)
                                                                                          cTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCagag                                                            <  8:281805/148‑1 (MQ=255)
                                                                                             tCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTCCTCAGGCCCACCGCGCGCCATTTAGCGCTTGTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCgg                                                         <  3:49967/148‑1 (MQ=255)
                                                                                               aGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGtt                                                       >  3:238313/1‑148 (MQ=255)
                                                                                                ggTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTa                                                     <  1:318008/149‑1 (MQ=255)
                                                                                                  tGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTAtt                                                   <  1:254241/149‑1 (MQ=255)
                                                                                                  tGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTAtt                                                   <  6:55885/149‑1 (MQ=255)
                                                                                                         cGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTGCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTtgt                                             >  6:56083/1‑148 (MQ=255)
                                                                                                            ccAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGat                                          >  7:146110/1‑148 (MQ=255)
                                                                                                                 gTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGt                                    <  6:56353/149‑1 (MQ=255)
                                                                                                                   tttCATGAGACCTATTCTGGCGATGAATGGCAGTCATTCCTTCAGATTAATCAGGCCAACCGCGCGCCATTTAGAGCGTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCGAGAGCGGTTTATGCTTGGTGATAATAGTg                                   >  6:254737/1‑148 (MQ=255)
                                                                                                                          cGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCa                            <  4:283227/148‑1 (MQ=255)
                                                                                                                          cGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCa                            <  8:17664/148‑1 (MQ=255)
                                                                                                                           gACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCTGTTTATTCTTTGTGATAATAGTGAAATCCAGa                          <  5:159205/149‑1 (MQ=255)
                                                                                                                             ccTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGAcc                        >  3:144975/1‑149 (MQ=255)
                                                                                                                             ccTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCAGTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCATAcc                        >  6:318338/1‑149 (MQ=255)
                                                                                                                                aTTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCAtt                                                                                                                >  6:331038/1‑58 (MQ=255)
                                                                                                                                aTTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCAtt                                                                                                                <  5:331038/58‑1 (MQ=255)
                                                                                                                                    tGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGa                 >  6:35539/1‑149 (MQ=255)
                                                                                                                                    tGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGa                 >  7:191172/1‑149 (MQ=255)
                                                                                                                                      tgaaggAATGCCTGCCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGAt                >  1:148643/6‑148 (MQ=255)
                                                                                                                                      tgaaggAATGCCTGCCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGAt                >  1:153329/6‑148 (MQ=255)
                                                                                                                                      gCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGAtt               <  6:38339/149‑1 (MQ=255)
                                                                                                                                        gATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTaa             <  1:136848/149‑1 (MQ=255)
                                                                                                                                         aggAATGCCTGCCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTaaa            <  2:148643/147‑1 (MQ=255)
                                                                                                                                         aggAATGCCTGCCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTaaa            <  2:153329/147‑1 (MQ=255)
                                                                                                                                         aTGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTaaa            <  2:277951/149‑1 (MQ=255)
                                                                                                                                              tGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGa                          >  7:328036/1‑130 (MQ=255)
                                                                                                                                              tGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGa                          <  8:328036/130‑1 (MQ=255)
                                                                                                                                                 caggcaTTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCAGATTAAAGGCACGc     >  6:111904/1‑148 (MQ=255)
                                                                                                                                                    gcaTTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTAAAGGCACGCTAc  >  1:219699/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTAAAGGCACGCTAC  >  NC_000913/1895341‑1895636

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