Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,627,828 C→A 12.4% 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,8280CA12.4% 298.5 / 14.8 119P172T (CCG→ACG) nadRtrifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (53/49);  new base A (2/14);  total (55/64)
Fisher's exact test for biased strand distribution p-value = 3.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

ATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATA  >  NC_000913/4627679‑4627976
                                                                                                                                                     |                                                                                                                                                    
aTCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATACGCACGGCTGGGATGTGTAGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGc                                                                                                                                                      >  1:322013/1‑150 (MQ=255)
  cGTCTGCGTTGGTAATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGACGGCATGCAGACGAATCCGCACGGCTGGGATGTGTGGAGCAACCCCATCAAAAAGTATATGGCTGAAAAAGGGATCCAGCCg                                                                                                                                                    >  1:14949/1‑150 (MQ=255)
   gTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCgg                                                                                                                                                   <  1:649179/150‑1 (MQ=255)
      ggCGTTGGTTAATGCAAACTTTTACATATCAGAAAAATCTTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTAGCCGCACGGCTGGGATGGGTGGAGCAACGGCATCAAAAAGTGTATGGCTGAAAAAGGGATCCAGCCGgatc                                                                                                                                                <  1:673522/149‑1 (MQ=255)
       gCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatct                                                                                                                                               <  1:1696675/150‑1 (MQ=255)
          ttGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgat                                                                                                                                            <  2:35468/150‑1 (MQ=255)
           tGGTTATTGCAAACTTTTAAATATCAGAAATAAAATCGCATTCATTCTTTCAACGAAGAGGGCATGGAGCCGTATCCGGACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatc                                                                                                                                           <  2:462371/150‑1 (MQ=255)
           tGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCATATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatc                                                                                                                                           <  1:771071/150‑1 (MQ=255)
            ggTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatct                                                                                                                                          <  2:446144/150‑1 (MQ=255)
            ggTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatct                                                                                                                                          <  1:1450158/150‑1 (MQ=255)
              ttATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTac                                                                                                                                        <  1:949934/150‑1 (MQ=255)
                aTTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTacac                                                                                                                                      <  1:1680481/150‑1 (MQ=255)
                   gCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGTCTGAAAAAGGGATTCAGCCGGATCTGATCTACACCTc                                                                                                                                   >  1:1389737/1‑150 (MQ=255)
                    cAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATTTGTGGAGCAAAGGCATCAAAAAGTTTATGGCTGAAAAAGGGAGACAGCCGGATCTGATCTACACATCg                                                                                                                                  >  1:245137/1‑150 (MQ=255)
                    cAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCg                                                                                                                                  <  1:983674/150‑1 (MQ=255)
                    cAAACTTTTAAATATCAGAAAAATATTCGCATACATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGATTATGTCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCg                                                                                                                                  <  1:119705/150‑1 (MQ=255)
                        catttAAATATCAGAAAAATATTCGCATTCATGCATTCAACGAAGAGGGCATGGAGCCTTATCCTCACGGCTCGGATGTGTGGAGCAACGGCATCAAAAAGATTATGGCTGAAAAACGGATCCAGCAGGAACTGATCTACACCTCGgaag                                                                                                                              >  2:146925/3‑150 (MQ=255)
                           ttAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTTTGGAGCAACGGCATCAAAAAGTTTATTTTTTTTAATGTTATGCAGACGTATGTGATCTACACCTCGGaagaag                                                                                                                           <  2:911124/150‑1 (MQ=255)
                             aaaTATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGcc                                                                                                                         <  1:425882/150‑1 (MQ=255)
                               atatCAGAAAAATATTCGCAGTCATGCTATCAACGAAGAGGGAATGGAGCCGTAACAGAACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGATAAAGGGATCCAGACGGATCTGATCTACAACTCGGAAGAATCCGa                                                                                                                       >  1:1001328/1‑150 (MQ=255)
                                 atCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTTTTTTTAATGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATg                                                                                                                     <  2:783942/150‑1 (MQ=255)
                                  tCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATgc                                                                                                                    >  2:1179818/1‑150 (MQ=255)
                                   cAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATgcg                                                                                                                   <  1:1077884/150‑1 (MQ=255)
                                   cAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGCAAGAAGCCGATgcc                                                                                                                   >  2:269682/1‑149 (MQ=255)
                                    aGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATgcgc                                                                                                                  >  1:1539927/1‑150 (MQ=255)
                                    aGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATgcgc                                                                                                                  >  1:1532813/1‑150 (MQ=255)
                                     gAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTTATTAAAATGTTATCCAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGcc                                                                                                                 <  2:507165/150‑1 (MQ=255)
                                       aaaaTATTCGCATTCATGCTGTCAACCAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCcac                                                                                                               <  1:1343219/150‑1 (MQ=255)
                                          atatTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCCGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGt                                                                                                            >  1:472526/1‑150 (MQ=255)
                                            atTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGtat                                                                                                          >  2:1570016/1‑150 (MQ=255)
                                              tCGCATTCATGCTTTCAACGAAGAGGGCAGGGAGCCGTATCCGCACGGCTGGCATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGAACCAGACGGCTCAGATCAACAACACGGAAGAAGCCGATGCGCCACAGTatat                                                                                                        >  2:1452231/1‑150 (MQ=255)
                                               cGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATg                                                                                                       <  1:35120/150‑1 (MQ=255)
                                                 cATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACATTATATGGa                                                                                                     >  1:1412504/1‑150 (MQ=255)
                                                 cATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGa                                                                                                     >  1:1126876/1‑150 (MQ=255)
                                                 cATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGa                                                                                                     <  1:1670510/150‑1 (MQ=255)
                                                 cATTCATGCTATCAACGAAGAGCGCATGGAGACGTATCCGCACGGCTGGGATGTGTTGAGCAACGGCATCAAAAAGTTTATGGCTCAAAAAGGGATCAAGCCGCATCTGATCTACACCTCGGAAGAAGCCGATGCCCCACAGAATATGGa                                                                                                     >  2:1401078/1‑150 (MQ=255)
                                                 cATTCATCCCTTCACCTCCGAGCGCATGGAGCCGCATCCGCACGGCTGGGATGTGTCGAGCAACGGCATCAAAAATTTTTTTTTTTTAAATGTTAGGCATACTTATGTGATCTACACCTCGGAAGAAGCCGTTGCGCAACAGTATATGGa                                                                                                     <  2:1517509/150‑1 (MQ=255)
                                                  aTTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCGGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTTTTTTTAATGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGaa                                                                                                    <  2:1766876/150‑1 (MQ=255)
                                                      aTGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAATGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATc                                                                                                <  1:1231694/150‑1 (MQ=255)
                                                       tGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCt                                                                                               >  1:105148/1‑150 (MQ=255)
                                                         cTTTCAACGAAGATGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTTTTTTTTTTAATGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTgg                                                                                             <  2:1260087/150‑1 (MQ=255)
                                                          tttCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTTTTGTTAATGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTggg                                                                                            <  2:1393226/150‑1 (MQ=255)
                                                          tttCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTggg                                                                                            <  1:1726939/150‑1 (MQ=255)
                                                             cAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCACGGAAGAAGCCGATGCGCCACAGTATATGGAACATCAGGGGaa                                                                                         >  2:1134068/1‑149 (MQ=255)
                                                             cAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGAGGTGTGGAGCAACGGCATCAAAAAGTTTTTGTTTTAAAAAGTGAGGGAGACGTATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAt                                                                                         <  2:1185662/150‑1 (MQ=255)
                                                              aaCGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATc                                                                                        >  1:1051377/1‑150 (MQ=255)
                                                              aaCGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATc                                                                                        >  2:779332/1‑150 (MQ=255)
                                                              aaCGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGAAGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATc                                                                                        >  1:1284864/1‑150 (MQ=255)
                                                                 gAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGTTTAAAAAGTGATCGAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCgag                                                                                     <  2:1525584/150‑1 (MQ=255)
                                                                     agGGCAGGGATCCGTATCCGCACCGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGTTTGTAAAAGTGATCGAGACGTATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACgg                                                                                 <  2:242979/150‑1 (MQ=255)
                                                                      gggCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGt                                                                                >  2:1802688/1‑150 (MQ=255)
                                                                      gggCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGt                                                                                >  2:1576115/1‑150 (MQ=255)
                                                                      gggCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCACCACAGTATATGGAACATCTGGGGATCGAGACGGt                                                                                >  2:1408176/1‑150 (MQ=255)
                                                                       ggCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTg                                                                               >  1:1796042/1‑150 (MQ=255)
                                                                         cATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCt                                                                             <  1:919776/150‑1 (MQ=255)
                                                                         cATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAACCCGATGCGCCACAGAATATGGAACATCTGCGGATCGAGACGGTGCt                                                                             >  2:397384/1‑150 (MQ=255)
                                                                          ttGGAGCCGTAGCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAATGGGATCGAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTg                                                                            <  2:385261/149‑1 (MQ=255)
                                                                             gAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTc                                                                         <  1:1406741/150‑1 (MQ=255)
                                                                             gAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTc                                                                         <  1:344767/150‑1 (MQ=255)
                                                                                  gaaTCCGCACGGCTGGGATGGGGGGAGCAACGGCAGCAAAAAGTTTATGGGTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATcc                                                                    <  1:156775/148‑1 (MQ=255)
                                                                                   tATCCGCACGGCTGGGATGCGTGCCGCAACGGCATCAAAAAGTTTATGGCGGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGGCGATCCg                                                                   <  1:497035/150‑1 (MQ=255)
                                                                                    aTCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTGTTTATAAAGTTATCCAGCCGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGa                                                                  <  2:474791/150‑1 (MQ=255)
                                                                                      ccGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGaaa                                                                >  2:1756032/1‑150 (MQ=255)
                                                                                       cGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAAc                                                               >  2:1574290/1‑150 (MQ=255)
                                                                                       cGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGAGCCGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAAc                                                               <  1:905835/150‑1 (MQ=255)
                                                                                         cACGGCTGGGATGTGTGGAGCAACGGCATCAAAAATTTTTTTTTTTTTAATGTTATGGAGTCGTTTGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGt                                                             <  2:129295/150‑1 (MQ=255)
                                                                                            ggCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTTAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTAcc                                                          <  2:628398/150‑1 (MQ=255)
                                                                                              cTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATTGAGACGGTGCTGGTCGATCCGAAACGTACCtt                                                        <  2:230983/150‑1 (MQ=255)
                                                                                              cTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCCCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCtt                                                        <  1:839095/150‑1 (MQ=255)
                                                                                                gggATGTGTGGAGCAACGGCATCAAAAAGTTTATTTTTTAAAAAGTTATCCAGCCGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTa                                                      <  2:1278810/150‑1 (MQ=255)
                                                                                                  gATGTGTGGAGCAACGGCATCAAAAATTTTATGTCTGAAAAAGGTATCCAGACGTATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATg                                                    <  2:385008/150‑1 (MQ=255)
                                                                                                    tgtgtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAg                                                  <  1:1045079/150‑1 (MQ=255)
                                                                                                      tgtgGAGCAACGGCATCAAAATGTTTTTTTTTTTTAATGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTa                                                <  2:1207781/150‑1 (MQ=255)
                                                                                                      tgtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTa                                                >  1:1330948/1‑150 (MQ=255)
                                                                                                       gtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTAt                                               >  1:1752933/1‑150 (MQ=255)
                                                                                                          gAGCAACGGCATAAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAg                                            <  1:235560/150‑1 (MQ=255)
                                                                                                           aGCAACGTCATCAAAAAGGTTATGGCTCAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAGCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGc                                           <  1:271137/150‑1 (MQ=255)
                                                                                                           aGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAATAAGCCGATGCGCCACAGTATATGGAGCATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTCCCTTTATGAGTATCAGc                                           >  1:303783/1‑150 (MQ=255)
                                                                                                             cAACGGCATCAAAAAGTTTATGGCTGAAATAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCgg                                         <  1:269682/150‑1 (MQ=255)
                                                                                                                  gCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTgcgc                                    <  2:105148/150‑1 (MQ=255)
                                                                                                                   cATCAAAACGTTTATGGCTTAAAAAGTGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAGCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCa                                   <  2:1532813/150‑1 (MQ=255)
                                                                                                                   cATCAAAAAGTTTATTTATTATAATGTTATCGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCa                                   <  2:1196205/150‑1 (MQ=255)
                                                                                                                   cATCAAAAAGTTTATTGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGGCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCa                                   <  1:727065/150‑1 (MQ=255)
                                                                                                                    aTCAAAAAGTTTTTTTTTGTAAAAGGGATCCAGCCGTATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAg                                  <  2:317486/150‑1 (MQ=255)
                                                                                                                      cAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGAt                                >  2:1478736/1‑150 (MQ=255)
                                                                                                                       aaaaaGTTTATGGCTGAAAAAGTGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATc                               <  1:723887/150‑1 (MQ=255)
                                                                                                                       aaaaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGCTCGAGACGGTGCTGGTCGATCCGAAACGTAACTTTATGAGTATCATCGGTGCGCAGATc                               >  1:325548/1‑150 (MQ=255)
                                                                                                                        aaaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATcc                              >  1:1608283/1‑150 (MQ=255)
                                                                                                                           aGTTTATGTTTGAAAAAGTTATCGAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCcgcg                           <  2:488300/150‑1 (MQ=255)
                                                                                                                               tATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGaaaa                       >  2:611020/1‑150 (MQ=255)
                                                                                                                                     tGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGtt                 >  1:1550755/1‑150 (MQ=255)
                                                                                                                                      gAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACCTCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTc                >  1:1566798/1‑150 (MQ=255)
                                                                                                                                      gAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTc                >  1:395001/1‑150 (MQ=255)
                                                                                                                                      gAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTc                <  1:1756032/150‑1 (MQ=255)
                                                                                                                                       aaaaaGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTcc               <  1:432336/150‑1 (MQ=255)
                                                                                                                                       aaaaaGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTcc               >  1:1527282/1‑150 (MQ=255)
                                                                                                                                          aaGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGTAACATCTGGGCATCGAGACGGTGCTGGTCTATACGAAACGTACCTTTATGAGTATCAGCGGTGCGCATCTCCGCGAAAACCGGTGCCGCt            >  2:125640/1‑150 (MQ=255)
                                                                                                                                              gATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGCGACGGGGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTg        >  2:169996/1‑150 (MQ=255)
                                                                                                                                              gATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTg        >  1:574977/1‑150 (MQ=255)
                                                                                                                                                    gCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGCGGATCGAGACGGTGCTGGTCCATCCGAAACGTACCTTAATGAGTATAAGCTGAGCGCAGATCCGCGAAAACCCGATCCGCTACTGGGAata  >  2:1446861/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATA  >  NC_000913/4627679‑4627976

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: