Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 733,154 T→C 5.6% L1191P (CTG→CCG)  rhsC → Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913733,1540TC 5.6% 213.6 / 2.4 71L1191P (CTG→CCG) rhsCRhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (36/31);  new base C (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.59e-01

AAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATC  >  NC_000913/733006‑733299
                                                                                                                                                    |                                                                                                                                                 
aaaaTTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCt                                                                                                                                                   <  4:169286/149‑1 (MQ=255)
       aCCTGTACCACTGCGACCATGGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGCGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGACCCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCgggg                                                                                                                                            <  7:401063/149‑1 (MQ=255)
       aCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCgggg                                                                                                                                            <  5:267414/149‑1 (MQ=255)
       aCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCggg                                                                                                                                             <  1:206858/148‑1 (MQ=255)
        ccTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGgc                                                                                                                                           <  7:265139/149‑1 (MQ=255)
            tACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCgggg                                                                                                                                            >  6:170476/1‑144 (MQ=255)
            tACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGgcagca                                                                                                                                       <  5:77925/149‑1 (MQ=255)
            tACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCAGCCTGCCgggg                                                                                                                                            <  5:170476/144‑1 (MQ=255)
               cACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTa                                                                                                                                    <  1:60201/149‑1 (MQ=255)
                aCTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTat                                                                                                                                   >  1:413769/1‑149 (MQ=255)
                    cGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTatgat                                                                                                                                >  4:182155/1‑148 (MQ=255)
                    cGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTatgat                                                                                                                                <  7:388088/148‑1 (MQ=255)
                         aTCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATgaggag                                                                                                                          <  1:203745/149‑1 (MQ=255)
                               gCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGc                                                                                                                    >  2:222040/1‑149 (MQ=255)
                               gCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGc                                                                                                                    <  2:373942/149‑1 (MQ=255)
                                    ccGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTa                                                                                                               <  2:413769/149‑1 (MQ=255)
                                       cTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGTGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTa                                                                                                            <  7:387343/149‑1 (MQ=255)
                                            ccTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAAcc                                                                                                       <  1:222040/149‑1 (MQ=255)
                                            ccTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACATGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAAcc                                                                                                       >  8:146298/1‑149 (MQ=255)
                                              ttATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGc                                                                                                     >  6:112915/1‑149 (MQ=255)
                                                    gCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGCTGAGTCCGGCCTGTATTACAACCGCCACCGc                                                                                               >  1:341736/1‑149 (MQ=255)
                                                    gCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGt                                                                                                                <  4:222381/132‑1 (MQ=255)
                                                    gCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGt                                                                                                                >  3:222381/1‑132 (MQ=255)
                                                    gCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCCGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGc                                                                                               >  6:261307/1‑149 (MQ=255)
                                                    gCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCCGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGc                                                                                               >  6:261305/1‑149 (MQ=255)
                                                        ggAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCtat                                                                                            >  5:353626/1‑148 (MQ=39)
                                                        ggAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCtat                                                                                            >  8:35812/1‑148 (MQ=39)
                                                         gAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCtatta                                                                                          <  7:35812/149‑1 (MQ=39)
                                                          aaGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTAc                                                                                                           >  6:266202/1‑131 (MQ=40)
                                                          aaGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTAc                                                                                                           <  5:266202/131‑1 (MQ=40)
                                                              gACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGAcc                                                                                     >  4:98435/1‑149 (MQ=39)
                                                                    aGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGgc                                                                                                                                           <  6:271074/89‑1 (MQ=255)
                                                                    aGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGgc                                                                                                                                           >  5:271074/1‑89 (MQ=255)
                                                                    aGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGc                                                                               <  3:352767/149‑1 (MQ=38)
                                                                    aGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGc                                                                               <  5:112915/149‑1 (MQ=38)
                                                                          gTACGCAGCATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTCCAACCGCCACCGCGATTATGCCCCGCTGCAGGGGc                                                                         >  1:202479/1‑149 (MQ=255)
                                                                          gTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGc                                                                         >  1:175319/1‑149 (MQ=34)
                                                                             cGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGat                                                                      <  6:353626/149‑1 (MQ=33)
                                                                                aGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatata                                                                   >  7:314430/1‑149 (MQ=34)
                                                                                   aTACGATGAATGGGGTAACCTGCTGAATGAAGAGACCCCGCATCAGCTGCAGCACCTTATCCGCCTGCCGCGGCAGCAGTATCATGAGGAGTCCGGCCTGAATTACAACCGCCACCGCTATTATGACCCGCTGCAGGCGCGATATATCa                                                                >  5:279515/1‑149 (MQ=11)
                                                                                          gAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATACGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGa                                                         >  5:370047/1‑149 (MQ=21)
                                                                                            aTGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATc                                                       <  2:175319/149‑1 (MQ=34)
                                                                                                ggTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAt                                                   >  4:371856/1‑149 (MQ=34)
                                                                                                 gTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCCGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAtt                                                  <  5:261305/149‑1 (MQ=18)
                                                                                                 gTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCCGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAtt                                                  <  5:261307/149‑1 (MQ=18)
                                                                                                   aaCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTgg                                                >  3:367962/1‑149 (MQ=32)
                                                                                                    aCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTgg                                                <  7:27623/148‑1 (MQ=32)
                                                                                                    aCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTggg                                               <  8:289513/149‑1 (MQ=17)
                                                                                                       tgctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCt                                             >  2:378654/1‑148 (MQ=32)
                                                                                                        gctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTg                                            >  5:85824/1‑148 (MQ=32)
                                                                                                             atgaaTAGAACCCGAATAATCTGCATGAGCTTATCCGCCTGCCGTGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCATGATCCGATTGGGCTGAAAggt                                      <  2:341736/149‑2 (MQ=11)
                                                                                                                aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGa                                   >  4:232775/1‑149 (MQ=34)
                                                                                                                 agagCACCCGCATCAGCTGCAGCAGCTAATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGAt                                  >  5:348407/1‑149 (MQ=21)
                                                                                                                 agagAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGAt                                  >  5:97902/1‑149 (MQ=34)
                                                                                                                   agaACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATg                                 >  7:22105/1‑148 (MQ=34)
                                                                                                                     aaCCCGCATCAGCTTCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCTCTCAGGATCCGATTGGGCTGAAAGGGGGATGGaa                              <  7:146298/149‑1 (MQ=21)
                                                                                                                              cAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTAt                       >  7:131383/1‑147 (MQ=34)
                                                                                                                               aGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCa                     <  8:131383/148‑1 (MQ=34)
                                                                                                                                gCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTAt                       <  2:346668/145‑1 (MQ=34)
                                                                                                                                gCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTAt                       >  1:346668/1‑145 (MQ=34)
                                                                                                                                gCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAg                    >  4:206477/1‑148 (MQ=34)
                                                                                                                                  tGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTAt                 >  1:65469/1‑149 (MQ=34)
                                                                                                                                  tGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATAACTCAGGATCCGATTGGGCTGAAAGGGGGATGGATTTTTTATCAGTAt                 >  2:50558/1‑149 (MQ=21)
                                                                                                                                     agcagcTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACACGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCg              >  7:146670/1‑149 (MQ=21)
                                                                                                                                      gcagcTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGa                                                                       >  3:112480/1‑91 (MQ=31)
                                                                                                                                      gcagcTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGa                                                                       <  4:112480/91‑1 (MQ=31)
                                                                                                                                       cagcTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGtt            <  1:378654/149‑1 (MQ=34)
                                                                                                                                       cagcTAATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGtt            <  6:348407/149‑1 (MQ=21)
                                                                                                                                             aTCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATcc      >  2:177148/1‑149 (MQ=34)
                                                                                                                                             aTCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATcc      >  8:242594/1‑149 (MQ=34)
                                                                                                                                                cGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTCTTACCACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGAt   <  2:202479/149‑1 (MQ=255)
                                                                                                                                                 gCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGAtc  >  1:166688/1‑149 (MQ=38)
                                                                                                                                                    |                                                                                                                                                 
AAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATC  >  NC_000913/733006‑733299

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: