Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,073,667 T→C 6.0% V710A (GTC→GCC)  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,073,6670TC 6.0% 200.7 / 2.8 66V710A (GTC→GCC) fluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base T (33/29);  new base C (2/2);  total (35/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 = 9.32e-01

TGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGG  >  NC_000913/2073523‑2073814
                                                                                                                                                |                                                                                                                                                   
tggAAGCCATTAACGGGGCCACCACGGAGGAAGGGGCCGTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTcccc                                                                                                                                                 <  2:339656/149‑1 (MQ=255)
  gAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGGCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCt                                                                                                                                               <  3:97695/149‑1 (MQ=255)
  gAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCt                                                                                                                                               >  3:216627/1‑149 (MQ=255)
     gCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTa                                                                                                                                            >  1:416314/1‑149 (MQ=255)
     gCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTa                                                                                                                                            <  3:177887/149‑1 (MQ=255)
      ccATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTAt                                                                                                                                           <  4:392204/149‑1 (MQ=255)
        aTTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGc                                                                                                                                         >  2:339634/1‑149 (MQ=255)
         ttAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGcc                                                                                                                                        >  8:86644/1‑149 (MQ=255)
           aaCGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTc                                                                                                                                      >  1:392022/1‑149 (MQ=255)
              ggTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCAt                                                                                                                                   >  8:273722/1‑149 (MQ=255)
                  ccaccaCGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCt                                                                                                                                <  6:166742/148‑1 (MQ=255)
                   caccacGGAGGAAGGGGACTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGCGAGCTGTTATCTGCGCAGTGAAAATGCTTATCGTGAAGAAGTCACCCTGTAATCCTCCATGCTGa                                                                                                                              >  2:88598/1‑149 (MQ=255)
                       acGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGAcaca                                                                                                                          <  2:416314/149‑1 (MQ=255)
                        cGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAg                                                                                                                         <  8:190817/149‑1 (MQ=255)
                                ggggCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGa                                                                                                                 <  5:188903/149‑1 (MQ=255)
                                ggggCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGa                                                                                                                 <  7:86644/149‑1 (MQ=255)
                                  ggCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGAc                                                                                                                >  8:330513/1‑148 (MQ=255)
                                    ccTTTGTGCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGAGGAGAGCTGGTATCTGCGCAGTGAACATGATTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGATGACACAGGCAATGGACTAt                                                                                                             <  4:179376/149‑1 (MQ=255)
                                        tGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGAc                                                                                                          <  3:117887/148‑1 (MQ=255)
                                             aGGGGAACAGGATGCAGGCCGTTGCATTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACCCAGGCAATGGACTATGACCGGAtt                                                                                                    <  1:88598/149‑1 (MQ=255)
                                               gggAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATtgt                                                                                                  >  4:345197/1‑149 (MQ=255)
                                                   aCAGGCTGCAGGCCGGTGCCTTTAACTACTCGCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCa                                                                                              <  7:208569/149‑1 (MQ=255)
                                                   aCAGGCTGCAGGCCGGTGCCTTTAACTACTCACTCAACCGGGACATTGATGAGAGCTGGTATCTTCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCa                                                                                              <  1:296842/149‑1 (MQ=255)
                                                       gCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCt                                                                                          >  2:351393/1‑149 (MQ=255)
                                                       gCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCt                                                                                          >  8:325876/1‑149 (MQ=255)
                                                         tGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGCCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGa                                                                                                                 >  5:214587/1‑124 (MQ=255)
                                                         tGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGCCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGa                                                                                                                 <  6:214587/124‑1 (MQ=255)
                                                              gCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAg                                                                                   >  4:389850/1‑149 (MQ=255)
                                                                  gTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCAt                                                                               <  4:216627/149‑1 (MQ=255)
                                                                      cTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGa                                                                           >  2:222427/1‑149 (MQ=255)
                                                                      cTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGCCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGa                                                                           >  3:278491/1‑149 (MQ=255)
                                                                        ttAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGAcc                                                                         >  8:390991/1‑149 (MQ=255)
                                                                          aaCTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCgg                                                                       >  4:148442/1‑149 (MQ=255)
                                                                            ctactCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGgtg                                                                     >  8:76253/1‑149 (MQ=255)
                                                                               ctCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTa                                                                   >  7:417412/1‑148 (MQ=255)
                                                                                  ccTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCGCCTGTATGCCTCCAAGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCCGCCATCAGACCGCTGTAAATg                                                               >  4:376236/1‑149 (MQ=255)
                                                                                    tCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGgt                                                                      >  5:35356/1‑140 (MQ=255)
                                                                                    tCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGgt                                                                      <  6:35356/140‑1 (MQ=255)
                                                                                    tCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGt                                                             >  4:28240/1‑149 (MQ=255)
                                                                                      aaCCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGt                                                             <  8:39179/147‑1 (MQ=255)
                                                                                          gggACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAaac                                                       <  3:345197/149‑1 (MQ=255)
                                                                                            gACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAaacaa                                                     <  2:392022/149‑1 (MQ=255)
                                                                                            gACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAaacaa                                                     <  7:330513/149‑1 (MQ=255)
                                                                                              cAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAacaaca                                                   <  1:222427/149‑1 (MQ=255)
                                                                                                gTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGa                                                                                                                              <  8:335125/72‑1 (MQ=255)
                                                                                                gTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGa                                                                                                                              >  7:335125/1‑72 (MQ=255)
                                                                                                        agCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGtct                                         >  2:390026/1‑149 (MQ=255)
                                                                                                          cTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGCCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCa                                       <  4:278491/149‑1 (MQ=255)
                                                                                                          cTGGTAGCTGCGCAGTGAAAATGCTTCTCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCa                                       <  6:79488/149‑1 (MQ=255)
                                                                                                                tCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTAAGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCAtt                                  >  5:106130/1‑148 (MQ=255)
                                                                                                                    cgcAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCa                                <  7:76253/146‑1 (MQ=255)
                                                                                                                    cgcAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAgg                              <  7:273722/148‑1 (MQ=255)
                                                                                                                    cgcAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAgg                              <  4:187119/148‑1 (MQ=255)
                                                                                                                     gcAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGc                            >  5:56187/1‑149 (MQ=255)
                                                                                                                       aGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCgg                          >  4:155485/1‑149 (MQ=255)
                                                                                                                               tGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTATGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCg                  >  1:75853/1‑149 (MQ=255)
                                                                                                                                  ttATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTTAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGt                >  5:73205/1‑148 (MQ=255)
                                                                                                                                  ttATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGt                >  8:42860/1‑148 (MQ=255)
                                                                                                                                    aTCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCAtctc                   <  5:221020/143‑1 (MQ=255)
                                                                                                                                    aTCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCAtctc                   >  6:221020/1‑143 (MQ=255)
                                                                                                                                    aTCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCa              <  1:351393/148‑1 (MQ=255)
                                                                                                                                     tCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCAc             >  7:219557/1‑148 (MQ=255)
                                                                                                                                         gCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGGAAATGGTGAAACCAACAGCTTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATaa        >  5:398403/1‑149 (MQ=255)
                                                                                                                                         gCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGGGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATaa        >  5:278189/1‑149 (MQ=255)
                                                                                                                                           aGAAGTCCCCCAGTATGCCACCATGCAGACACAGGCAATCGACTATGACCGGATTGTGGCAGGCTCCCGCAGACATCAGACCGGTGTAAATGGTGAAAACAAAAGCGTCCGTCTCAGCATTCAGGGCGGTCAACTCGGTCAAGAGAAc       >  4:137131/1‑148 (MQ=255)
                                                                                                                                               gTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATgg  <  3:148442/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TGGAAGCCATTAACGGTGCCACCACGGAGGAAGGGGCCTTTGTCCAGGGGAACAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGG  >  NC_000913/2073523‑2073814

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