Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,685,166 T→C 6.1% L303P (CTG→CCG)  tus → DNA replication terminus site‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,685,1660TC 6.1% 195.1 / 3.1 65L303P (CTG→CCG) tusDNA replication terminus site‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base C (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

CCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGCCCTTCTTTATGCGCTTTTTT  >  NC_000913/1685019‑1685313
                                                                                                                                                   |                                                                                                                                                   
ccTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGc                                                                                                                                                   >  1:1377712/1‑150 (MQ=255)
         caccacTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATg                                                                                                                                          <  2:638550/150‑1 (MQ=255)
           ccacTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGtt                                                                                                                                        >  2:1186844/1‑150 (MQ=255)
           ccacTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGtt                                                                                                                                        <  1:823902/150‑1 (MQ=255)
            cacTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTg                                                                                                                                       >  1:509366/1‑150 (MQ=255)
               tGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAg                                                                                                                                    <  2:168540/150‑1 (MQ=255)
                  ttGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGAtt                                                                                                                                 <  2:1322924/150‑1 (MQ=255)
                    gCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTaa                                                                                                                               >  2:1403159/1‑150 (MQ=255)
                     cACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAAc                                                                                                                              >  2:818059/1‑150 (MQ=255)
                     cACTGATTAATCGGGATAATGGCGCGGGCGTGACGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAAc                                                                                                                              >  2:1173534/1‑150 (MQ=255)
                             aaTCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTaa                                                                                                                               >  2:777030/1‑141 (MQ=255)
                             aaTCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTaa                                                                                                                               <  1:777030/141‑1 (MQ=255)
                                cGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTc                                                                                                                   >  1:817026/1‑150 (MQ=255)
                                cGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTc                                                                                                                   <  2:951636/150‑1 (MQ=255)
                                 gggATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCa                                                                                                                  <  2:664200/150‑1 (MQ=255)
                                       aTGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGc                                                                                                            >  1:307365/1‑150 (MQ=255)
                                        tGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGcc                                                                                                           <  1:340540/150‑1 (MQ=255)
                                        tGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGACGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGcc                                                                                                           >  2:566827/1‑150 (MQ=255)
                                         ggCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCg                                                                                                          >  2:646077/1‑150 (MQ=255)
                                           cgcgGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGAc                                                                                                        <  1:1248682/150‑1 (MQ=255)
                                            gcgGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCAGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTTCCGAcc                                                                                                       >  2:1394926/1‑150 (MQ=255)
                                               ggCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATc                                                                                                    >  2:256315/1‑150 (MQ=255)
                                                  gTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATTTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGt                                                                                                 <  2:1373552/150‑1 (MQ=255)
                                                     ccGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCt                                                                                              >  2:53346/1‑150 (MQ=255)
                                                        gACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGc                                                                                           <  1:38266/150‑1 (MQ=255)
                                                          cGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCGGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCg                                                                                         <  1:617200/150‑1 (MQ=255)
                                                              ggTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTAcc                                                                                     >  1:219118/1‑150 (MQ=255)
                                                                    ttgttAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTg                                                                               <  1:1398309/150‑1 (MQ=255)
                                                                    ttgttAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTg                                                                               <  1:744875/150‑1 (MQ=255)
                                                                      gttAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTc                                                                             >  1:472182/1‑150 (MQ=255)
                                                                       ttAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCa                                                                            <  1:256315/150‑1 (MQ=255)
                                                                          aaTTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAAc                                                                         <  1:1186844/150‑1 (MQ=255)
                                                                               cGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGc                                                                    >  2:1097777/1‑150 (MQ=255)
                                                                                 aTGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCtt                                                                  <  2:1305886/150‑1 (MQ=255)
                                                                                  tGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTc                                                                 <  2:307365/150‑1 (MQ=255)
                                                                                   gCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCCGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTg                                                                                                  >  1:131107/1‑116 (MQ=255)
                                                                                   gCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCCGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTg                                                                                                  <  2:131107/116‑1 (MQ=255)
                                                                                   gCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCCGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTg                                                                                                  <  2:136227/116‑1 (MQ=255)
                                                                                   cccGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCCGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTg                                                                                                  >  1:136227/2‑116 (MQ=255)
                                                                                    ccGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGc                                                               >  2:1120326/1‑150 (MQ=255)
                                                                                          aTGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGAt                                                         >  2:281427/1‑150 (MQ=255)
                                                                                             tgCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGTCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTAGCTAAGATAcc                                                      <  2:612279/150‑1 (MQ=255)
                                                                                                 gcaCCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACGATCTGTTCTGGCCGTAGCGAGCTGTGATACTCGGCTTCGCTAAGATACGCCAg                                                  <  2:509366/150‑1 (MQ=255)
                                                                                                 gcaCCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGtt                                                                                                <  1:1202247/104‑1 (MQ=255)
                                                                                                 gcaCCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGtt                                                                                                >  2:1202247/1‑104 (MQ=255)
                                                                                                      gTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCaa                                             >  2:188582/1‑150 (MQ=255)
                                                                                                         ataAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAggg                                          >  2:175875/1‑150 (MQ=255)
                                                                                                            aaCCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCg                                       >  2:1049122/1‑150 (MQ=255)
                                                                                                            aaCCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCg                                       >  1:858920/1‑150 (MQ=255)
                                                                                                             aCCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGc                                      <  1:254236/150‑1 (MQ=255)
                                                                                                             aCCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGc                                      >  1:444775/1‑150 (MQ=255)
                                                                                                                 cAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCt                                  <  2:858920/150‑1 (MQ=255)
                                                                                                                       cAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGTCTTTCATACTGCCGACAATCTGTTCTGGCCGTACCCAGCTGTAAAACTCGTCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAg                            <  1:1394926/150‑1 (MQ=255)
                                                                                                                        aGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAgg                           <  1:646077/150‑1 (MQ=255)
                                                                                                                               tCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTGAGccc                    <  2:1377712/150‑1 (MQ=255)
                                                                                                                                     gATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATCCTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGCCCTTCttt              <  1:365670/150‑1 (MQ=255)
                                                                                                                                     gATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGCCcttct                >  1:640076/1‑148 (MQ=255)
                                                                                                                                     gATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGCCcttct                <  2:640076/148‑1 (MQ=255)
                                                                                                                                       tCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGCCCTTCTTTAt            >  1:939707/1‑150 (MQ=255)
                                                                                                                                        cATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGCCCTTCTTTATg           <  2:950108/150‑1 (MQ=255)
                                                                                                                                           tCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGCCCTTCTTTATgcgc        >  2:920146/1‑150 (MQ=255)
                                                                                                                                            ccACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGCCCTTCTTTATGCGCt       >  2:361197/1‑150 (MQ=255)
                                                                                                                                            ccACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGACGCAGCTTTTAAGGTCAGCCCTTCTTTATGCGCt       >  1:768013/1‑150 (MQ=255)
                                                                                                                                            ccACGGCTGAAAATGTATGTTGAAGATTAACGCAAGGATTTCATACTGACGACCATCTGTTCTGGCCGTACCCAGATGTCAACCTCGCATTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAATGTCAGCCATTATTTATGctct       >  2:1021856/1‑150 (MQ=255)
                                                                                                                                                 gCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGCCCTTCTTTATGCGCtttttt  <  1:53346/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAATGGCGCGGGCGTGCCGGACGTTGGTGAGTTGTTAAATTACGATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGCCCTTCTTTATGCGCTTTTTT  >  NC_000913/1685019‑1685313

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