Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,072,813 G→T 5.9% L425L (CTG→CTT flu → CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,072,8130GT 5.9% 237.4 / 2.7 68L425L (CTG→CTTfluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base T (2/2);  total (34/34)
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.84e-01

CACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTG  >  NC_000913/2072666‑2072954
                                                                                                                                                   |                                                                                                                                             
caccacCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTg                                                                                                                                               <  5:39273/148‑1 (MQ=255)
 accaccGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGa                                                                                                                                              <  6:53193/148‑1 (MQ=255)
   caccGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAAc                                                                                                                                            >  6:54659/1‑148 (MQ=255)
    accGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACg                                                                                                                                           >  5:322700/1‑148 (MQ=255)
        tATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCgg                                                                                                                                       >  3:357893/1‑148 (MQ=255)
          tCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTg                                                                                                                                     <  3:302338/148‑1 (MQ=255)
           ccATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTgc                                                                                                                                    >  1:383420/1‑147 (MQ=255)
            cATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTg                                                                                                                                     >  4:69075/1‑146 (MQ=255)
            cATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTg                                                                                                                                     <  3:69075/146‑1 (MQ=255)
               gggAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGt                                                                                                                                      <  8:379547/142‑1 (MQ=255)
               gggAAATGGCGGTGTACTGCGGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGt                                                                                                                                      >  7:379547/1‑142 (MQ=255)
                  aaaTGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGcc                                                                                                                             >  6:103963/1‑148 (MQ=255)
                    aTGGCTGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCAc                                                                                                                           >  4:135567/1‑148 (MQ=255)
                       gCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGCCGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGa                                                                                                                        >  1:29961/1‑148 (MQ=255)
                             tACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCAc                                                                                                                  <  2:56794/148‑1 (MQ=255)
                                    ggCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTg                                                                                                                                               >  5:50741/1‑112 (MQ=255)
                                    ggCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTg                                                                                                                                               <  6:50741/112‑1 (MQ=255)
                                      ccGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATgg                                                                                                         <  5:66023/148‑1 (MQ=255)
                                         aTTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCgg                                                                                                      <  4:357893/148‑1 (MQ=255)
                                           tCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGAc                                                                                                    >  1:251047/1‑148 (MQ=255)
                                                        tCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGcc                                                                                         <  2:29961/146‑1 (MQ=255)
                                                        tCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGcc                                                                                         <  6:322700/146‑1 (MQ=255)
                                                                     gAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTTAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTg                                                                          >  8:262987/1‑148 (MQ=255)
                                                                     gAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTTAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTg                                                                          >  8:269714/1‑148 (MQ=255)
                                                                      aGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTgg                                                                         <  7:21969/148‑1 (MQ=255)
                                                                      aGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTgg                                                                         <  2:251047/148‑1 (MQ=255)
                                                                          aCGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCggg                                                                     >  6:214858/1‑148 (MQ=255)
                                                                           cGGAAAGGCATTCAGTATCGGAGGCGGTCAGTCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGc                                                                    <  5:54659/148‑1 (MQ=255)
                                                                           cGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCg                                                                       <  7:131566/145‑1 (MQ=255)
                                                                           cGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCg                                                                       >  8:131566/1‑145 (MQ=255)
                                                                            ggAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACgg                                                                                                                         >  7:114305/1‑94 (MQ=255)
                                                                            ggAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACgg                                                                                                                         <  8:114305/94‑1 (MQ=255)
                                                                            ggAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGca                                                                   <  3:135567/148‑1 (MQ=255)
                                                                             gAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTaa                                                                                                             <  8:158542/105‑1 (MQ=255)
                                                                             gAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTaa                                                                                                             >  7:158542/1‑105 (MQ=255)
                                                                                 ggCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGcaccacc                                                              <  5:103963/148‑1 (MQ=255)
                                                                                  gCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGcaccacca                                                             >  5:299453/1‑148 (MQ=255)
                                                                                     ttCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGc                                                          <  5:15460/148‑1 (MQ=255)
                                                                                      tCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGcacc                                                                 <  1:166262/140‑1 (MQ=255)
                                                                                      tCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGcacc                                                                 >  2:166262/1‑140 (MQ=255)
                                                                                         gTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGaa                                                      >  4:409190/1‑148 (MQ=255)
                                                                                          tATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAAt                                                     <  3:409190/148‑1 (MQ=255)
                                                                                            tCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATa                                                    <  5:362880/147‑1 (MQ=255)
                                                                                             cGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGt                                                                                                 <  2:406148/101‑1 (MQ=255)
                                                                                             cGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGt                                                                                                 >  1:406148/1‑101 (MQ=255)
                                                                                                  gCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCAc                                                                                                                           >  2:369281/1‑70 (MQ=255)
                                                                                                  gCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCAc                                                                                                                           <  1:369281/70‑1 (MQ=255)
                                                                                                       cAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATAc                                        >  5:283618/1‑148 (MQ=255)
                                                                                                                     gATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCgg                           >  1:166444/1‑147 (MQ=255)
                                                                                                                     gATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCgg                           >  6:28510/1‑147 (MQ=255)
                                                                                                                     gATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCgg                           >  6:164162/1‑147 (MQ=255)
                                                                                                                     gATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCgg                           >  5:236859/1‑147 (MQ=255)
                                                                                                                      aTGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCgg                           >  6:93337/1‑146 (MQ=255)
                                                                                                                      aTGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCgg                           <  5:93337/146‑1 (MQ=255)
                                                                                                                           ggAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACg                    >  7:183300/1‑148 (MQ=255)
                                                                                                                            gAAAAAGGCAGTTCATTCACGCTTAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACgg                   <  7:262987/148‑1 (MQ=255)
                                                                                                                            gAAAAAGGCAGTTCATTCACGCTTAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACgg                   <  7:269714/148‑1 (MQ=255)
                                                                                                                            gAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACgg                   <  6:236859/148‑1 (MQ=255)
                                                                                                                                aaGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTAcc                                   >  2:124574/1‑128 (MQ=255)
                                                                                                                                aaGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTAcc                                   <  1:124574/128‑1 (MQ=255)
                                                                                                                                aaGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGaa               <  5:214858/148‑1 (MQ=255)
                                                                                                                                    cAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGaacaac           >  3:145047/1‑148 (MQ=255)
                                                                                                                                     aGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACAGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGaacaac           >  5:62444/1‑147 (MQ=255)
                                                                                                                                       ttcattcaCGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGAt        <  2:166444/148‑1 (MQ=255)
                                                                                                                                       ttcattcaCGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGAt        <  5:164162/148‑1 (MQ=255)
                                                                                                                                       ttcattcaCGCTGAACGCCGGTGACACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGa         >  7:408153/1‑147 (MQ=255)
                                                                                                                                        tcattcaCGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGaac              <  5:215463/141‑1 (MQ=255)
                                                                                                                                        tcattcaCGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGaac              >  6:215463/1‑141 (MQ=255)
                                                                                                                                            tcaCGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTAc                                    <  7:64023/115‑1 (MQ=255)
                                                                                                                                            tcaCGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTAc                                    >  8:64023/1‑115 (MQ=255)
                                                                                                                                              aCGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTg  <  8:183300/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTG  >  NC_000913/2072666‑2072954

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