Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 194,124 C→A 40.2% I364I (ATC→ATA rseP → inner membrane zinc RIP metalloprotease; RpoE activator, by degrading RseA; cleaved signal peptide endoprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273194,1240CA40.2% 336.8 / 286.8 260I364I (ATC→ATArsePinner membrane zinc RIP metalloprotease; RpoE activator, by degrading RseA; cleaved signal peptide endoprotease
Reads supporting (aligned to +/- strand):  ref base C (88/67);  new base A (54/51);  total (142/118)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 263 total aligned reads displayed.

TATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGATCGAAAAGAT  >  CP009273/193975‑194273
                                                                                                                                                     |                                                                                                                                                     
tataAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATc                                                                                                                                                       >  2:2129783/1‑150 (MQ=255)
 ataAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGACAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCg                                                                                                                                                      >  1:4372717/1‑150 (MQ=255)
   aaaGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGcc                                                                                                                                                    <  1:53760/150‑1 (MQ=255)
   aaaGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGcc                                                                                                                                                    >  2:3240538/1‑150 (MQ=255)
   aaaGTTGTACGCCAGTATGGGCCGTTCAACGCCATAGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATATCCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGcc                                                                                                                                                    <  1:1514142/150‑1 (MQ=255)
        tGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGCAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAggg                                                                                                                                               <  2:1752470/150‑1 (MQ=255)
        tGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAggg                                                                                                                                               >  1:4599171/1‑150 (MQ=255)
        tGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAggg                                                                                                                                               >  1:922690/1‑150 (MQ=255)
        tGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAggg                                                                                                                                               <  1:3101888/150‑1 (MQ=255)
         gTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAgggg                                                                                                                                              >  1:4307370/1‑150 (MQ=255)
         gTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAgggg                                                                                                                                              >  1:3469689/1‑150 (MQ=255)
            cGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCt                                                                                                                                           >  2:849530/1‑150 (MQ=255)
               cAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGTCAAAACGTGGCAGCTGATGAAGCTGACGGTCTGTATGCTGGGAAAATTGATCTCCGGTGATGTGAAACTGAAAAACCTCATTGGGACGATCTCTATAGCCAAGGGGGCTggg                                                                                                                                        <  2:2539870/150‑1 (MQ=255)
                aGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGAGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGa                                                                                                                                       >  2:4120515/1‑150 (MQ=255)
                 gTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGAt                                                                                                                                      >  2:2963454/1‑150 (MQ=255)
                  tATGGGCCGGTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATg                                                                                                                                     >  2:3895981/1‑150 (MQ=255)
                    tGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGAc                                                                                                                                   >  1:4142659/1‑150 (MQ=255)
                       gCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGc                                                                                                                                <  1:4167013/150‑1 (MQ=255)
                          gTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGa                                                                                                                             >  1:2963214/1‑150 (MQ=255)
                            tCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAAc                                                                                                                           >  1:619914/1‑150 (MQ=255)
                            tCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAAc                                                                                                                           <  2:2835004/150‑1 (MQ=255)
                             cAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACt                                                                                                                          >  1:2811112/1‑150 (MQ=255)
                             cAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACt                                                                                                                          >  1:2799238/1‑150 (MQ=255)
                             cAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACt                                                                                                                          >  1:3195473/1‑150 (MQ=255)
                             cAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACt                                                                                                                          <  1:4073197/150‑1 (MQ=255)
                              aaCGCCATCGTCGAAGCCACGGACACAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTc                                                                                                                         <  1:1016648/150‑1 (MQ=255)
                              aaCGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTc                                                                                                                         >  1:3395871/1‑150 (MQ=255)
                              aaCGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTc                                                                                                                         >  1:1082802/1‑150 (MQ=255)
                              aaCGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTc                                                                                                                         >  1:3878709/1‑150 (MQ=255)
                              aaCGCCATAGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTc                                                                                                                         >  1:3042505/1‑150 (MQ=255)
                               aCGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCg                                                                                                                        <  2:767902/150‑1 (MQ=255)
                                cGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCgg                                                                                                                       <  2:1906841/150‑1 (MQ=255)
                                   cATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGgt                                                                                                                    >  1:3812468/1‑150 (MQ=255)
                                    aTCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGgtt                                                                                                                   <  1:3498451/150‑1 (MQ=255)
                                     tcgtcgAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGgttg                                                                                                                  <  1:2051728/150‑1 (MQ=255)
                                      cgtcgAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGgttgt                                                                                                                 >  1:1683288/1‑150 (MQ=255)
                                       gtcgAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGgttgtt                                                                                                                >  1:998921/1‑150 (MQ=255)
                                        tcgAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGttt                                                                                                               >  1:2859253/1‑150 (MQ=255)
                                         cgAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTtta                                                                                                              >  1:3224691/1‑150 (MQ=255)
                                          gAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTttat                                                                                                             <  2:131700/150‑1 (MQ=255)
                                           aaGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTttatt                                                                                                            >  1:849682/1‑150 (MQ=255)
                                            aGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTttatta                                                                                                           >  1:2909469/1‑150 (MQ=255)
                                              ccACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTAcc                                                                                                         >  2:3199890/1‑150 (MQ=255)
                                               cACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCt                                                                                                        >  2:3966458/1‑150 (MQ=255)
                                                aCGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTg                                                                                                       >  1:1711105/1‑150 (MQ=255)
                                                  ggACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGcc                                                                                                     <  1:1384363/150‑1 (MQ=255)
                                                        aaCGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCt                                                                                               >  2:3711806/1‑150 (MQ=255)
                                                         aCGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCtt                                                                                              <  2:1669770/150‑1 (MQ=255)
                                                              gCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCt                                                                                         >  2:4429566/1‑150 (MQ=255)
                                                               cAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCtt                                                                                        <  1:244921/150‑1 (MQ=255)
                                                                aGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGAGGTGAAACTGAACACCCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCGTGCGCtta                                                                                       <  2:1446815/150‑1 (MQ=255)
                                                                  cTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCttatt                                                                                     <  2:2667378/150‑1 (MQ=255)
                                                                   tgatgaAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCttatta                                                                                    >  1:1452869/1‑150 (MQ=255)
                                                                      tgaAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCg                                                                                 >  1:19260/1‑150 (MQ=255)
                                                                        aaGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTg                                                                               >  2:241195/1‑150 (MQ=255)
                                                                          gCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGaa                                                                             >  2:3072545/1‑150 (MQ=255)
                                                                           cTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAAc                                                                            <  2:3028429/150‑1 (MQ=255)
                                                                           cTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAAc                                                                            <  2:573576/150‑1 (MQ=255)
                                                                                    aGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACGGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATa                                                                   <  1:256705/150‑1 (MQ=255)
                                                                                    aGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAATTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATa                                                                   >  2:3854425/1‑150 (MQ=255)
                                                                                    aGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATa                                                                   <  2:3792975/150‑1 (MQ=255)
                                                                                     gTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATaa                                                                  <  2:3007998/150‑1 (MQ=255)
                                                                                     gTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATaa                                                                  >  1:2745693/1‑150 (MQ=255)
                                                                                         gCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTaa                                                              >  1:4417956/1‑150 (MQ=255)
                                                                                          cTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAAc                                                             >  1:715674/1‑150 (MQ=255)
                                                                                          cTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAAc                                                             >  2:4361372/1‑150 (MQ=255)
                                                                                            gtgAATATTGATCACCTGTGATGTGTAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGATGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGATCTTAGGGATAATTAACCt                                                           <  1:89727/148‑1 (MQ=255)
                                                                                               aaaaTTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGtt                                                        <  2:1883098/150‑1 (MQ=255)
                                                                                               aaaaTTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGtt                                                        >  2:654077/1‑150 (MQ=255)
                                                                                                aaaTTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGttt                                                       <  2:762004/150‑1 (MQ=255)
                                                                                                 aaTTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTc                                                      <  1:3672732/150‑1 (MQ=255)
                                                                                                  aTTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTcc                                                     >  1:3684001/1‑150 (MQ=255)
                                                                                                        cACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGcc                                               >  2:2753775/1‑150 (MQ=255)
                                                                                                        cACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGcc                                               >  2:2889408/1‑150 (MQ=255)
                                                                                                          ccGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTAATAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCg                                             >  1:4409168/1‑150 (MQ=255)
                                                                                                                  gtgAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGAc                                     >  1:1570682/1‑150 (MQ=255)
                                                                                                                      caCTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGACGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACgggg                                 <  1:1235214/149‑1 (MQ=255)
                                                                                                                       aCTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACggggg                                >  1:3736592/1‑150 (MQ=255)
                                                                                                                          gAACAACCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCa                             >  2:443957/1‑150 (MQ=255)
                                                                                                                             caacCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATct                          <  1:2129492/150‑1 (MQ=255)
                                                                                                                             caacCTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTAGTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATct                          <  2:1618154/150‑1 (MQ=255)
                                                                                                                                 cTCAGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATctgctg                      >  2:4569157/1‑150 (MQ=255)
                                                                                                                                    aGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTc                   >  2:4095768/1‑150 (MQ=255)
                                                                                                                                    aGTGGGCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTc                   <  2:809704/150‑1 (MQ=255)
                                                                                                                                      tGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCt                 >  2:3667154/1‑150 (MQ=255)
                                                                                                                                        ggCCGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTg               >  1:781834/1‑150 (MQ=255)
                                                                                                                                         gCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGc              <  2:500637/150‑1 (MQ=255)
                                                                                                                                         gCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGc              <  2:500691/150‑1 (MQ=255)
                                                                                                                                          ccGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCg             >  1:2277346/1‑150 (MQ=255)
                                                                                                                                           cGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGa            >  2:1094937/1‑150 (MQ=255)
                                                                                                                                           cGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGa            >  1:56919/1‑150 (MQ=255)
                                                                                                                                           cGATCTCTATAGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGa            >  1:653308/1‑150 (MQ=255)
                                                                                                                                            gATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGAt           >  1:1661821/1‑150 (MQ=255)
                                                                                                                                            gATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGAt           >  2:3680740/1‑150 (MQ=255)
                                                                                                                                            gATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGAt           >  1:4340940/1‑150 (MQ=255)
                                                                                                                                            gATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGAt           >  2:3966049/1‑150 (MQ=255)
                                                                                                                                            gATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGAt           <  2:249588/150‑1 (MQ=255)
                                                                                                                                             aTCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGATc          >  2:2203356/1‑150 (MQ=255)
                                                                                                                                                  tATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGATCGaaaa     <  2:2743602/150‑1 (MQ=255)
                                                                                                                                                  tATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGATCGaaaa     <  2:3889296/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATAAAGTTGTACGCCAGTATGGGCCGTTCAACGCCATCGTCGAAGCCACGGACAAAACGTGGCAGCTGATGAAGCTGACGGTCAGTATGCTGGGAAAATTGATCACCGGTGATGTGAAACTGAACAACCTCAGTGGGCCGATCTCTATCGCCAAGGGGGCTGGGATGACAGCGGAACTCGGGGTTGTTTATTACCTGCCGTTTCTTGCGCTTATTAGCGTGAACTTAGGGATAATTAACCTGTTTCCGTTGCCCGTACTTGACGGGGGGCATCTGCTGTTCCTTGCGATCGAAAAGAT  >  CP009273/193975‑194273

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