Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 594,950 T→G 6.7% H165P (CAT→CCT)  cusR ← response regulator in two‑component regulatory system with CusS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913594,9500TG 6.7% 395.8 / 18.1 119H165P (CAT→CCT) cusRresponse regulator in two‑component regulatory system with CusS
Reads supporting (aligned to +/- strand):  ref base T (54/57);  new base G (6/2);  total (60/59)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

Only 100 of 119 total aligned reads displayed.

CTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCG  >  NC_000913/594808‑595095
                                                                                                                                               |                                                                                                                                                 
cTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaa                                                                                                                                             <  3:665551/149‑1 (MQ=255)
 tGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaag                                                                                                                                            >  1:847651/1‑149 (MQ=255)
 tGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaag                                                                                                                                            <  4:344372/149‑1 (MQ=255)
 tGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaag                                                                                                                                            <  4:360380/149‑1 (MQ=255)
 tGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaag                                                                                                                                            <  3:85958/149‑1 (MQ=255)
   aaTTAGCTTCGGCTCAAAGTCGTTGTCGATTTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaagg                                                                                                                                           <  1:44784/149‑1 (MQ=255)
    aTTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaaggaa                                                                                                                                         <  2:685384/149‑1 (MQ=255)
    aTTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaaggaa                                                                                                                                         <  2:841868/149‑1 (MQ=255)
     ttAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaaggaag                                                                                                                                        >  3:146677/1‑149 (MQ=255)
       aGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGgaagaa                                                                                                                                      <  4:735404/149‑1 (MQ=255)
       aGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGgaagaa                                                                                                                                      <  1:625677/149‑1 (MQ=255)
        gCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaaggaag                                                                                                                                        >  4:342755/1‑146 (MQ=255)
        gCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaaggaag                                                                                                                                        <  3:342755/146‑1 (MQ=255)
         cTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACt                                                                                                                                    >  4:803029/1‑149 (MQ=255)
           tCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTcc                                                                                                                                  >  2:399261/1‑149 (MQ=255)
               cTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagca                                                                                                                              >  3:358846/1‑149 (MQ=255)
               cTCAAAGGCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagca                                                                                                                              >  2:828604/1‑149 (MQ=255)
                tCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagcag                                                                                                                             >  1:154764/1‑149 (MQ=255)
                 cAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCaga                                                                                                                            <  2:202899/149‑1 (MQ=255)
                 cAAAGGCGTTGTCGATTTGGCCGCGCAGCCGCTTTACCGCCACATCAATAGCCTTGGTATCGCTCTCAAAATTCATTTCCCATACGTGCGAGGCAATAAGTGAGTGGTGCAGCACTTCTCCCTGATGGCGAAGGAAGAACTCCAGCaga                                                                                                                            >  2:41181/1‑149 (MQ=255)
                  aaaGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCGATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCagag                                                                                                                           >  4:540388/1‑149 (MQ=255)
                  aaaGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCagag                                                                                                                           >  4:903313/1‑149 (MQ=255)
                     gTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTaa                                                                                                                        >  4:479215/1‑149 (MQ=255)
                      tCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTaaa                                                                                                                       <  3:513711/149‑1 (MQ=255)
                      tCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTaaa                                                                                                                       <  3:753018/149‑1 (MQ=255)
                             cGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTa                                                                                                                <  4:462474/149‑1 (MQ=255)
                             cGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTa                                                                                                                <  3:714111/149‑1 (MQ=255)
                              gATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTAc                                                                                                               <  3:819204/149‑1 (MQ=255)
                                 ttttGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTg                                                                                                             >  3:425419/1‑149 (MQ=255)
                                       cgcgCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCaaa                                                                                                       <  3:405838/149‑1 (MQ=255)
                                       cgcgCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCaaa                                                                                                       <  4:146677/149‑1 (MQ=255)
                                         cgcAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGt                                                                                                     <  1:378323/149‑1 (MQ=255)
                                          gcAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTg                                                                                                    <  2:458892/149‑1 (MQ=255)
                                                gCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATgcgc                                                                                              >  4:767456/1‑149 (MQ=255)
                                                gCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATgcgc                                                                                              >  1:274472/1‑149 (MQ=255)
                                                 cTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGcgcg                                                                                             <  2:707264/149‑1 (MQ=255)
                                                  ttCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGt                                                                                            >  2:526172/1‑149 (MQ=255)
                                                   tCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTg                                                                                           <  1:617305/149‑1 (MQ=255)
                                                       cGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGc                                                                                       <  1:84728/149‑1 (MQ=255)
                                                         ccACATCAATAGCATTGGTATCGCTGTCAAGATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTg                                                                                     >  2:211843/1‑149 (MQ=255)
                                                          cacaTCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGc                                                                                    <  4:217329/149‑1 (MQ=255)
                                                            caTCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCgg                                                                                  <  3:619104/149‑1 (MQ=255)
                                                            caTCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCgg                                                                                  <  2:919223/149‑1 (MQ=255)
                                                              tCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGt                                                                                <  2:237666/149‑1 (MQ=255)
                                                                      aTTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGc                                                                        >  1:738591/1‑149 (MQ=255)
                                                                      aTTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGc                                                                        >  3:650138/1‑149 (MQ=255)
                                                                       ttGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCg                                                                       >  3:555545/1‑149 (MQ=255)
                                                                            aTCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGa                                                                  <  3:426453/149‑1 (MQ=255)
                                                                             tCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTa                                                                                                                <  3:348170/102‑1 (MQ=255)
                                                                             tCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTa                                                                                                                >  4:348170/1‑102 (MQ=255)
                                                                               gCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACga                                                               <  4:561910/149‑1 (MQ=255)
                                                                               gCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACga                                                               <  4:55516/149‑1 (MQ=255)
                                                                                   tCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCTAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACAGCAGCAGAGTAAACTCTTTACTTGTCAAAGTGATGCGCATGCCGCTGCGGGTGACTTTGCGGCTGACGAGATc                                                           >  4:654648/1‑149 (MQ=255)
                                                                                   tCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATc                                                           >  2:874508/1‑149 (MQ=255)
                                                                                     aaaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACgaga                                                             >  4:847442/1‑145 (MQ=255)
                                                                                     aaaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACgaga                                                             <  3:847442/145‑1 (MQ=255)
                                                                                     aaaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGa                                                         <  4:450009/149‑1 (MQ=255)
                                                                                     aaaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGa                                                         >  1:397157/1‑149 (MQ=255)
                                                                                     aaaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGa                                                         >  2:434775/1‑149 (MQ=255)
                                                                                     aaaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGa                                                         >  2:858405/1‑149 (MQ=255)
                                                                                     aaaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGa                                                         <  1:520803/149‑1 (MQ=255)
                                                                                     aaaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGa                                                         >  3:904223/1‑149 (MQ=255)
                                                                                      aaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGAc                                                        >  2:387878/1‑149 (MQ=255)
                                                                                      aaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGAc                                                        <  4:358846/149‑1 (MQ=255)
                                                                                       aaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGAcc                                                       >  3:308609/1‑149 (MQ=255)
                                                                                       aaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGAcc                                                       <  3:159856/149‑1 (MQ=255)
                                                                                          tCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATc                                                    >  1:777769/1‑149 (MQ=255)
                                                                                                cccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCg                                              >  4:456690/1‑149 (MQ=255)
                                                                                                cccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCg                                              <  1:399261/149‑1 (MQ=255)
                                                                                                 ccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCgg                                             >  2:738960/1‑149 (MQ=255)
                                                                                                 ccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCgg                                             >  2:118365/1‑149 (MQ=255)
                                                                                                 ccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCgg                                             >  2:546251/1‑149 (MQ=255)
                                                                                                  cATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGc                                            >  1:337486/1‑149 (MQ=255)
                                                                                                            aGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAAc                                  >  1:853237/1‑149 (MQ=255)
                                                                                                             ggCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACt                                 <  2:154764/149‑1 (MQ=255)
                                                                                                               cGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGa                               <  4:687287/149‑1 (MQ=255)
                                                                                                                gATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAGGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGAc                              >  3:603673/1‑149 (MQ=255)
                                                                                                                gATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGAc                              >  2:915637/1‑149 (MQ=255)
                                                                                                                gATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGAc                              >  1:89967/1‑149 (MQ=255)
                                                                                                                  tAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACtt                            >  3:111423/1‑149 (MQ=255)
                                                                                                                    agcgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTc                          <  4:472983/149‑1 (MQ=255)
                                                                                                                     gcgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCg                         >  4:216710/1‑149 (MQ=255)
                                                                                                                      cgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGa                        <  4:650138/149‑1 (MQ=255)
                                                                                                                      cgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGa                        <  1:434775/149‑1 (MQ=255)
                                                                                                                       gagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGAt                       >  4:254034/1‑149 (MQ=255)
                                                                                                                        agcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATa                      <  4:550048/149‑1 (MQ=255)
                                                                                                                         gcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATaa                     <  4:603673/149‑1 (MQ=255)
                                                                                                                              ttcttccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAcc                >  1:151908/7‑149 (MQ=255)
                                                                                                                              gcagcaCTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAcc                >  3:9162/1‑149 (MQ=255)
                                                                                                                                cttccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAccgc              >  1:74109/5‑149 (MQ=255)
                                                                                                                                 ttccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACt                             >  1:232844/4‑133 (MQ=255)
                                                                                                                                  tcgttCGCCATCAGGGCGAAGGAAGAACTCCAACAGAGTAGACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGAGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAccgccg            >  2:518650/4‑149 (MQ=255)
                                                                                                                                  tccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAccgccg            >  3:453709/3‑149 (MQ=255)
                                                                                                                                  tccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAccgccg            >  3:343532/3‑149 (MQ=255)
                                                                                                                                      tCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGcccc        <  4:861166/149‑1 (MQ=255)
                                                                                                                                          ccTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCg                         <  2:850488/128‑1 (MQ=255)
                                                                                                                                          ccTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCg                         >  1:850488/1‑128 (MQ=255)
                                                                                                                                          ccTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCgcgc    >  2:171434/1‑149 (MQ=255)
                                                                                                                                           cTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGcgcc   <  4:827591/149‑1 (MQ=255)
                                                                                                                                            tGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGcgccg  >  4:495238/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCG  >  NC_000913/594808‑595095

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