Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,193,041 T→C 7.0% K218R (AAA→AGA)  ymfE ← e14 prophage; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,193,0410TC 7.0% 282.6 / 9.2 85K218R (AAA→AGA) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (39/40);  new base C (2/4);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTAT  >  CP009273/1192893‑1193186
                                                                                                                                                    |                                                                                                                                                 
cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAtttt                                                                                                                                                  >  1:2199733/1‑150 (MQ=255)
    cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCg                                                                                                                                              >  1:4042996/1‑150 (MQ=255)
    cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCg                                                                                                                                              >  1:4043010/1‑150 (MQ=255)
       aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgact                                                                                                                                           <  1:1532743/150‑1 (MQ=255)
          tatCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttga                                                                                                                                        >  1:4073999/1‑150 (MQ=255)
           atCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttgac                                                                                                                                       >  2:232760/1‑150 (MQ=255)
               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTc                                                                                                                                   >  2:820221/1‑150 (MQ=255)
                   ggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGtt                                                                                                                               >  2:1268897/1‑150 (MQ=255)
                     cATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCt                                                                                                                             >  1:922560/1‑150 (MQ=255)
                       tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTgg                                                                                                                           <  2:2873451/150‑1 (MQ=255)
                          tATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCt                                                                                                                        >  2:3071739/1‑150 (MQ=255)
                          tATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCt                                                                                                                        >  2:109445/1‑150 (MQ=255)
                          tATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCt                                                                                                                        >  2:109431/1‑150 (MQ=255)
                           aTGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTCTTCCGACTTGACTTTCAGTTCTGGGCtt                                                                                                                       <  2:1414409/150‑1 (MQ=255)
                           aTGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTCTTCCGACTTGACTTTCAGTTCTGGGCtt                                                                                                                       <  2:917742/150‑1 (MQ=255)
                           aTGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTCTTCCGACTTGACTTTCAGTTCTGGGCtt                                                                                                                       <  2:917756/150‑1 (MQ=255)
                           aTGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTCTTCCGACTTGACTTTCAGTTCTGGGCtt                                                                                                                       <  2:917948/150‑1 (MQ=255)
                              aGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGTTTTTTACGACTTGACTTTCAGTTCTGGGCTTAAc                                                                                                                    >  2:2843117/1‑150 (MQ=255)
                               ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACt                                                                                                                   >  2:480652/1‑150 (MQ=255)
                               ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACt                                                                                                                   <  1:1820115/150‑1 (MQ=255)
                                      ctgcATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaag                                                                                                            >  2:2062993/1‑150 (MQ=255)
                                      ctgcATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaag                                                                                                            <  1:3512468/150‑1 (MQ=255)
                                      ctgcATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaag                                                                                                            <  1:3512462/150‑1 (MQ=255)
                                       tgcATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagc                                                                                                           >  2:2275606/1‑150 (MQ=255)
                                          aTTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagcaaa                                                                                                        <  1:3374642/150‑1 (MQ=255)
                                           ttCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagcaaag                                                                                                       <  2:73768/150‑1 (MQ=255)
                                           ttCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagcaaag                                                                                                       <  2:73665/150‑1 (MQ=255)
                                           ttCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagcaaag                                                                                                       <  2:432573/150‑1 (MQ=255)
                                           ttCTCTGAGAGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagcaaag                                                                                                       <  2:73877/150‑1 (MQ=255)
                                            tctcTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTCTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGa                                                                                                      >  2:1794596/1‑150 (MQ=255)
                                              tctGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGata                                                                                                    <  2:3254289/150‑1 (MQ=255)
                                                       atatTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactc                                                                                           <  2:2457710/150‑1 (MQ=255)
                                                         atTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactcac                                                                                         >  1:2288815/1‑150 (MQ=255)
                                                              cTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTgaga                                                                                    >  2:2813071/1‑150 (MQ=255)
                                                              cTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTgaga                                                                                    >  2:2192566/1‑150 (MQ=255)
                                                              cTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTgaga                                                                                    >  2:2191972/1‑150 (MQ=255)
                                                                     cATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATaa                                                                             <  1:2485586/150‑1 (MQ=255)
                                                                      atatTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTATTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATaaa                                                                            >  2:3427519/1‑150 (MQ=255)
                                                                      atatTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATaaa                                                                            <  1:1500822/150‑1 (MQ=255)
                                                                      atatTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATaaa                                                                            >  1:2637044/1‑150 (MQ=255)
                                                                      atatTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTCTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATaaa                                                                            >  1:553423/1‑150 (MQ=255)
                                                                       tattatAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATaaaa                                                                           <  2:1737123/150‑1 (MQ=255)
                                                                          tataGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTa                                                                        >  1:2037689/1‑150 (MQ=255)
                                                                           ataGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTAc                                                                       >  2:474459/1‑150 (MQ=255)
                                                                           ataGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTAc                                                                       >  2:2557501/1‑150 (MQ=255)
                                                                            taGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACa                                                                      <  1:3426660/150‑1 (MQ=255)
                                                                            taGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGAGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACa                                                                      <  2:2375115/150‑1 (MQ=255)
                                                                              gTTCAAATCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACaaa                                                                    <  2:2724332/150‑1 (MQ=255)
                                                                              gTTCAAAGCAATTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACaaa                                                                    <  2:2724327/150‑1 (MQ=255)
                                                                              gTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACaaa                                                                    >  2:2945644/1‑150 (MQ=255)
                                                                               ttCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACaaaa                                                                   >  2:713934/1‑150 (MQ=255)
                                                                                 cAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACaaaaaa                                                                 <  2:3290906/150‑1 (MQ=255)
                                                                                  aaaGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAg                                                                >  1:3271178/1‑150 (MQ=255)
                                                                                   aaGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGa                                                               <  2:2762484/150‑1 (MQ=255)
                                                                                      cAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATa                                                            >  2:3792111/1‑150 (MQ=255)
                                                                                       aaGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATaa                                                           >  2:2024968/1‑150 (MQ=255)
                                                                                           aCCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGAc                                                       <  2:3304924/150‑1 (MQ=255)
                                                                                            ccATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACa                                                      >  1:3051081/1‑150 (MQ=255)
                                                                                             cATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTTACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAg                                                     >  2:1745825/1‑150 (MQ=255)
                                                                                             cATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAg                                                     >  2:1745817/1‑150 (MQ=255)
                                                                                               ttCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCa                                                   >  2:1623342/1‑150 (MQ=255)
                                                                                                 cAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCaaa                                                 <  1:505637/150‑1 (MQ=255)
                                                                                                 cAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCaaa                                                 <  1:505613/150‑1 (MQ=255)
                                                                                                    ttAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGaa                                              <  2:3218841/150‑1 (MQ=255)
                                                                                                    ttAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGaa                                              >  2:2158306/1‑150 (MQ=255)
                                                                                                            ttttttAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAAc                                      <  2:3543518/149‑1 (MQ=255)
                                                                                                            aTTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAAc                                      <  2:3279153/150‑1 (MQ=255)
                                                                                                            aTTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAAc                                      <  2:3279186/150‑1 (MQ=255)
                                                                                                               tttAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATg                                   <  1:3822300/150‑1 (MQ=255)
                                                                                                                ttAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGa                                  >  2:1964920/1‑150 (MQ=255)
                                                                                                                 tAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGaa                                 >  2:1483896/1‑150 (MQ=255)
                                                                                                                  aGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGaaa                                >  2:554688/1‑150 (MQ=255)
                                                                                                                  aGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGaaa                                >  1:262079/1‑150 (MQ=255)
                                                                                                                  aGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGaaa                                >  1:311705/1‑150 (MQ=255)
                                                                                                                    tGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGaaaaa                              <  1:1670805/150‑1 (MQ=255)
                                                                                                                     gAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAat                             <  2:4043010/150‑1 (MQ=255)
                                                                                                                     gAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAat                             <  2:4042996/150‑1 (MQ=255)
                                                                                                                     gAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAat                             <  1:1153183/150‑1 (MQ=255)
                                                                                                                        aaaTTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATc                          <  2:3377690/150‑1 (MQ=255)
                                                                                                                         aaTTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATcc                         <  1:2063064/150‑1 (MQ=255)
                                                                                                                            tCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCAtt                      >  2:3336599/1‑150 (MQ=255)
                                                                                                                                  attttatCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCAtttgtttg                <  1:4129362/150‑1 (MQ=255)
                                                                                                                                   ttttatCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGc               <  1:871897/150‑1 (MQ=255)
                                                                                                                                   ttttatCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGc               <  1:1864954/150‑1 (MQ=255)
                                                                                                                                           tCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAg       <  1:2888310/150‑1 (MQ=255)
                                                                                                                                                gATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGattat  <  2:2655795/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTAT  >  CP009273/1192893‑1193186

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: ‑