Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,158 G→C 16.2% intergenic (‑141/+30) yedQ ← / ← yodC putative membrane‑anchored diguanylate cyclase/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,028,1580GC16.2% 383.9 / 56.7 141intergenic (‑141/+30)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (48/70);  new base C (11/12);  total (59/82)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 160 total aligned reads displayed.

CCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTG  >  NC_000913/2028059‑2028257
                                                                                                   |                                                                                                   
ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtg                                                                                                     >  2:2150086/1‑100 (MQ=255)
ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtg                                                                                                     <  2:1262693/100‑1 (MQ=255)
ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtg                                                                                                     >  2:547221/1‑100 (MQ=255)
ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtg                                                                                                     <  1:2907847/100‑1 (MQ=255)
ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtg                                                                                                     <  2:3628431/100‑1 (MQ=255)
ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtg                                                                                                     >  2:1866233/1‑100 (MQ=255)
 cTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtgt                                                                                                    <  2:1402893/100‑1 (MQ=255)
  tGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTAATTGCAGGCATTAGCACACCAGGCCGGAtgtg                                                                                                   >  1:3193565/1‑100 (MQ=255)
    gCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                 >  2:454210/1‑100 (MQ=255)
    gCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                 >  1:2474160/1‑100 (MQ=255)
        tAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg                                                                                             <  2:3584157/100‑1 (MQ=255)
         aGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTA‑‑ACATCAGGCCGGATCACATCCGGCCt                                                                                          <  1:2919938/100‑1 (MQ=37)
          gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                           >  2:555874/1‑100 (MQ=255)
          gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                           >  2:3400597/1‑100 (MQ=255)
           cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt                                                                                          <  2:923847/100‑1 (MQ=38)
           cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt                                                                                          <  2:3589106/100‑1 (MQ=39)
           cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                           <  2:2056270/99‑1 (MQ=255)
           cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                          <  2:2768313/100‑1 (MQ=255)
           cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                          <  1:32525/100‑1 (MQ=255)
            gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTg                                                                                         <  2:2745086/100‑1 (MQ=255)
            gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTg                                                                                         <  2:2907171/100‑1 (MQ=255)
            gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTg                                                                                         <  1:3821709/100‑1 (MQ=255)
             tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                        >  2:4294212/1‑100 (MQ=255)
              aCCGCAAAACAGCGATAATTTTAGAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt                                                                                          <  1:632582/97‑1 (MQ=38)
               ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                      <  2:3171410/100‑1 (MQ=255)
               ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                      <  2:70306/100‑1 (MQ=255)
                cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGt                                                                                     <  2:224409/100‑1 (MQ=255)
                 gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGtt                                                                                    >  2:1901152/1‑100 (MQ=255)
                  cAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGACCCGGCCTGATGTTa                                                                                   <  2:1791397/100‑1 (MQ=255)
                   aaaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTaa                                                                                  <  2:1764929/100‑1 (MQ=255)
                   aaaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTaa                                                                                  <  1:1805997/100‑1 (MQ=255)
                   aaaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTaa                                                                                  >  1:569435/1‑100 (MQ=255)
                    aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                 <  2:1117267/100‑1 (MQ=255)
                     aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                <  1:253995/100‑1 (MQ=255)
                     aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                <  2:4298545/100‑1 (MQ=255)
                     aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                <  2:798670/100‑1 (MQ=255)
                     aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                <  2:3018558/100‑1 (MQ=255)
                         gCGATAATTTTATAATGGTTGCATGATAAGACGGGACTATTATTCATTGCAGGCATTAGCACACCAGGCCTGATGTGATCCGGCCTGATGTTAATGAt                                                                              <  1:3695044/98‑1 (MQ=255)
                           gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAccc                                                                          >  2:2130633/1‑100 (MQ=255)
                           gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAccc                                                                          >  1:3826794/1‑100 (MQ=255)
                             tAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCtt                                                                        <  2:1959528/100‑1 (MQ=255)
                                ttttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCTGGCATTAGGACACCAGGCAGGATGAGATCCGGCATGATGTTAATGATACCCTTTca                                                                     >  1:2972817/1‑100 (MQ=255)
                                 tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac                                                                    <  2:237171/100‑1 (MQ=255)
                                  ttATAATGGTGGAATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCAGATGTTAATGATACCCATTaaca                                                                   >  1:891402/1‑100 (MQ=255)
                                   tataATGGTGGCATGATAAGACGCGACTATCATTCACTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCacac                                                                  <  2:3097810/100‑1 (MQ=255)
                                    ataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCGTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACAct                                                                 <  1:2828265/100‑1 (MQ=255)
                                     taATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt                                                                >  1:2914517/1‑100 (MQ=255)
                                      aaTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActtn                                                               <  1:41046/100‑2 (MQ=255)
                                      aaTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActtc                                                               <  2:1142278/100‑1 (MQ=255)
                                          gtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTcc                                                           <  2:1072966/100‑1 (MQ=255)
                                           tggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTAAATGATACCCTTTCACACTTCTTCCg                                                          <  1:572111/100‑1 (MQ=255)
                                            ggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGc                                                         <  2:2019891/100‑1 (MQ=255)
                                             gCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCa                                                        <  1:4241634/100‑1 (MQ=255)
                                                 gATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAAc                                                    <  2:4403499/100‑1 (MQ=255)
                                                 gATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAAc                                                    <  2:2499764/100‑1 (MQ=255)
                                                  aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACg                                                   >  1:2552417/1‑100 (MQ=255)
                                                  aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACg                                                   >  2:1500763/1‑100 (MQ=255)
                                                      gACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActc                                               >  2:2534006/1‑100 (MQ=255)
                                                       aCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                              <  1:2533993/100‑1 (MQ=255)
                                                       aCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                              <  1:1877358/100‑1 (MQ=255)
                                                       aCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                              <  1:3788056/100‑1 (MQ=255)
                                                            aCTAGCATTCATTGCAGGCATTAGCACGCCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccctcc                                         <  1:895007/100‑1 (MQ=255)
                                                                 cattcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCgg                                    <  1:2062428/100‑1 (MQ=255)
                                                                  attcattGCAGGCATTAGCACACCAGGACGGATGTGATGCGGCCTGATGTTAATGATACCATTTCACAATTCTTCCGCACCACGACTCCCCTCCCCCg                                     >  1:4065372/1‑98 (MQ=255)
                                                                    tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAAc                                 <  1:3569225/100‑1 (MQ=255)
                                                                     cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACa                                >  1:2937413/1‑100 (MQ=255)
                                                                     cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACa                                >  1:2328208/1‑100 (MQ=255)
                                                                      attGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACaa                               >  2:582795/1‑100 (MQ=255)
                                                                         gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                            <  1:4317083/100‑1 (MQ=255)
                                                                           aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCg                          <  1:1729995/100‑1 (MQ=255)
                                                                            ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg                         <  1:2229590/100‑1 (MQ=255)
                                                                              cATTAGCACACCATGCCGGATGTGATCCGGCCTGATGTAAATGATACCCGTCCACACTTCTTCCGCAGAACGACTCCCATCCCCCGGAACAAGCTCgg                         <  1:1180236/98‑1 (MQ=255)
                                                                                ttAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACCCTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTc                     <  1:401620/100‑1 (MQ=255)
                                                                                 tAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACTCTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                    <  2:2226848/100‑1 (MQ=255)
                                                                                 tAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                    <  2:4043564/100‑1 (MQ=255)
                                                                                 attaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                    >  1:201597/9‑100 (MQ=25)
                                                                                  ttaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                   >  2:657241/8‑100 (MQ=25)
                                                                                  aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                   >  2:254639/1‑100 (MQ=255)
                                                                                  aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                   >  1:4281889/1‑100 (MQ=255)
                                                                                  aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGAAACCCTTTCACACTTCTTCCGCAGAACGCCTCCCCTCCCCCGGAACAACATCGGTTTCAt                   >  1:3760204/1‑100 (MQ=255)
                                                                                   gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                  >  1:2510603/1‑100 (MQ=255)
                                                                                   gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                  <  1:784684/100‑1 (MQ=255)
                                                                                   gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                  >  2:3775018/1‑100 (MQ=255)
                                                                                     acacCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaa                <  1:1570837/100‑1 (MQ=255)
                                                                                      catcAGGCCGGATCACATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaa                >  2:1191285/4‑99 (MQ=38)
                                                                                       acCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaaa              <  1:3430833/100‑1 (MQ=255)
                                                                                        tcAGGCCGGATCACATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAg             >  1:705522/2‑100 (MQ=255)
                                                                                         cacgCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGc            <  1:2849247/97‑1 (MQ=255)
                                                                                         cAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGc            >  2:2850132/1‑100 (MQ=255)
                                                                                          aGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCt           >  2:2914889/1‑100 (MQ=255)
                                                                                           ggCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACCCTTCTTCCGCAGAACGACTCCCCTCACCCGGAACAAGCTCGGTTTCATGAAAAGCtt          >  1:501782/1‑100 (MQ=255)
                                                                                            gCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTc         <  1:216906/100‑1 (MQ=255)
                                                                                             ccGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTcc        >  1:1196221/1‑100 (MQ=255)
                                                                                             ccGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTcc        <  2:450098/100‑1 (MQ=255)
                                                                                             ccGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTcc        >  2:3369606/1‑100 (MQ=255)
                                                                                                   gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCTCCTCCCCCGGAACACGCTCGGTTCCATGAAAAGCTTCCCCCt    >  1:1840036/1‑98 (MQ=255)
                                                                                                   gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg  <  1:523355/100‑1 (MQ=255)
                                                                                                   gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg  <  2:792040/100‑1 (MQ=255)
                                                                                                   gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg  <  2:1608098/100‑1 (MQ=255)
                                                                                                   gtgATCCGGCCTGATGTTAATGATACCCTTTCACAATTCTTCCGCAGGACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg  <  1:1962928/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTG  >  NC_000913/2028059‑2028257

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