Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,922,524 A→T 5.5% V149D (GTT→GAT)  pphA ← serine/threonine‑specific protein phosphatase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,922,5240AT 5.5% 238.8 / 2.9 73V149D (GTT→GAT) pphAserine/threonine‑specific protein phosphatase 1
Reads supporting (aligned to +/- strand):  ref base A (34/35);  new base T (2/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TTGCCAATATCCACGCGATGTCGCAACGGTGTATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAGCGCAATAAACCAGTCGCCG  >  NC_000913/1922377‑1922670
                                                                                                                                                   |                                                                                                                                                  
ttGCCAATATCCACTCTATGTCGAAACGGTGTATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGTCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATaaa                                                                                                                                                    <  5:58619/148‑1 (MQ=255)
  gCCAATATCCACGCGATGTCGCAACGGTGTATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAAc                                                                                                                                                   <  5:194255/147‑1 (MQ=255)
     aaTATCCACGCGATGTCGCAACGGTGTATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGACCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAAGCAACCTCCTTTTGCCATTCATAAAcatc                                                                                                                                                >  3:401030/1‑147 (MQ=255)
       tatCCACGCGATGTCGCAACGGTGTATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACatcatc                                                                                                                                             <  1:314291/148‑1 (MQ=255)
        atCCACGCGATGTCGCAACGGTGTATGACCAAACCAGAAATCATCAGCACCTGTAAATCCCTGCCCTTTTTGGCGTTCACCTAATCCCGCGCGGCTCCACAAGCCCTGATGCACATCAAAGTCATTTTGCCATTCATAAACATCATCt                                                                                                                                            >  5:141190/1‑148 (MQ=255)
          ccACGCGATGTCGCAACGGTGTATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTg                                                                                                                                           >  1:185506/1‑147 (MQ=255)
          ccACGCGATGTCGCAACGGTGTATGACAAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGAGTTCACCTAATCGCGAGAGGCTCCAAAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTg                                                                                                                                           >  4:347635/1‑147 (MQ=255)
             cgcgATGTCGCAACGGTGTATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGAt                                                                                                                                        >  5:178533/1‑147 (MQ=255)
               cgATGTCGCAACGGTGTATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAAt                                                                                                                                     <  2:411123/148‑1 (MQ=255)
                gATGTCGCAACTGTGTATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGATCCACAAGACCTGATGCATATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATc                                                                                                                                    <  3:250167/148‑1 (MQ=255)
                    tCGCAACGGTGTATGTCCTAACCTGATATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGc                                                                                                                                 <  6:397920/147‑1 (MQ=255)
                     cGCAACGGTGTATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCAt                                                                                                                               >  1:223377/1‑148 (MQ=255)
                      gCAACGGTGTATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATg                                                                                                                              <  7:148848/148‑1 (MQ=255)
                        aaCGGTGTATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAg                                                                                                                            <  5:283630/148‑1 (MQ=255)
                          cGGTGTATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCa                                                                                                                          <  3:2719/148‑1 (MQ=255)
                               tATGACCACACCAGAAATGATCAGCAACTGTAATTCCCTGCCCTTTTTGGCGTTCAACTAATCGCGAGCGGCTCCACAAGACCTGATGCTAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAAc                                                                                                                     <  6:141190/148‑1 (MQ=255)
                               tATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAAc                                                                                                                     <  1:308781/148‑1 (MQ=255)
                                      aaaCCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTg                                                                                                              >  4:25155/1‑148 (MQ=255)
                                       aaCCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATATCATCATCTGGATAATCGTCATGAGCAATAACATGCTTGcc                                                                                                            >  4:199741/1‑148 (MQ=255)
                                       aaCCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATATCATCATCTGGATAATCGGCATGAGCAATAACATGCTTGcc                                                                                                            >  4:193499/1‑148 (MQ=255)
                                        aCCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGcc                                                                                                            >  7:32005/1‑148 (MQ=255)
                                        aCCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGcc                                                                                                            <  5:202128/148‑1 (MQ=255)
                                          cAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCgg                                                                                                          <  4:143813/148‑1 (MQ=255)
                                                   tCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACtgt                                                                                                 <  3:25155/148‑1 (MQ=255)
                                                        aCCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATg                                                                                                                  >  7:392420/1‑126 (MQ=255)
                                                        aCCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATg                                                                                                                  <  8:392420/126‑1 (MQ=255)
                                                         ccTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACtt                                                                                           <  8:53040/148‑1 (MQ=255)
                                                           tGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATTAACATCATCTGGATAATCGGCATGAGCAGTAACAGGCTTGCCGGTGCGCCTGTGTACTTc                                                                                          >  8:124303/1‑147 (MQ=255)
                                                           tGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTc                                                                                          >  4:236829/1‑147 (MQ=255)
                                                            gTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAAc                                                                                                                     >  7:411657/1‑119 (MQ=255)
                                                            gTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAAc                                                                                                                     <  8:411657/119‑1 (MQ=255)
                                                            gTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGAAATAACATGCTTGCCGGTGCGACTGTGTACTTc                                                                                          >  1:138562/1‑146 (MQ=255)
                                                             tAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCa                                                                                         <  8:32005/146‑1 (MQ=255)
                                                             tAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCa                                                                                         <  8:60669/146‑1 (MQ=255)
                                                              aaTTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAg                                                                                       >  5:218689/1‑147 (MQ=255)
                                                              aaTTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCATTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAg                                                                                       >  5:106710/1‑147 (MQ=255)
                                                                  cccTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATa                                                                                   <  8:80641/147‑1 (MQ=255)
                                                                   ccTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATaa                                                                                  <  6:218689/147‑1 (MQ=255)
                                                                         cTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAgg                                                                               <  1:32924/144‑1 (MQ=255)
                                                                         cTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAgg                                                                               >  2:32924/1‑144 (MQ=255)
                                                                         cTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAg                                                                                <  4:169801/143‑1 (MQ=255)
                                                                         cTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAg                                                                                >  3:169801/1‑143 (MQ=255)
                                                                           ttttGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAAc                                                                                                                     <  6:202255/104‑1 (MQ=255)
                                                                           ttttGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAAc                                                                                                                     >  5:202255/1‑104 (MQ=255)
                                                                                       ccTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATg                                                                                                                              >  1:74664/1‑83 (MQ=255)
                                                                                       ccTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATg                                                                                                                              <  2:74664/83‑1 (MQ=255)
                                                                                         tAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGc                                                                                                                                 <  6:371380/78‑1 (MQ=255)
                                                                                         tAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGc                                                                                                                                 >  5:371380/1‑78 (MQ=255)
                                                                                          aaTCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCg                                                                                                      <  1:167225/104‑1 (MQ=255)
                                                                                          aaTCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCg                                                                                                      >  2:167225/1‑104 (MQ=255)
                                                                                           aTCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACAttatccagc                                                         <  4:75653/148‑8 (MQ=255)
                                                                                           aTCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTcc                                                           <  2:134520/146‑1 (MQ=255)
                                                                                                       tCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGc                                             >  3:21885/1‑148 (MQ=255)
                                                                                                            aaGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGtt                                        >  5:392569/1‑148 (MQ=255)
                                                                                                            aaGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGACCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGtt                                        >  1:441373/1‑148 (MQ=255)
                                                                                                               aCCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGTCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCt                                     <  3:347635/148‑1 (MQ=255)
                                                                                                                ccTGATGCAAATCAACGTACTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAATGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTg                                    <  3:236829/148‑1 (MQ=255)
                                                                                                                 cTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGt                                   <  2:185506/148‑1 (MQ=255)
                                                                                                                    aTGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACt                                                                                            >  1:181421/1‑88 (MQ=255)
                                                                                                                    aTGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACt                                                                                            <  2:181421/88‑1 (MQ=255)
                                                                                                                    aTGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGa                                >  4:225249/1‑148 (MQ=255)
                                                                                                                         aaTCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATc                           <  4:21885/148‑1 (MQ=255)
                                                                                                                          aTCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCt                          >  8:44861/1‑148 (MQ=255)
                                                                                                                             aaCGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGcc                       >  8:333580/1‑148 (MQ=255)
                                                                                                                                gTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAgc                    >  7:185579/1‑148 (MQ=255)
                                                                                                                                gTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAgc                    >  2:26066/1‑148 (MQ=255)
                                                                                                                                gTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAgc                    >  4:6331/1‑148 (MQ=255)
                                                                                                                                 tCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAgcg                   >  5:222011/1‑148 (MQ=255)
                                                                                                                                        gCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAGCGCAATaa             <  2:77307/147‑1 (MQ=255)
                                                                                                                                         ccATTCATATCATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAGCGCAATAAAcc          <  3:199741/148‑1 (MQ=255)
                                                                                                                                         ccATTCATATCATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAGCGCAATAAAcc          <  3:193499/148‑1 (MQ=255)
                                                                                                                                             tCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAGCGCAATAAAcc          <  3:377921/145‑1 (MQ=255)
                                                                                                                                             tCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAGCGCAATAAAcc          >  4:377921/1‑145 (MQ=255)
                                                                                                                                                tAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAGCGCAATAAACCAGTcgc    >  8:73331/1‑148 (MQ=255)
                                                                                                                                                 aaaCATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAGCGCAATAAACCAGTcgcc   >  8:302739/1‑148 (MQ=255)
                                                                                                                                                 aaaCATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAATTCTGACATTTTTACAGGGGCGGTTTGGCTTGTTTCTGTTGATTATGTGGCAGCGCAATAAACCATTcggcg  >  8:39511/1‑148 (MQ=255)
                                                                                                                                                  aaCATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAGCGCAATAAACCAGTcgccg  >  2:229668/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGCCAATATCCACGCGATGTCGCAACGGTGTATGACCAAACCAGAAATGATCAGCACCTGTAATTCCCTGCCCTTTTTGGCGTTCACCTAATCGCGAGCGGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAGCGCAATAAACCAGTCGCCG  >  NC_000913/1922377‑1922670

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