Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,328,452 C→T 5.3% R440C (CGT→TGT)  crfC → clamp‑binding sister replication fork colocalization protein, dynamin‑related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,328,4520CT 5.3% 257.0 / 3.5 75R440C (CGT→TGT) crfCclamp‑binding sister replication fork colocalization protein, dynamin‑related
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base T (2/2);  total (39/36)
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.08e-01

GGGAAGATTCGCTGTTCGCCCAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCG  >  NC_000913/4328315‑4328593
                                                                                                                                         |                                                                                                                                             
gggAAGATTCGCTGTTCGCCCAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACggg                                                                                                                                    >  1:147936/1‑149 (MQ=255)
gggAAGATTCGCTGTTCGCCCAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACggg                                                                                                                                    >  3:416586/1‑149 (MQ=255)
 ggAAGATTCGCTGTTCGCCCAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGt                                                                                                                                   >  2:255396/1‑149 (MQ=255)
 ggAAGATTCGCTGTTCGCCCAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGt                                                                                                                                   <  4:202908/149‑1 (MQ=255)
  gAAGATTCGCTGTTCGCCCAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGtt                                                                                                                                  >  1:515187/1‑149 (MQ=255)
   aaGATTCGCTGTTCGCCCAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTa                                                                                                                                 <  2:457337/149‑1 (MQ=255)
       ttCGCTGTTCGCCCAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACg                                                                                                                             >  2:531617/1‑149 (MQ=255)
        tCGCTGTTCGCCCAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGt                                                                                                                            <  2:504543/149‑1 (MQ=255)
           cTGTTCGCCCAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGc                                                                                                                         <  4:173212/149‑1 (MQ=255)
             gTTCGCCCAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCtt                                                                                                                       >  2:574498/1‑149 (MQ=255)
              ttCGCCCAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTtg                                                                                                                      >  2:187882/1‑149 (MQ=255)
                 gCCCAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGa                                                                                                                   >  2:221322/1‑149 (MQ=255)
                   ccAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAAc                                                                                                                 >  3:605211/1‑149 (MQ=255)
                     aGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTTGTGCGCACGGGTTAAACGTCGCTTGTGAACaa                                                                                                               >  3:582874/1‑149 (MQ=255)
                     aGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTTGTGCGCACGGGTTAAACGTCGCTTGTGAACaa                                                                                                               >  1:310761/1‑149 (MQ=255)
                         aaTTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGc                                                                                                           <  1:574498/149‑1 (MQ=255)
                          aTTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCg                                                                                                          <  3:177290/149‑1 (MQ=255)
                            tCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGc                                                                                                        >  2:483060/1‑149 (MQ=255)
                                       ttCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATccac                                                                                             >  3:104152/1‑149 (MQ=255)
                                       ttCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATccac                                                                                             <  3:157144/149‑1 (MQ=255)
                                       ttCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATccac                                                                                             >  3:382055/1‑149 (MQ=255)
                                         cATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATccacca                                                                                           <  2:84242/149‑1 (MQ=255)
                                             ccGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATc                                                                                       >  1:172963/1‑149 (MQ=255)
                                                cTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATTTCCACCAGATCgaa                                                                                    <  3:46518/149‑1 (MQ=255)
                                                cTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCgaa                                                                                    <  4:617880/149‑1 (MQ=255)
                                                cTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCgaa                                                                                    <  2:147936/149‑1 (MQ=255)
                                                cTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCgaa                                                                                    <  1:480396/149‑1 (MQ=255)
                                                 ttACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCgaag                                                                                   >  4:412680/1‑149 (MQ=255)
                                                    cGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGaaa                                                                                <  4:408504/149‑1 (MQ=255)
                                                         aCGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTg                                                                           <  2:515187/149‑1 (MQ=255)
                                                         aCGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTg                                                                           <  2:342218/149‑1 (MQ=255)
                                                          cGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGc                                                                          <  4:163892/149‑1 (MQ=255)
                                                           gCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCa                                                                         <  2:482987/149‑1 (MQ=255)
                                                           gCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCa                                                                         >  2:303800/1‑149 (MQ=255)
                                                           gCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCa                                                                         <  4:416080/149‑1 (MQ=255)
                                                             cTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGc                                                                       >  2:613746/1‑149 (MQ=255)
                                                                 ttGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTAtt                                                                   <  1:187882/149‑1 (MQ=255)
                                                                        ctctGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTc                                                            >  3:580733/1‑149 (MQ=255)
                                                                           tGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAAt                                                         <  4:580733/149‑1 (MQ=255)
                                                                             cGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCa                                                       <  3:496245/149‑1 (MQ=255)
                                                                               aTCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAgg                                                     <  3:615187/149‑1 (MQ=255)
                                                                                tCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGgc                                                    >  3:406424/1‑149 (MQ=255)
                                                                                   gCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCa                                                 <  4:104152/149‑1 (MQ=255)
                                                                                      gcgcATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGt                                              >  2:613496/1‑149 (MQ=255)
                                                                                        gcATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGa                                            <  2:419570/149‑1 (MQ=255)
                                                                                             aaCTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCgaagaa                                                                                 <  4:196984/107‑1 (MQ=255)
                                                                                             aaCTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCgaagaa                                                                                 >  3:196984/1‑107 (MQ=255)
                                                                                                tGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCgag                                    >  2:247045/1‑149 (MQ=255)
                                                                                                  ttGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGAt                                  >  1:137411/1‑149 (MQ=255)
                                                                                                    gAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTa                                >  1:522679/1‑149 (MQ=255)
                                                                                                      aTTACGCGCAGCAGGCGCGGGAATACCTGGATTTTTGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTaaa                              <  4:582874/149‑1 (MQ=255)
                                                                                                      aTTACGCGCAGCAGGCGCGGGAATACCTGGATTTTTGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTaaa                              <  2:310761/149‑1 (MQ=255)
                                                                                                      aTTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTaaa                              <  4:605211/149‑1 (MQ=255)
                                                                                                      aTTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTaaa                              <  4:290504/149‑1 (MQ=255)
                                                                                                          cgcgCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATg                          >  2:585144/1‑149 (MQ=255)
                                                                                                                gcagGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCg                    >  2:339153/1‑149 (MQ=255)
                                                                                                                 cagGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGa                   <  2:522679/149‑1 (MQ=255)
                                                                                                                  agGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAg                  <  1:483060/149‑1 (MQ=255)
                                                                                                                  agGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAg                  <  1:255396/149‑1 (MQ=255)
                                                                                                                   ggCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGc                 >  2:265147/1‑149 (MQ=255)
                                                                                                                   ggCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGc                 >  4:291305/1‑149 (MQ=255)
                                                                                                                    gcgcGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCt                >  3:305190/1‑149 (MQ=255)
                                                                                                                      gcgGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTgg              <  1:303800/149‑1 (MQ=255)
                                                                                                                          gAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCt          >  2:264376/1‑149 (MQ=255)
                                                                                                                          gAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCt          >  4:35661/1‑149 (MQ=255)
                                                                                                                          gAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCt          <  2:599885/149‑1 (MQ=255)
                                                                                                                            aTACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCATCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGa        >  4:396266/1‑149 (MQ=255)
                                                                                                                            aTACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTaaa                              >  2:76103/1‑127 (MQ=255)
                                                                                                                            aTACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTaaa                              <  1:76103/127‑1 (MQ=255)
                                                                                                                             tACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGAc       >  1:226366/1‑149 (MQ=255)
                                                                                                                                cTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTc    >  1:85530/1‑149 (MQ=255)
                                                                                                                                cTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTc    >  4:518161/1‑149 (MQ=255)
                                                                                                                                 tGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTcc   <  1:221322/149‑1 (MQ=255)
                                                                                                                                 tGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTcc   >  2:315777/1‑149 (MQ=255)
                                                                                                                                  ggATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTAGCCCTGACCTCCg  >  4:528296/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                             
GGGAAGATTCGCTGTTCGCCCAGCCAATTCAGGCGTTGCTTCATGCCGCTTACGCTAACGCCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCG  >  NC_000913/4328315‑4328593

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