Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,492,267 T→G 5.6% intergenic (+28/+67) cysNC → / ← imm sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase/colicin‑E7 immunity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,492,2670TG 5.6% 229.1 / 4.6 71intergenic (+28/+67)cysNC/immsulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase/colicin‑E7 immunity protein
Reads supporting (aligned to +/- strand):  ref base T (38/29);  new base G (2/2);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

TTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTACCA  >  NC_002947/1492126‑1492412
                                                                                                                                              |                                                                                                                                                 
ttCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCggg                                                                                                                                            >  7:226939/1‑149 (MQ=255)
   ccGTTCGATGTGCCGCTGGATGCTGATCTTGGGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCCGAAGCGGGGtt                                                                                                                                         <  3:283490/149‑1 (MQ=255)
    cGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGttt                                                                                                                                        >  4:309358/1‑149 (MQ=255)
        cGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGtttttt                                                                                                                                     <  6:50449/148‑1 (MQ=255)
         gATGTGCCGCTGTATGCTTATCTTGTGCTCTAGACCCAGGCAAGCAGCGTGGACGATGGCGTGAACCAGGTGCTGGACGTGATTCGTCAGCGTGGCGCGAGCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTc                                                                                                                                    <  5:40877/148‑1 (MQ=255)
           tgtgCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCCGAAGCGGGGttttccttgc                                                                                                                                 <  1:35979/149‑7 (MQ=255)
            gtgCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTAtt                                                                                                                                >  3:4317/1‑149 (MQ=255)
            gtgCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTAtt                                                                                                                                >  7:152860/1‑149 (MQ=255)
             tgCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTAtt                                                                                                                                <  1:288884/148‑1 (MQ=255)
                    ggATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGtt                                                                                                                        >  8:192280/1‑149 (MQ=255)
                        gCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGCCGTGAAACAGCTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGTCGGGGTTTTTTCTATAGCTCAGTTCgga                                                                                                                    >  6:320403/1‑148 (MQ=255)
                         cTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCggg                                                                                                                    <  1:122505/148‑1 (MQ=255)
                          tGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGt                                                                                                                   >  8:162449/1‑148 (MQ=255)
                           gATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTcggc                                                                                                                                               >  1:60571/1‑119 (MQ=255)
                           gATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTcggc                                                                                                                                               <  2:60571/119‑1 (MQ=255)
                           gATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTg                                                                                                                  >  3:57063/1‑148 (MQ=255)
                                   gATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGccc                                                                                                         >  7:234822/1‑149 (MQ=255)
                                     tCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCt                                                                                                        >  7:95192/1‑148 (MQ=255)
                                     tCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCt                                                                                                        <  8:95192/148‑1 (MQ=255)
                                       gACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCggg                                                                                                                    <  7:54024/134‑1 (MQ=255)
                                       gACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCggg                                                                                                                    >  8:54024/1‑134 (MQ=255)
                                        acacCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCaga                                                                                                    <  7:122101/149‑1 (MQ=255)
                                        acacCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCaga                                                                                                    <  7:281378/149‑1 (MQ=255)
                                          acCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGt                                                                                                  >  6:230468/1‑149 (MQ=255)
                                            ccAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGtt                                                                                                 >  2:62557/1‑148 (MQ=255)
                                            ccAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGtt                                                                                                 >  8:68002/1‑148 (MQ=255)
                                            ccAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGtt                                                                                                 >  5:18191/1‑148 (MQ=255)
                                                   accaGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCTGGGCCt                                                                                         <  8:152860/149‑1 (MQ=255)
                                                   accaGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCt                                                                                         >  5:100939/1‑149 (MQ=255)
                                                     caGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTg                                                                                       <  2:120278/149‑1 (MQ=255)
                                                            gACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCgtgt                                                                                <  4:57063/149‑1 (MQ=255)
                                                               gAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCTATGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTATATTGATCAGTTCGGGTCAAAAGGCCATCAGAGTTCAGGGCCTTGATCgtgt                                                                                >  1:90365/1‑147 (MQ=255)
                                                               gAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTCTTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAg                                                                             <  1:83020/149‑1 (MQ=255)
                                                                 aGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCt                                                                           >  1:8471/1‑149 (MQ=255)
                                                                       gAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCtt                                                                     >  6:230114/1‑149 (MQ=255)
                                                                          actggggCTGGTCGTGCTGCGTGAGCGTGTCGGGTTCTTAGCGGTAGCCTGCTAGGTAAACCCAGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGa                                                                   <  4:4317/142‑1 (MQ=255)
                                                                          aCAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGa                                                                   >  8:132032/1‑148 (MQ=255)
                                                                          aCAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGa                                                                   >  8:136030/1‑148 (MQ=255)
                                                                             ggTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTATCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAAcc                                                               <  7:71710/149‑1 (MQ=255)
                                                                              gTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCg                                                              >  2:52431/1‑149 (MQ=255)
                                                                                 cTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTAt                                                                                                                                 >  2:306674/1‑79 (MQ=255)
                                                                                 cTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTAt                                                                                                                                 <  1:306674/79‑1 (MQ=255)
                                                                                 cTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTc                                                            <  5:230114/148‑1 (MQ=255)
                                                                                    gACGTGCTGCTTCAGCGTTGCGCTATCTAATCGATTGCCTGCAAGGAAAACCGCGCTTCGGCTGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCtt                                                          <  5:320403/147‑1 (MQ=255)
                                                                                    gACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGtt                                                                                                 >  1:219327/1‑108 (MQ=255)
                                                                                    gACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGtt                                                                                                 <  2:219327/108‑1 (MQ=255)
                                                                                    gACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTa                                                        >  4:141929/1‑149 (MQ=255)
                                                                                    gACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTa                                                        <  7:192280/149‑1 (MQ=255)
                                                                                      cGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTAcc                                                      <  2:90365/149‑1 (MQ=255)
                                                                                       gTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCg                                                     >  2:215875/1‑149 (MQ=255)
                                                                                         gctgcGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGtt                                                   >  4:39096/1‑149 (MQ=255)
                                                                                          ctgcGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAggg                                                                                            <  2:191781/107‑1 (MQ=255)
                                                                                           tgcGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGACTTTACAGttcg                                                 >  6:258415/1‑147 (MQ=255)
                                                                                                 aGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAgc                                            <  5:258415/148‑1 (MQ=255)
                                                                                                       gcgcGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCAtt                                     <  7:132032/149‑1 (MQ=255)
                                                                                                       gcgcGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCAtt                                     <  7:136030/149‑1 (MQ=255)
                                                                                                       gcgcGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGACTAGCCGCCAtt                                     >  1:57832/1‑149 (MQ=255)
                                                                                                        cgcgATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGGTCAGGGCCTTGATCGTGTTAGCTGGGCTTGAAACCGTCATTACAGTTGGCTAGCCGCCATTg                                    >  6:14294/1‑149 (MQ=255)
                                                                                                                 aGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGa                               <  1:52431/145‑1 (MQ=255)
                                                                                                                 aGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTAt                           >  6:14808/1‑149 (MQ=255)
                                                                                                                    gcTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTc                        <  1:215875/149‑1 (MQ=255)
                                                                                                                        gCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTa                      <  3:141929/147‑1 (MQ=255)
                                                                                                                         ccTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTa                   <  5:230468/149‑1 (MQ=255)
                                                                                                                               aaGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACtttt              >  4:117615/1‑149 (MQ=255)
                                                                                                                                 aaaaaaCCCCGCCGAAGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACttttt             >  2:57414/3‑148 (MQ=255)
                                                                                                                                  gAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACttttt             >  8:50015/1‑147 (MQ=255)
                                                                                                                                   aaaaCCCCGCCGAAGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAg          >  6:238514/1‑149 (MQ=255)
                                                                                                                                        cccGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTa     >  7:15966/1‑149 (MQ=255)
                                                                                                                                        cccGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTa     >  7:11308/1‑149 (MQ=255)
                                                                                                                                        cccGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTa     >  5:79436/1‑149 (MQ=255)
                                                                                                                                           gCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTACCa  >  7:202002/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
TTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTACCA  >  NC_002947/1492126‑1492412

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