Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,895,487 A→T 6.7% Q228L (CAG→CTG)  pabB → aminodeoxychorismate synthase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,895,4870AT 6.7% 188.9 / 3.6 60Q228L (CAG→CTG) pabBaminodeoxychorismate synthase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (24/32);  new base T (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.18e-01

CGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTAAAGGCACGC  >  NC_000913/1895340‑1895633
                                                                                                                                                   |                                                                                                                                                  
cGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggcag                                                                                                                                                   <  8:44956/149‑1 (MQ=255)
cGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggcag                                                                                                                                                   <  7:287961/149‑1 (MQ=255)
 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                                                                                                                                    >  4:158164/1‑148 (MQ=255)
                gCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCtt                                                                                                                                    >  2:231723/1‑148 (MQ=255)
                gCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCtt                                                                                                                                    >  6:82898/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                                                <  6:55885/149‑1 (MQ=255)
                                                                                                   tGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTAtt                                                <  1:254241/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                                                                                                             <  5:331038/58‑1 (MQ=255)
                                                                                                                                 aTTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCAtt                                                                                                             >  6:331038/1‑58 (MQ=255)
                                                                                                                                     tGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGa              >  7:191172/1‑149 (MQ=255)
                                                                                                                                     tGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGa              >  6:35539/1‑149 (MQ=255)
                                                                                                                                       tgaaggAATGCCTGCCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGAt             >  1:153329/6‑148 (MQ=255)
                                                                                                                                       tgaaggAATGCCTGCCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGAt             >  1:148643/6‑148 (MQ=255)
                                                                                                                                       gCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGAtt            <  6:38339/149‑1 (MQ=255)
                                                                                                                                         gATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTaa          <  1:136848/149‑1 (MQ=255)
                                                                                                                                          aggAATGCCTGCCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTaaa         <  2:153329/147‑1 (MQ=255)
                                                                                                                                          aggAATGCCTGCCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTaaa         <  2:148643/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)
                                                                                                                                                   |                                                                                                                                                  
CGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTAAAGGCACGC  >  NC_000913/1895340‑1895633

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: