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

AAGTCCCAGCAGAAGCGGAAGAGATGGGAGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAGGGAGAGGAT  >  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)
                         gggAGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATAcc                                                                                                                       >  1:1038986/1‑149 (MQ=255)
                          ggAGGTGAAAAAGTGGAAGAAGGGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑Gataatag                                                                                                                      <  2:109687/149‑7 (MQ=255)
                            aGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACGCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGAAATGGGAGGC‑GAATCCCTgg                                                                                                                    >  2:98730/1‑149 (MQ=255)
                            aGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTgg                                                                                                                    >  1:841690/1‑149 (MQ=255)
                                   aaaGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGcaacaa                                                                                                             <  2:145199/149‑1 (MQ=255)
                                   aaaGTGGAAGAAGGGGCAGGAGGGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGcaacaa                                                                                                             <  2:639761/149‑1 (MQ=255)
                                      gTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGc                                                                                                          >  2:1279512/1‑149 (MQ=255)
                                      gTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGGTGTGTCAGGTACTGGCTGCAATGGGCGGG‑GAATACCTGGGCAACAACGc                                                                                                          >  2:1279590/1‑149 (MQ=255)
                                       tGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCa                                                                                                         <  2:975308/149‑1 (MQ=255)
                                         gaagaaGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAgg                                                                                                       <  2:85838/149‑1 (MQ=255)
                                         gaagaaGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAgg                                                                                                       <  2:423303/149‑1 (MQ=255)
                                                  gcagcaGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGCACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTAcagca                                                                                              <  1:784091/149‑1 (MQ=255)
                                                   cagcagAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTGAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTAcagcag                                                                                             <  1:526606/149‑1 (MQ=255)
                                                   cagcagAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTGAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTAcagcag                                                                                             <  1:919539/149‑1 (MQ=255)
                                                         agCAAGCGCCCAACGATGAAAGAGAGATGGGAGCTACGGAGTGGGTTGACCCGAACCCGAAGGAATACCGCGTGGGCAAGTTGTGTCAGGTACCGGATGCAA‑AGGCGGCGGGTTAGCTGAGCCGCAACGCAGGGCTAGCACGGAAGATc                                                                                       >  2:116171/1‑149 (MQ=255)
                                                          gCAAGCGCCCAACGATGAAAACGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGAGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCa                                                                                      <  1:1011310/149‑1 (MQ=255)
                                                              gcgcCCAACGTGGAAAAAGCGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCAGGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGGTACAGCAGAAGATCACtgt                                                                                  >  2:1408371/1‑149 (MQ=255)
                                                              gcgcCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACtgt                                                                                  >  1:392966/1‑149 (MQ=255)
                                                              gcgcCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTAGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACtgt                                                                                  >  2:1047157/1‑149 (MQ=255)
                                                                gcCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATAACTgcgt                                                                                >  2:1261506/1‑149 (MQ=255)
                                                                 cccAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGGCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGtt                                                                               >  2:303547/1‑149 (MQ=255)
                                                                 cccAACGATGAAAAAGAGAGGGGAGCTAAGCAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGtt                                                                               <  2:1045838/149‑1 (MQ=255)
                                                                   cAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTgg                                                                             >  2:391606/1‑149 (MQ=255)
                                                                   cAACGACGAAAAAGAGGTGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GCCTACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTgg                                                                             <  1:292084/149‑1 (MQ=255)
                                                                    aaCGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGa                                                                            <  1:311039/149‑1 (MQ=255)
                                                                       gATGGAAAAGAGATGGGGGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATaa                                                                         <  2:281042/149‑1 (MQ=255)
                                                                       gATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATaa                                                                         <  1:826821/149‑1 (MQ=255)
                                                                           aaaaaGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAacgac                                                                     <  1:814721/149‑1 (MQ=255)
                                                                            gaaaGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCAGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGACGATCACTGTGTTGGATAacgacg                                                                    <  2:311940/148‑1 (MQ=255)
                                                                            aaaaGAGATGGGAGGTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAacgacg                                                                    <  2:79978/149‑1 (MQ=255)
                                                                            aaaaGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAacgacg                                                                    <  1:676335/149‑1 (MQ=255)
                                                                            aaaaGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTAAAGCAGAAGTTCACTGTGTTGGATAacgacg                                                                    <  1:761790/149‑1 (MQ=255)
                                                                             aaaGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACATCGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACgg                                                                   <  1:1295267/149‑1 (MQ=255)
                                                                                 agaTGGGAGCTAAGGAGTGGGTTACCCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAAc                                                               <  1:528051/149‑1 (MQ=255)
                                                                                 agaTGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAAc                                                               <  1:448095/149‑1 (MQ=255)
                                                                                 agaTGGGAGCGAAGGGGTGGGTTAACCGGAACCCGAAAAAAGACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAAc                                                               <  1:51587/149‑1 (MQ=255)
                                                                                       gggCTAAGCGGTGGGTTAACCGGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGGCGGCAACCGGAAg                                                         <  2:1013099/147‑1 (MQ=255)
                                                                                            aaGGAGTGGTTTACCCCAAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAAt                                                    <  1:327191/149‑1 (MQ=255)
                                                                                            aaGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAAt                                                    <  2:430289/149‑1 (MQ=255)
                                                                                              ggAGTGGGTTAAGCCGACCCCGACAAAACACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCACGGTTACAGCAGAAGACCACTCTGTTGGATAACGACGGCAACCGGAAGCCAATCa                                                  <  1:419292/149‑1 (MQ=255)
                                                                                               gAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAg                                                 >  1:933201/1‑149 (MQ=255)
                                                                                                 agggtttAACCCGAACCCGAAAAAATACACCGTCGGCATGTTTTGTCAAGTACTTGCTGCAATGAGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAAAAGATCACTTTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTa                                               <  2:290436/144‑1 (MQ=255)
                                                                                                  tGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTaa                                              <  1:391606/149‑1 (MQ=255)
                                                                                                    ggTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACg                                            >  2:1121072/1‑149 (MQ=255)
                                                                                                          cccgaacccgaaAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGt                                      <  2:1300655/149‑1 (MQ=255)
                                                                                                           ccgaacccgaaAAAATACAGCGTCGGCATGTTATGTCAGGTACTGGCTGCAATGGGCGGA‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGtt                                     >  2:1043867/1‑149 (MQ=255)
                                                                                                            ctaacccgaaAAAAGACAGCGTCGGCAGGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGttt                                    <  1:98730/147‑1 (MQ=255)
                                                                                                                cccgaaAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTacc                                <  2:1109144/149‑3 (MQ=255)
                                                                                                                 ccgaaAAAATACAGCGTGGGCATGTTGTGTCAGGGACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCCAACAGTAACGGCGCGTTTTATCg                               >  2:397186/1‑149 (MQ=255)
                                                                                                                  cgaaAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGa                              >  2:1209236/1‑149 (MQ=255)
                                                                                                                      aaaaTACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGa                          >  1:667455/1‑149 (MQ=255)
                                                                                                                            cAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGCCCGATTGAGc                    >  2:28724/1‑149 (MQ=255)
                                                                                                                               cGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTCACGGCGCGTTTTATTGCCTGATAGAGCCgg                 >  1:593516/1‑149 (MQ=255)
                                                                                                                                gTCGGCAAGTTTCGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGc                <  2:513652/149‑1 (MQ=255)
                                                                                                                                     cATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAg           <  1:1417156/149‑1 (MQ=255)
                                                                                                                                       tgttgtGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGACGCCAATCAGTAACGGCGCGGTTTATCGACTGATTGAGCAGGCCAAAggg         >  2:10957/1‑149 (MQ=255)
                                                                                                                                         ttgtGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAGTCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAGGgag       >  2:393692/1‑149 (MQ=255)
                                                                                                                                         ttgtGTCAGGTAATGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGAGAACTATGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGATTTATCGACAGATTGAGCAGGGCCAACGgag       >  1:77644/1‑149 (MQ=255)
                                                                                                                                          tgggtCAGGTACTGGCTGCCATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAGGgaga      <  2:48658/146‑1 (MQ=255)
                                                                                                                                              tCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGGTCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGCTTTATCGACTCATTGAGCAGGCCCAAGGGAGAGgat  >  1:898036/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
AAGTCCCAGCAGAAGCGGAAGAGATGGGAGGTGAAAAAGTGGAAGAAGTGGCAGCAGAGCAAGCGCCCAACGATGAAAAAGAGATGGGAGCTAAGGAGTGGGTTAACCCGAACCCGAAAAAATACAGCGTCGGCATGTTGTGTCAGGTACTGGCTGCAATGGGCGGC‑GAATACCTGGGCAACAACGCAGGGTTACAGCAGAAGATCACTGTGTTGGATAACGACGGCAACCGGAAGCCAATCAGTAACGGCGCGTTTTATCGACTGATTGAGCAGGCCAAAGGGAGAGGAT  >  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: ‑