Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,627,828 C→A 8.6% P172T (CCG→ACG)  nadR → trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,627,8280CA 8.6% 287.3 / 5.4 104P172T (CCG→ACG) nadRtrifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (47/46);  new base A (2/9);  total (49/55)
Fisher's exact test for biased strand distribution p-value = 5.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

CGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAA  >  NC_000913/4627681‑4627974
                                                                                                                                                   |                                                                                                                                                  
cGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCg                                                                                                                                                  >  2:1426955/1‑150 (MQ=255)
   cTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGTTCCAGCCGgat                                                                                                                                               >  1:733938/1‑150 (MQ=255)
      cGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctg                                                                                                                                            >  1:1689013/1‑150 (MQ=255)
        ttGGTTATTGCAAACTTTTAAATATCAGAAAAATCTTCGCATTCATGCTTTCAACGAAGAGGGCATCGAGCCGTATCCGCACGGCTGGGATGTGTGGAGAAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgat                                                                                                                                          <  1:744484/150‑1 (MQ=255)
        ttGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgat                                                                                                                                          <  2:744634/150‑1 (MQ=255)
        ttGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgat                                                                                                                                          <  1:1306393/150‑1 (MQ=255)
        ttGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgat                                                                                                                                          <  1:980280/150‑1 (MQ=255)
        ttGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGATGGGATGTCTGGAGCAACGGCATCAAAAAGATTATGGCTGAAAAAGAGATACAGACCgatctgat                                                                                                                                          >  1:24161/1‑150 (MQ=255)
         tGGTTATTGCTAACTTTTAAATATCAGAAAAATATTCGCAGTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCAAGGCTGGGATTTGAGGAGCAACGGCATCAAAATGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatc                                                                                                                                         <  1:98584/150‑1 (MQ=255)
         tGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatc                                                                                                                                         <  2:795731/150‑1 (MQ=255)
         tGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatc                                                                                                                                         <  1:499277/150‑1 (MQ=255)
         tGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatc                                                                                                                                         <  1:1582844/150‑1 (MQ=255)
         tGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatc                                                                                                                                         <  1:788456/150‑1 (MQ=255)
         tGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatc                                                                                                                                         <  2:1607048/150‑1 (MQ=255)
         gggctATTGCAAACTTTTAAATATCAGAAAATTATTCGCACTCATGCTCTCAACGAAGACGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatc                                                                                                                                         <  2:333514/146‑1 (MQ=255)
          ggTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatct                                                                                                                                        <  2:1466776/150‑1 (MQ=255)
          ggTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatct                                                                                                                                        <  1:523121/150‑1 (MQ=255)
          ggTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAGGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatct                                                                                                                                        <  1:529691/150‑1 (MQ=255)
              aTTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTacac                                                                                                                                    >  1:722173/1‑150 (MQ=255)
              aTTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTacac                                                                                                                                    >  1:684455/1‑150 (MQ=255)
               ttGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTacaca                                                                                                                                   >  1:199752/1‑149 (MQ=255)
               ttGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACAcc                                                                                                                                   <  1:524179/150‑1 (MQ=255)
                 gCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTc                                                                                                                                 >  2:448164/1‑150 (MQ=255)
                   aaaCTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCgg                                                                                                                               <  1:1121298/150‑1 (MQ=255)
                    aaCTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGga                                                                                                                              <  1:421373/150‑1 (MQ=255)
                         ttAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATCCAGCCGGATCTGATCTACACCTCGGaagaag                                                                                                                         <  1:631396/150‑1 (MQ=255)
                            aaTATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTTTGGATTTTAAAGGGATGCAGACGGATGTGATCTACACCTCGGAAGAAGCCg                                                                                                                      <  2:1487061/150‑1 (MQ=255)
                            aaTATCAGAAAAATAATCGCATTCATGCTTTCAACGAAGAGGGCCTGGAGCAGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATAAAAAAGTTTTTTTTTTTTAATGTTACGGAGACGTATGTGATCTACACCTCGGAAGAAGCCg                                                                                                                      <  2:260383/150‑1 (MQ=255)
                             atatCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATATACACCTCGGAAGAAGCCGa                                                                                                                     >  1:1055423/1‑150 (MQ=255)
                             atatCAGAAAAATATACGCATTCAAGCTATCAACGAAGAGCTCAGCGAGCCGTATCCTCACTGCTGGTATGTGTGTAGCAACGCCATCAATAAGTTTTTTTTTTTTAATGTTATTGAGACGTATGAGATCTACACCTCGGAAGAAGACGa                                                                                                                     <  2:92314/150‑1 (MQ=255)
                              tatCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGAt                                                                                                                    >  1:497503/1‑150 (MQ=255)
                               agcAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGACGTATCCGCACGGCAGGGATGTGTGGAACAACGGCATCAAAAAGTTTATAGCTGAAAAATGTATCCATCCGGATCTGATCTACACATATTAAGAAGcctaag                                                                                                                   >  2:122414/3‑146 (MQ=255)
                                tCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGGGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATgc                                                                                                                  <  1:1404357/150‑1 (MQ=255)
                                  aGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGAGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGACGGATCTGATCAACACCTCGGAAGAAGCCGATgctc                                                                                                                >  1:267458/1‑148 (MQ=255)
                                   gAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTTTTTTTTTTAATGTTAGGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGcc                                                                                                               <  2:796567/150‑1 (MQ=255)
                                   gAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGcc                                                                                                               >  2:1550829/1‑150 (MQ=255)
                                     aaaaTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCcac                                                                                                             <  1:1537941/150‑1 (MQ=255)
                                      aaaTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTATATGGCTGAAAAAGGGATCCCGCCGGCTCTGATCTACACCTCGTAAGAAGCCGAAGCGCcaca                                                                                                            >  1:441004/1‑150 (MQ=255)
                                       aaTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAg                                                                                                           <  1:1377782/150‑1 (MQ=255)
                                             cGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATg                                                                                                     <  1:283369/150‑1 (MQ=255)
                                             cGCATTAATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATg                                                                                                     >  2:1633507/1‑150 (MQ=255)
                                                 ttCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTTCTGAAAAAGTTATCCAGACGGATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAAc                                                                                                 <  2:1478294/150‑1 (MQ=255)
                                                    aTGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGCTGCGCCACAGTATATGGAACATc                                                                                              >  2:1435360/1‑150 (MQ=255)
                                                         ttCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTgggg                                                                                         <  1:1213863/150‑1 (MQ=255)
                                                           cAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAt                                                                                       >  2:650681/1‑150 (MQ=255)
                                                             aCGCAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGGGTGGAGCAACGGCATCAAAAAGTTTATGTCTGTAAAAGTGATCCAGCCGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCg                                                                                     <  2:475228/150‑1 (MQ=255)
                                                             aCGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTTTTTTTTTTAATGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCg                                                                                     <  2:749854/150‑1 (MQ=255)
                                                             aCGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTATATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGCTGCGCCACAGTATATGGAACATCTGGGGATCg                                                                                     >  1:39096/1‑150 (MQ=255)
                                                               gAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCgag                                                                                   <  2:1658479/150‑1 (MQ=255)
                                                                 agagGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGAc                                                                                 <  1:443715/150‑1 (MQ=255)
                                                                 agagGGCATGGAGCCGTATCCGCACGGCTGGGAGGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGAc                                                                                 <  1:7411/150‑1 (MQ=255)
                                                                    gggCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGGATATGGAACATCTGGGGATCGAGACggg                                                                              >  2:1652933/1‑149 (MQ=255)
                                                                     ggCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTg                                                                             >  2:1088119/1‑150 (MQ=255)
                                                                      gCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGc                                                                            >  2:449057/1‑150 (MQ=255)
                                                                          ggAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGTCATCAAAAAGTTTTTTTTTTTAAATGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGt                                                                        <  2:813564/150‑1 (MQ=255)
                                                                            aGCAGTATCCGCACGGCTGGGATGTGTGGAGCAACGTCATCAAAAAGTTTTTTTTTTTTAATGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCg                                                                      <  2:24161/150‑1 (MQ=255)
                                                                                gtttgcGCCCGGCTGGGAGGTGTGGAGCACCGGCATCATAATGTTTTTGGTTTTTTAAGTTATGGCGCCGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATcc                                                                  <  2:875923/145‑1 (MQ=255)
                                                                                gaaTCCGCTCGGCGGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATcc                                                                  <  1:450093/148‑1 (MQ=255)
                                                                                gTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATcc                                                                  <  1:518135/150‑1 (MQ=255)
                                                                                gTATCCGCACGGCTGGGAGGTGTGGAGCAACGGCATCAAAAAGTTTATTGCTGAAAATGGTATCCAGCCTGATGTGATATACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATcc                                                                  <  2:858714/150‑1 (MQ=255)
                                                                                    gcGCACGGCTTGGATGTGTGGAGCAAGGTCATCAAATAGTTTATGTTTTATAATGTTATGGAGACGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGaaa                                                              <  2:1077986/149‑1 (MQ=255)
                                                                                    ccGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGaaa                                                              >  2:1003730/1‑150 (MQ=255)
                                                                                    ccGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGaaa                                                              >  1:29587/1‑150 (MQ=255)
                                                                                            cTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCtt                                                      >  2:956411/1‑150 (MQ=255)
                                                                                                gATGTGTGGAGCAACGGCATCAAAAAGTTTATGGATGATAAAGTGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATg                                                  <  2:500236/150‑1 (MQ=255)
                                                                                                       ggAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGCCGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCa                                           >  1:874547/1‑150 (MQ=255)
                                                                                                         aGCAACGTCATCAAATATTTTATGGATGAAAAAGGTATCCAGCCGGATCTGATCTACACCTCGGAAGATGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGc                                         <  2:1700864/150‑1 (MQ=255)
                                                                                                            aaCGGCATCAAAAAGTTTATGGCTTAAAAAGTGATCCAGCAGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGt                                      <  2:1635662/150‑1 (MQ=255)
                                                                                                            aaCGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGACACAGTATATGGAACATCTGGGGATCGAGAAGGTGCTGGTCGATCCGAAACGTACCTATATGAGTATCAGCGGt                                      >  1:340077/1‑150 (MQ=255)
                                                                                                                 cATCAAAAAGTTTATGTCTGAAAAAGGTATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCa                                 <  2:1232996/150‑1 (MQ=255)
                                                                                                                  aTCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAg                                >  2:437300/1‑150 (MQ=255)
                                                                                                                     aaaaaTTTTATGGTTTATAATGGTATGGAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATc                             <  2:684455/150‑1 (MQ=255)
                                                                                                                     aaaaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTTATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATc                             <  1:1426955/150‑1 (MQ=255)
                                                                                                                     aaaaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGACGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATc                             >  1:873867/1‑150 (MQ=255)
                                                                                                                     aaaaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATc                             <  1:1693740/150‑1 (MQ=255)
                                                                                                                       aaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCcg                           <  1:1453245/150‑1 (MQ=255)
                                                                                                                        aaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCcgc                          <  1:708029/150‑1 (MQ=255)
                                                                                                                         aGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGAACCGAAACGTACCTTGATGAGTATCAGCGGTGCGCAGATCcgcg                         >  1:299819/1‑150 (MQ=255)
                                                                                                                                ggCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAGCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGACTATCAGCGGTGCGCAGATCCGCGAAAAccc                  <  1:1633507/150‑1 (MQ=255)
                                                                                                                                  cTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCATCGGTGCGCAGATCCGCGAAAACCCGt                >  2:32701/1‑150 (MQ=255)
                                                                                                                                  cTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGt                >  1:621450/1‑150 (MQ=255)
                                                                                                                                  cTGAAAAAGGGAACCAGCAGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCAAAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCATTATGAGAATCAGCGGTGCGCAGCTCCGCGAAAACCCGt                >  1:869263/1‑150 (MQ=255)
                                                                                                                                   tGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAGCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGAGCCGCGAAAACCCGtt               <  1:1012835/150‑1 (MQ=255)
                                                                                                                                     aaaaaGGGATCCAGCCGGATCTGATCTACACCTCGGAAGCAGCCGATGCGCCACAGTATATGGACCATCTGGGGATCGAGACGGTGCTGGTCGATCCGAACCGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTcc             >  1:1143611/1‑150 (MQ=255)
                                                                                                                                     aaaaaGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTcc             >  1:1324498/1‑150 (MQ=255)
                                                                                                                                     aaaaaGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTcc             >  1:1272890/1‑150 (MQ=255)
                                                                                                                                      aaaaGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCg            >  1:1523201/1‑150 (MQ=255)
                                                                                                                                      aaaaGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCg            >  2:786426/1‑150 (MQ=255)
                                                                                                                                       gatcagATCCGGCTGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATCGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGc           >  1:713668/6‑150 (MQ=255)
                                                                                                                                       aaaGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACATTATATGGAACATCTGGGGATCGAGCCGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTcaga           >  2:684353/1‑147 (MQ=255)
                                                                                                                                          gggATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTAc        >  1:947880/1‑150 (MQ=255)
                                                                                                                                           ggATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACt       >  1:1530719/1‑150 (MQ=255)
                                                                                                                                           ggATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACt       >  1:1466374/1‑150 (MQ=255)
                                                                                                                                           ggATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACt       >  1:144844/1‑150 (MQ=255)
                                                                                                                                           ggATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCAACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGCTCGATCCGAAACGTACCTTAAAGAGTATCAGCGGTGCGCAGATCCGCGAACACCCGTTCCGCTACt       >  1:1323078/1‑150 (MQ=255)
                                                                                                                                            gATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTg      <  1:577457/150‑1 (MQ=255)
                                                                                                                                              tCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTggg    >  2:807215/1‑150 (MQ=255)
                                                                                                                                                cAGCCGGATCTGATCTACCCCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAGCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGaa  <  1:810073/150‑1 (MQ=255)
                                                                                                                                                cAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATTCTCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGATAACCCGTTCCGCTACTGGGaa  <  1:803119/150‑1 (MQ=255)
                                                                                                                                                cAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGaa  >  2:70600/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAA  >  NC_000913/4627681‑4627974

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