Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,898,845 G→A 5.6% R140H (CGT→CAT)  yoaD → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase, regulator of cellulose production

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,898,8450GA 5.6% 226.9 / 2.5 71R140H (CGT→CAT) yoaDputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase, regulator of cellulose production
Reads supporting (aligned to +/- strand):  ref base G (36/31);  new base A (2/2);  total (38/33)
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.15e-01

CAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGC  >  NC_000913/1898698‑1898990
                                                                                                                                                   |                                                                                                                                                 
cAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGt                                                                                                                                                  >  4:267744/1‑149 (MQ=255)
cAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGt                                                                                                                                                  <  3:230749/149‑1 (MQ=255)
   aaCTCCAAACTGGGCGATCCATTTGCCTGGTGCAATACGGAACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATCTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTccc                                                                                                                                               <  3:238278/149‑1 (MQ=255)
     cTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCt                                                                                                                                             <  3:352571/149‑1 (MQ=255)
      tCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGCTCGACCGTGCCCt                                                                                                                                             >  6:244794/1‑148 (MQ=255)
      tCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCt                                                                                                                                             >  7:82623/1‑148 (MQ=255)
        cAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGCCATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTCACGATCGACCGTGCCCTGAt                                                                                                                                          >  5:205960/1‑149 (MQ=255)
            cTGTGCGATCCATTGTCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTaaa                                                                                                                                      <  5:98945/149‑1 (MQ=255)
            cTGTGCGATCCATTGGCCTGGGGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGCTTATCGCAATGTGCCCGTCGTGGACATTCTGGCGGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGAGCGACCGTGCCCTGATTaaa                                                                                                                                      <  1:324383/149‑1 (MQ=255)
                  gATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGt                                                                                                                                >  5:162460/1‑149 (MQ=255)
                   aTCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTc                                                                                                                               >  7:173648/1‑149 (MQ=255)
                    tCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTcc                                                                                                                              >  7:180874/1‑149 (MQ=255)
                        ttGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCGGGCTGCACTTCCTGCACCGCAACCCCTTTTACGCCTGACGATCGACCGTGCCCGGATTAAAGGCAGTCCGGtg                                                                                                                          >  5:309741/1‑148 (MQ=255)
                           gCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGtttt                                                                                                                        >  8:62261/1‑148 (MQ=255)
                              tGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCATGCCCTGATTAAAGGCAGTCCGGTTTTGAtt                                                                                                                    >  8:215117/1‑149 (MQ=255)
                                  gCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAAt                                                                                                                >  4:192162/1‑149 (MQ=255)
                                   cAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAAt                                                                                                                <  8:82623/148‑1 (MQ=255)
                                      gACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGAc                                                                                                            >  7:287917/1‑149 (MQ=255)
                                      gACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGAc                                                                                                            >  6:288065/1‑149 (MQ=255)
                                      gACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCATGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGAc                                                                                                            >  5:84934/1‑149 (MQ=255)
                                                  tATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCgg                                                                                                 <  6:162460/148‑1 (MQ=255)
                                                   aTTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCggg                                                                                                >  5:74988/1‑148 (MQ=255)
                                                    ttGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCa                                                                                              <  8:180874/149‑1 (MQ=255)
                                                          cAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGAt                                                                                                                     >  4:116291/1‑120 (MQ=255)
                                                          cAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGAt                                                                                                                     <  3:116291/120‑1 (MQ=255)
                                                           aGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGttt                                                                                                                         <  6:36760/115‑1 (MQ=255)
                                                           aGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGttt                                                                                                                         >  5:36760/1‑115 (MQ=255)
                                                           aGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCaa                                                                                       >  4:337871/1‑149 (MQ=255)
                                                               tttttGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCt                                                                                   >  5:163860/1‑149 (MQ=255)
                                                                   tGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTggg                                                                                <  8:287917/148‑1 (MQ=255)
                                                                   tGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTggg                                                                                <  8:173648/148‑1 (MQ=255)
                                                                   tGGGTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGAGCGTGCCCTGATTAAACGCAGTACGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCTATGCTggg                                                                                <  7:190157/148‑1 (MQ=255)
                                                                         tCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATg                                                                          >  1:325966/1‑148 (MQ=255)
                                                                            cAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCACGACCACTTCTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATg                                                                          >  8:62141/1‑145 (MQ=255)
                                                                            cAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCACGACCACTTCTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATg                                                                          <  7:62141/145‑1 (MQ=255)
                                                                               tgtgCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAgatga                                                                   <  4:105868/149‑1 (MQ=255)
                                                                                 tgCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAgatgat                                                                  <  5:288065/148‑1 (MQ=255)
                                                                                  gCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTa                                                                <  6:302630/149‑1 (MQ=255)
                                                                                   cccGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTa                                                                <  2:325966/148‑1 (MQ=255)
                                                                                    ccGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAAc                                                              <  7:214342/149‑1 (MQ=255)
                                                                                      gtcgtGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACAt                                                            <  2:95243/149‑1 (MQ=255)
                                                                                             aCATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCATGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTa                                                     <  7:215117/149‑1 (MQ=255)
                                                                                               aTTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACt                                                   >  1:389614/1‑149 (MQ=255)
                                                                                                     gCTGAACTTCCTGCACCGCAACGACTTTTACGCCTGCCGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTGGATTCAATGGACGCCTGCAGCGTCCAGGAGCTATGATTGGCTCATGGAGATCATAAACATCGACTTACTGCCgg                                              >  4:278217/1‑148 (MQ=255)
                                                                                                     gCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGAc                                                <  6:163860/146‑1 (MQ=255)
                                                                                                     gCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACgg                                              <  3:337871/148‑1 (MQ=255)
                                                                                                     gCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACgg                                              >  7:34993/1‑148 (MQ=255)
                                                                                                       tGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGAc                                                        <  1:43792/136‑1 (MQ=255)
                                                                                                       tGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGAc                                                        >  2:43792/1‑136 (MQ=255)
                                                                                                         aaCTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTa                                                     <  5:252214/137‑1 (MQ=255)
                                                                                                         aaCTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTa                                                     >  6:252214/1‑137 (MQ=255)
                                                                                                               cTGCACCGCAACCACTTTTACGCAGGAAGATCGACCGCGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCGGGGGCCATGGAGATGATTAACATCGACTTTCTGACGGCAATGCTGCtt                                   <  4:217698/149‑1 (MQ=255)
                                                                                                                tGCACCGCAACCTCTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTGAATGGACGCCTGCAGCggg                                                                                                <  8:336524/87‑1 (MQ=255)
                                                                                                                tGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCggg                                                                                                >  7:336524/1‑87 (MQ=255)
                                                                                                                  cACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGTAATGCTGCTTGAg                                >  1:357069/1‑149 (MQ=255)
                                                                                                                  cACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGa                                 <  8:125710/148‑1 (MQ=255)
                                                                                                                                ttACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGc                   >  2:18654/1‑148 (MQ=255)
                                                                                                                                  aCGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGGTTTGATTCAAAGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGACATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCaaa                >  6:105161/1‑149 (MQ=255)
                                                                                                                                     ccTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCa             <  3:192162/149‑1 (MQ=255)
                                                                                                                                      cTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGCGCCACAACTGCCGCAAATCAg            >  7:75545/1‑149 (MQ=255)
                                                                                                                                      cTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAg            >  5:388783/1‑149 (MQ=255)
                                                                                                                                      cTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAg            >  7:181224/1‑149 (MQ=255)
                                                                                                                                        gACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAg            >  8:50684/1‑147 (MQ=255)
                                                                                                                                          cgttCGACCGTGCCCTGATTAAAGGCAGGCCTGGTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGGCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAGTAgc        <  6:309741/146‑1 (MQ=255)
                                                                                                                                          cGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAGTAgc        <  8:181224/149‑1 (MQ=255)
                                                                                                                                           gATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTATCATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAGTAgcg       >  7:131352/1‑149 (MQ=255)
                                                                                                                                           gATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAGTAgcg       >  7:128575/1‑149 (MQ=255)
                                                                                                                                            aTCGACCATGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCTTGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAGTAgcgc      <  6:84934/149‑1 (MQ=255)
                                                                                                                                               gACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAg   >  8:31332/1‑149 (MQ=255)
                                                                                                                                                aCCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGc  <  3:267744/149‑1 (MQ=255)
                                                                                                                                                aCCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGc  >  4:109928/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCGTCGTGGACATTCTGGCTGAACTTCCTGCACCGCAACCACTTTTACGCCTGACGATCGACCGTGCCCTGATTAAAGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGC  >  NC_000913/1898698‑1898990

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: