Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,565,055 C→G 5.8% Q179E (CAG→GAG)  yffS → CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,565,0550CG 5.8% 220.9 / 5.3 69Q179E (CAG→GAG) yffSCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (26/39);  new base G (2/2);  total (28/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

AAGTCCCAGCAGAAGCGGAAGAGATGGGAGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAGGGAGAGGAT  >  NC_000913/2564912‑2565202
                                                                                                                                               |                                                                                                                                                    
aaGTCCCAGCAGAAGCGGAAGAGATGGGAGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAggga                                                                                                                                                 >  2:1417156/1‑147 (MQ=255)
    cccAGCAGAAGCGGAAGAGATGGGAGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTGAGGTACTgg                                                                                                                                             >  2:919539/1‑149 (MQ=255)
    cccAGCAGAAGCGGAAGAGATGGGAGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTGAGGTACTgg                                                                                                                                             >  2:526606/1‑149 (MQ=255)
      cagcagAAGCGGAAGAGATGGGAGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAACGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCAAGTTGTGTCAGGTACTGGCt                                                                                                                                           >  1:676305/1‑149 (MQ=255)
          agAAGCGGAAGAGAGGGGAGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCaa                                                                                                                                       <  2:1275151/149‑1 (MQ=255)
                 gAAGAGATGGGGGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGggcgg                                                                                                                                <  1:948030/149‑1 (MQ=255)
                  aaGAGATGGGAGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGggcggc                                                                                                                               <  2:87952/149‑1 (MQ=255)
                  aaGAGATGGGAGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGggcggc                                                                                                                               <  2:331072/149‑1 (MQ=255)
                  aaGAGATGGGAGGTGAAAAAGGGGAAGAAGGGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGggcggc                                                                                                                               <  2:954747/149‑1 (MQ=255)
                         gggAGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATAcc                                                                                                                       >  1:1038986/1‑149 (MQ=255)
                          ggAGGTGAAAAAGTGGAAGAAGGGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGataatag                                                                                                                      <  2:109687/149‑7 (MQ=255)
                            aGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACGCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGAAATGGGAGGCGAATCCCTgg                                                                                                                    >  2:98730/1‑149 (MQ=255)
                            aGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTgg                                                                                                                    >  1:841690/1‑149 (MQ=255)
                                   aaaGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGcaacaa                                                                                                             <  2:145199/149‑1 (MQ=255)
                                   aaaGTGGAAGAAGGGGCAGGAGGGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGcaacaa                                                                                                             <  2:639761/149‑1 (MQ=255)
                                      gTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGc                                                                                                          >  2:1279512/1‑149 (MQ=255)
                                      gTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGGTGTGTCAGGTACTGGCTGCAATGGGCGGGGAATACCTGGGCAACAACGc                                                                                                          >  2:1279590/1‑149 (MQ=255)
                                       tGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCa                                                                                                         <  2:975308/149‑1 (MQ=255)
                                         gaagaaGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAgg                                                                                                       <  2:85838/149‑1 (MQ=255)
                                         gaagaaGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAgg                                                                                                       <  2:423303/149‑1 (MQ=255)
                                                  gcagcaGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGCACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTAcagca                                                                                              <  1:784091/149‑1 (MQ=255)
                                                   cagcagAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTGAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTAcagcag                                                                                             <  1:526606/149‑1 (MQ=255)
                                                   cagcagAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTGAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTAcagcag                                                                                             <  1:919539/149‑1 (MQ=255)
                                                         agCAAGCGCCCAACGATGAAAGAGAGATGGGAGCTACGGAGTGGGTTGACCCGAACCCGAAGGAATACCGCGTGGGCAAGTTGTGTCAGGTACCGGATGCAAAGGCGGCGGGTTAGCTGAGCCGCAACGCAGGGCTAGCACGGAAGATc                                                                                       >  2:116171/1‑149 (MQ=255)
                                                          gCAAGCGCCCAACGATGAAAACGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGAGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCa                                                                                      <  1:1011310/149‑1 (MQ=255)
                                                              gcgcCCAACGTGGAAAAAGCGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCAGGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGGTACAGCAGAAGATCACtgt                                                                                  >  2:1408371/1‑149 (MQ=255)
                                                              gcgcCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACtgt                                                                                  >  1:392966/1‑149 (MQ=255)
                                                              gcgcCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTAGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACtgt                                                                                  >  2:1047157/1‑149 (MQ=255)
                                                                gcCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATAACTgcgt                                                                                >  2:1261506/1‑149 (MQ=255)
                                                                 cccAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGGCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGtt                                                                               >  2:303547/1‑149 (MQ=255)
                                                                 cccAACGATGAAAAAGAGAGGGGAGCTAAGCAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGtt                                                                               <  2:1045838/149‑1 (MQ=255)
                                                                   cAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTgg                                                                             >  2:391606/1‑149 (MQ=255)
                                                                   cAACGACGAAAAAGAGGTGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGCCTACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTgg                                                                             <  1:292084/149‑1 (MQ=255)
                                                                    aaCGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGa                                                                            <  1:311039/149‑1 (MQ=255)
                                                                       gATGGAAAAGAGATGGGGGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATaa                                                                         <  2:281042/149‑1 (MQ=255)
                                                                       gATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATaa                                                                         <  1:826821/149‑1 (MQ=255)
                                                                           aaaaaGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAacgac                                                                     <  1:814721/149‑1 (MQ=255)
                                                                            gaaaGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCAGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGACGATCACTGTGTTGGATAacgacg                                                                    <  2:311940/148‑1 (MQ=255)
                                                                            aaaaGAGATGGGAGGTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAacgacg                                                                    <  2:79978/149‑1 (MQ=255)
                                                                            aaaaGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAacgacg                                                                    <  1:676335/149‑1 (MQ=255)
                                                                            aaaaGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTAAAGCAGAAGTTCACTGTGTTGGATAacgacg                                                                    <  1:761790/149‑1 (MQ=255)
                                                                             aaaGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACATCGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACgg                                                                   <  1:1295267/149‑1 (MQ=255)
                                                                                 agaTGGGAGCTAAGGAGTGGGTTACCCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAAc                                                               <  1:528051/149‑1 (MQ=255)
                                                                                 agaTGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAAc                                                               <  1:448095/149‑1 (MQ=255)
                                                                                 agaTGGGAGCGAAGGGGTGGGTTAACCGGAACCCGAAAAAAGACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAAc                                                               <  1:51587/149‑1 (MQ=255)
                                                                                       gggCTAAGCGGTGGGTTAACCGGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGGCGGCAACCGGAAg                                                         <  2:1013099/147‑1 (MQ=255)
                                                                                            aaGGAGTGGTTTACCCCAAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAAt                                                    <  1:327191/149‑1 (MQ=255)
                                                                                            aaGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAAt                                                    <  2:430289/149‑1 (MQ=255)
                                                                                              ggAGTGGGTTAAGCCGACCCCGACAAAACACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCACGGTTACAGCAGAAGACCACTCTGTTGGATAACGACGGCAACCGGAAGCCAATCa                                                  <  1:419292/149‑1 (MQ=255)
                                                                                               gAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAg                                                 >  1:933201/1‑149 (MQ=255)
                                                                                                 agggtttAACCCGAACCCGAAAAAATACACCGTCGGCATGTTTTGTCAAGTACTTGCTGCAATGAGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAAAAGATCACTTTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTa                                               <  2:290436/144‑1 (MQ=255)
                                                                                                  tGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTaa                                              <  1:391606/149‑1 (MQ=255)
                                                                                                    ggTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACg                                            >  2:1121072/1‑149 (MQ=255)
                                                                                                          cccgaacccgaaAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGt                                      <  2:1300655/149‑1 (MQ=255)
                                                                                                           ccgaacccgaaAAAATACAGCGTCGGCATGTTATGTCAGGTACTGGCTGCAATGGGCGGAGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGtt                                     >  2:1043867/1‑149 (MQ=255)
                                                                                                            ctaacccgaaAAAAGACAGCGTCGGCAGGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGttt                                    <  1:98730/147‑1 (MQ=255)
                                                                                                                cccgaaAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTacc                                <  2:1109144/149‑3 (MQ=255)
                                                                                                                 ccgaaAAAATACAGCGTGGGCATGTTGTGTCAGGGACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCCAACAGTAACGGCGCGTTTTATCg                               >  2:397186/1‑149 (MQ=255)
                                                                                                                  cgaaAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGa                              >  2:1209236/1‑149 (MQ=255)
                                                                                                                      aaaaTACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGa                          >  1:667455/1‑149 (MQ=255)
                                                                                                                            cAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGCCCGATTGAGc                    >  2:28724/1‑149 (MQ=255)
                                                                                                                               cGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTCACGGCGCGTTTTATTGCCTGATAGAGCCgg                 >  1:593516/1‑149 (MQ=255)
                                                                                                                                gTCGGCAAGTTTCGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGc                <  2:513652/149‑1 (MQ=255)
                                                                                                                                     cATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAg           <  1:1417156/149‑1 (MQ=255)
                                                                                                                                       tgttgtGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGACGCCAATCAGTAACGGCGCGGTTTATCGACTGATTGAGCAGGCCAAAggg         >  2:10957/1‑149 (MQ=255)
                                                                                                                                         ttgtGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAGTCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAGGgag       >  2:393692/1‑149 (MQ=255)
                                                                                                                                         ttgtGTCAGGTAATGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGAGAACTATGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGATTTATCGACAGATTGAGCAGGGCCAACGgag       >  1:77644/1‑149 (MQ=255)
                                                                                                                                          tgggtCAGGTACTGGCTGCCATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAGGgaga      <  2:48658/146‑1 (MQ=255)
                                                                                                                                              tCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGGTCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGCTTTATCGACTCATTGAGCAGGCCCAAGGGAGAGgat  >  1:898036/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
AAGTCCCAGCAGAAGCGGAAGAGATGGGAGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGCGAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAGGGAGAGGAT  >  NC_000913/2564912‑2565202

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