Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 594,950 T→G 5.3% 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 5.3% 396.0 / 12.5 114H165P (CAT→CCT) cusRresponse regulator in two‑component regulatory system with CusS
Reads supporting (aligned to +/- strand):  ref base T (53/55);  new base G (3/3);  total (56/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 114 total aligned reads displayed.

GGTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCC  >  NC_000913/594805‑595094
                                                                                                                                                 |                                                                                                                                                
ggTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGc                                                                                                                                               <  3:448713/149‑1 (MQ=255)
 gTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGc                                                                                                                                               <  1:914316/148‑1 (MQ=255)
 gTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCg                                                                                                                                              >  3:864292/1‑149 (MQ=255)
   cTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaa                                                                                                                                            <  3:4626/149‑1 (MQ=255)
   cTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCga                                                                                                                                             <  2:126842/148‑1 (MQ=255)
    tGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaag                                                                                                                                           <  4:463415/149‑1 (MQ=255)
     gAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaag                                                                                                                                           <  2:966921/148‑1 (MQ=255)
       atgaGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGACGCTTCACCGCCACATCAATAGCATTGGAGTTGCGGTCAAAGTTCCTGTGCCATGCCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaagga                                                                                                                                         <  1:298016/145‑1 (MQ=255)
       aTTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaaggaa                                                                                                                                        <  3:355582/149‑1 (MQ=255)
       aTTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaaggaa                                                                                                                                        >  3:793431/1‑149 (MQ=255)
        ttAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaaggaag                                                                                                                                       <  4:249052/149‑1 (MQ=255)
         tAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaaggaag                                                                                                                                       >  1:156176/1‑148 (MQ=255)
          aGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGgaagaa                                                                                                                                     <  4:108717/149‑1 (MQ=255)
          aGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGgaaga                                                                                                                                      <  2:213738/148‑1 (MQ=255)
           gCTTCGGCTCAAAGTCGTTTTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATTAATAGCATTGGAATCGCTGACAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAAc                                                                                                                                    >  4:661263/1‑149 (MQ=255)
           gCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGgaagaa                                                                                                                                     <  2:488431/148‑1 (MQ=255)
                  cTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaag                                                                                                                                           <  1:595785/135‑1 (MQ=255)
                   tCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagcag                                                                                                                            >  3:330805/1‑149 (MQ=255)
                   tCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagca                                                                                                                             <  1:526517/148‑1 (MQ=255)
                    cAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACt                                                                                                                                   <  1:467554/141‑1 (MQ=255)
                    cAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACt                                                                                                                                   >  2:467554/1‑141 (MQ=255)
                    cAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagcag                                                                                                                            >  2:750031/1‑148 (MQ=255)
                    cAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagcag                                                                                                                            <  1:273942/148‑1 (MQ=255)
                    cAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagcag                                                                                                                            <  1:628793/148‑1 (MQ=255)
                     aaaGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCagag                                                                                                                          >  3:131419/1‑149 (MQ=255)
                       aGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGt                                                                                                                         >  2:663101/1‑148 (MQ=255)
                          cGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAGTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTaaa                                                                                                                      >  1:873409/1‑148 (MQ=255)
                            ttGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAActc                                                                                                                   >  3:247486/1‑149 (MQ=255)
                             tGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAActct                                                                                                                  <  4:864292/149‑1 (MQ=255)
                             tGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAActc                                                                                                                   >  1:766102/1‑148 (MQ=255)
                               tCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCtt                                                                                                                 >  1:848038/1‑148 (MQ=255)
                                cGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCttt                                                                                                                >  1:709510/1‑148 (MQ=255)
                                cGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCttt                                                                                                                <  2:716439/148‑1 (MQ=255)
                                  aTTTTGCCGCGCAGCCGCTTCACCGCCACATCAATGGCATTGGTATCGCTGTCAAACTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCGGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTAc                                                                                                              >  2:190361/1‑148 (MQ=255)
                                   ttttGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACt                                                                                                             <  1:484209/148‑1 (MQ=255)
                                   ttttGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTg                                                                                                            <  4:431000/149‑1 (MQ=255)
                                         cgcgCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCaa                                                                                                       <  2:848038/148‑1 (MQ=255)
                                           cgcAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGt                                                                                                    >  4:341295/1‑149 (MQ=255)
                                            gcAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCaaa                                                                                                      >  4:71096/1‑146 (MQ=255)
                                            gcAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCaaa                                                                                                      <  3:71096/146‑1 (MQ=255)
                                              aGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGAt                                                                                                 <  4:716927/149‑1 (MQ=255)
                                               gCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGAt                                                                                                 <  2:234603/148‑1 (MQ=255)
                                               gCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATg                                                                                                <  4:131419/149‑1 (MQ=255)
                                                  gCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATgcgc                                                                                             >  3:504826/1‑149 (MQ=255)
                                                   cTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATgcgc                                                                                             >  1:83620/1‑148 (MQ=255)
                                                     tCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGt                                                                                           >  1:121394/1‑148 (MQ=255)
                                                            cacaTCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCCGATGGCGAAGGAAGAACTACAGCAGAGTAAACTCTTTACTGTTCAAAGTGATGCGCGTGCCGCTg                                                                                    >  2:859049/1‑148 (MQ=255)
                                                              caTCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCgg                                                                                 <  3:487420/149‑1 (MQ=255)
                                                                  aaTAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGa                                                                             <  4:302290/149‑1 (MQ=255)
                                                                   aTAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTAc                                                                                                              >  2:120918/1‑115 (MQ=255)
                                                                   aTAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTAc                                                                                                              <  1:120918/115‑1 (MQ=255)
                                                                     aGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACt                                                                           >  1:89892/1‑148 (MQ=255)
                                                                        aTTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGc                                                                       <  4:502767/149‑1 (MQ=255)
                                                                        aTTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGc                                                                       <  4:538069/149‑1 (MQ=255)
                                                                         ttGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGc                                                                       <  2:775400/148‑1 (MQ=255)
                                                                          tGGTATCGCTTTCAAAATTCATATCCCATACCTGCGTGGCGATAAGCGAGCGGGGCAGCACTTCGACCTGATGGCGAAGGAAGAACTCCAGCCGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGAATTTGCg                                                                      >  2:682018/1‑148 (MQ=255)
                                                                            gTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATgg                                                                                                                                                <  3:39168/72‑1 (MQ=255)
                                                                             tATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCagag                                                                                                                          >  1:529919/1‑93 (MQ=255)
                                                                             tATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCagag                                                                                                                          <  2:529919/93‑1 (MQ=255)
                                                                             tATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGACt                                                                   >  2:590130/1‑148 (MQ=255)
                                                                              aTCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGAt                                                                                                 >  4:431039/1‑117 (MQ=255)
                                                                              aTCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGAt                                                                                                 <  3:431039/117‑1 (MQ=255)
                                                                              aTCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTg                                                                  <  2:497066/148‑1 (MQ=255)
                                                                                    gTCAAGATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACgaga                                                            >  2:684290/1‑148 (MQ=255)
                                                                                     tCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGGGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATc                                                          <  4:550516/149‑1 (MQ=255)
                                                                                     tCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGAt                                                           >  2:876664/1‑148 (MQ=255)
                                                                                      cAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATc                                                          >  2:321503/1‑148 (MQ=255)
                                                                                            tCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCAt                                                    >  2:963982/1‑148 (MQ=255)
                                                                                               tGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCaa                                                 >  2:439649/1‑148 (MQ=255)
                                                                                               tGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCaa                                                 <  2:121394/148‑1 (MQ=255)
                                                                                                gTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCaaa                                                >  2:493430/1‑148 (MQ=255)
                                                                                                gTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAAt                                               >  4:641470/1‑149 (MQ=255)
                                                                                                 tCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATc                                              >  4:37352/1‑149 (MQ=255)
                                                                                                 tCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGAGAAGGAAGAACTCCAGCAGAGTAAGCTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAAt                                               <  2:709510/148‑1 (MQ=255)
                                                                                                   ccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAActct                                                                                                                  <  3:152001/79‑1 (MQ=255)
                                                                                                   ccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAActct                                                                                                                  >  4:152001/1‑79 (MQ=255)
                                                                                                   ccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGa                                                        <  3:360148/137‑1 (MQ=255)
                                                                                                   ccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGa                                                        >  4:360148/1‑137 (MQ=255)
                                                                                                   ccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCgg                                            >  3:702839/1‑149 (MQ=255)
                                                                                                   ccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTCCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCgg                                            >  4:501763/1‑149 (MQ=255)
                                                                                                      tACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCGCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGGCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCa                                          >  2:15620/1‑148 (MQ=255)
                                                                                                      tACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCa                                          <  1:904999/148‑1 (MQ=255)
                                                                                                      tACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCa                                          <  1:765016/148‑1 (MQ=255)
                                                                                                         cTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAAcc                                       <  1:420568/148‑1 (MQ=255)
                                                                                                          tGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCt                                      >  1:377775/1‑148 (MQ=255)
                                                                                                          tGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGGGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCt                                      <  1:663101/148‑1 (MQ=255)
                                                                                                                     aaGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACtt                           <  1:321503/148‑1 (MQ=255)
                                                                                                                      agcgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTc                         <  3:341295/149‑1 (MQ=255)
                                                                                                                      agcgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTc                         <  4:389450/149‑1 (MQ=255)
                                                                                                                        cgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCg                        <  1:321226/148‑1 (MQ=255)
                                                                                                                          agcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGTTGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGAt                      <  2:873409/148‑1 (MQ=255)
                                                                                                                               ggCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAAt                   >  4:877738/1‑146 (MQ=255)
                                                                                                                               ggCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAAt                   <  3:877738/146‑1 (MQ=255)
                                                                                                                                    tccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACg                                                               >  1:808918/3‑97 (MQ=255)
                                                                                                                                    tccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACg                                                               <  2:808918/95‑1 (MQ=255)
                                                                                                                                       ttCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCgccgcc         >  2:112683/1‑148 (MQ=255)
                                                                                                                                       ttCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGccc        >  4:641720/1‑149 (MQ=255)
                                                                                                                                          gCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCcgc     >  3:411399/1‑149 (MQ=255)
                                                                                                                                            ccTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCcgcg    >  2:320630/1‑148 (MQ=255)
                                                                                                                                              tGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGcgcc  >  2:629302/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCC  >  NC_000913/594805‑595094

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