Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,950,972 T→G 5.5% intergenic (‑122/+35) fliL ← / ← fliK Flagellar protein FliL/Flagellar hook‑length control protein FliK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,950,9720TG 5.5% 237.4 / 2.6 66intergenic (‑122/+35)fliL/fliKFlagellar protein FliL/Flagellar hook‑length control protein FliK
Reads supporting (aligned to +/- strand):  ref base T (21/41);  new base G (2/2);  total (23/43)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

CACCGAAAGGCCGATCGCCAGCAACAGGCCGACGACAGCCAGCAGGATCAGCTTGAGTTTGCCTTTAGTGGCGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGAAAGGTTGCTGCCCTGCGCCTGCCCCTGCTGTTGTTGCTGCCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCC  >  AE015451/4950749‑4951205
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
cACCGAAAGGCCGATCGCCAGCAACAGGCCGACGACAGCCAGCAGGATCAGCTTGAGTTTGCCTTTAGTGGCGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaaga                                                                                                                                                                                                                               <  1:669471/236‑1 (MQ=255)
     aaaGGCCGATCGCCAGCAACAGGCCGACGACAGCCAGCAGGATCAGCTTGAGTTTGCCTTTAGTGGCGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagcca                                                                                                                                                                                                                          <  1:176004/236‑1 (MQ=255)
             aTCGCCAGCAACAGGCCGACGACAGCCAGCAGGATCAGCTTGAGTTTGCCTTTAGTGGCGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAg                                                                                                                                                                                                                    <  2:826804/234‑1 (MQ=255)
             aTCGCCAGCAACAGGCCGACGACAGCCAGCAGGATCAGCTTGAGTTTGCCTTTAGTGGCGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGc                                                                                                                                                                                                                  <  1:1228316/236‑1 (MQ=255)
              tCGCCAGCAACAGGCCGACGACAGCCAGCAGGATCAGCTTGAGTTTGCCTTTAGTGGCGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGc                                                                                                                                                                                                                  <  1:550462/235‑1 (MQ=255)
                    gCAACAGGCCGACGACAGCCAGCAGGATCAGCTTGAGTTTGCCTTTAGTGGCGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAg                                                                                                                                                                                                            <  2:770748/235‑1 (MQ=255)
                        cAGGCCGACGACAGCCAGCAGGATCAGCTTGAGTTTGCCTTTAGTGGCGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATc                                                                                                                                                                                                       <  1:557223/236‑1 (MQ=255)
                               acgacAGCCAGCAGGATCAGCTTGAGTTTGCCTTTAGTGGCGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                                                                                                                <  1:288455/236‑1 (MQ=255)
                                          cagGATCAGCTTGAGTTTGCCTTTAGTGGCGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGAcc                                                                                                                                                                                      <  1:1279394/235‑1 (MQ=255)
                                                   cTTGAGTTTGCCTTTAGTGGCGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGt                                                                                                                                                                               <  2:1128044/233‑1 (MQ=255)
                                                     tGAGTTTGTCTTTAGTGGCGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCAAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCg                                                                                                                                                                          <  1:1224198/236‑1 (MQ=255)
                                                           tGCCTTTAGTGGCGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAAGCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAg                                                                                                                                                                                                                    <  1:291700/188‑1 (MQ=255)
                                                                       cGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTAt                                                                                                                                                                                                        <  1:537520/188‑1 (MQ=255)
                                                                                         cGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCaa                                                                                                                                                                                                               >  2:724498/1‑163 (MQ=255)
                                                                                         cGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATc                                                                                                                                       >  1:317020/1‑235 (MQ=255)
                                                                                              ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATc                                                                                                                                                                                                       <  1:491357/166‑1 (MQ=255)
                                                                                              ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGAcc                                                                                                                                   <  2:1185996/234‑1 (MQ=255)
                                                                                              ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCg                                                                                                                                  <  1:521601/235‑1 (MQ=255)
                                                                                                   cATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGGCGCCGACCACCACCTGCTGCGCCAAAg                                                                                                                                                    >  2:192783/1‑212 (MQ=255)
                                                                                                       ccAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAg                                                                                                                                                                                                     <  1:68235/159‑1 (MQ=255)
                                                                                                        cAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCt                                                                                                                                          >  2:737484/1‑217 (MQ=255)
                                                                                                         aataatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCTAACAg                                                                                                                                                                                                            >  2:932172/1‑150 (MQ=255)
                                                                                                           taatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACa                                                                                                                                                                       <  1:534634/185‑1 (MQ=255)
                                                                                                               ccGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCt                                                                                                                   <  1:8400/233‑1 (MQ=255)
                                                                                                                 gtcgtcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCaa                                                                                                                                                      <  1:1003468/196‑1 (MQ=255)
                                                                                                                 ggcggccATGGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAGGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaaga                                                                                                                                                                                                                               <  2:277939/118‑1 (MQ=255)
                                                                                                                          cGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGgcg                                                                                                                                            >  2:272576/1‑197 (MQ=255)
                                                                                                                           gggggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCg                                                                                                                                  <  1:783074/206‑1 (MQ=255)
                                                                                                                             gggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaaga                                                                                                                                                                                                                               >  2:1132481/1‑111 (MQ=255)
                                                                                                                             gggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCgaaga                                                                                                                                                                                                                               <  1:1132486/111‑1 (MQ=255)
                                                                                                                              ggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATGGGCCGCACTGTCGCCGACAACCACCTGCTGCTc                                                                                                                                                         >  2:607576/1‑180 (MQ=255)
                                                                                                                                ttCAAAGGAGAATGGCTTGGCTTTAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCg                                                                                                                                  <  1:315641/201‑1 (MQ=255)
                                                                                                                                  cAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAc                                                                                             >  1:1146365/1‑236 (MQ=255)
                                                                                                                                   aaaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGt                                                                                                                                                                               <  2:1061744/153‑1 (MQ=255)
                                                                                                                                   aaaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGt                                                                                                                                                                               >  1:1061749/1‑153 (MQ=255)
                                                                                                                                            tggcttggctTGAGCAAGTGTTATGCCAGATTTTGCAGGGGAGGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCt                                                                                                                   <  1:74016/204‑1 (MQ=255)
                                                                                                                                             ggcttggcttGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCa                                                                                                                                                                                                                     >  1:821344/1‑105 (MQ=255)
                                                                                                                                             ggcttggcttGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCa                                                                                                                                                                                                                     <  2:821339/105‑1 (MQ=255)
                                                                                                                                             ggcttggcttGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGa                                                                                                                                     >  1:389560/1‑185 (MQ=255)
                                                                                                                                             ggcttggcttGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGc                                                                                                            <  1:1278449/210‑1 (MQ=255)
                                                                                                                                                 tggcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGcc                                                                                                                                                 <  1:939829/169‑1 (MQ=255)
                                                                                                                                                    cttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCTAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTCGTAATCGACCGCACTGTCGCCGACAaccac                                                                                                                                                                  >  2:290527/1‑149 (MQ=255)
                                                                                                                                                      tgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGAAAGGTTGCTGCCCt                                                                         <  1:272580/236‑1 (MQ=255)
                                                                                                                                                        aGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCg                                                                                                                                  <  2:695694/177‑1 (MQ=255)
                                                                                                                                                           aaGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCTAAAGGCCGGGCGCTATCGACCg                                                                                                                                  >  1:627534/1‑174 (MQ=255)
                                                                                                                                                           aaGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCTAAAGGCCGGGCGCTATCGACCg                                                                                                                                  <  2:627529/174‑1 (MQ=255)
                                                                                                                                                           aaGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGAcc                                                                                                                                   >  1:1248694/1‑173 (MQ=255)
                                                                                                                                                                ttATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcccc                                                                                                  >  2:612618/1‑201 (MQ=255)
                                                                                                                                                                ttATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGAAAGGTTGCTGCCCTGCGCCTGccc                                                               >  1:595850/1‑236 (MQ=255)
                                                                                                                                                                 tATGCCAGATTTTGCAGGGGAGGGATTCCAGCAGGTTTTGCGGTGCCGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGAAAgg                                                                                   <  1:178019/215‑1 (MQ=255)
                                                                                                                                                                           tttGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGAAAGGTTGCTGCCCTGCGCCTGCCCCTGCtgttg                                                      >  2:607617/1‑234 (MQ=255)
                                                                                                                                                                                          ttCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGccc                                                                                                                                                                                                                                      >  1:690361/1‑43 (MQ=255)
                                                                                                                                                                                          ttCCAGCAGGTTTTGCGGGGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCCCTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGcccc                                                                                                                        <  2:427936/153‑1 (MQ=255)
                                                                                                                                                                                                cagGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGAAAGGTTGCTGCCCTGCGCCTGCCCCTGCTGTTGTTGCTGCCCGCGCGACTGGTCa                                <  1:192786/235‑1 (MQ=255)
                                                                                                                                                                                                     tttGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAATGCCGGGCGCTATCGACCg                                                                                                                                  >  2:1243099/1‑132 (MQ=255)
                                                                                                                                                                                                     tttGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGAAAGGTTGCTGCCCTGCGcc                                                                    <  1:814613/194‑1 (MQ=255)
                                                                                                                                                                                                       gcctgggCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGACAGGTTGCTGCCCTGCGCCTg                                                                  >  1:838478/7‑194 (MQ=255)
                                                                                                                                                                                                         cGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCgg                                                                                        <  1:724503/170‑1 (MQ=255)
                                                                                                                                                                                                              cTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCg                                                                                                                                  <  1:724618/123‑1 (MQ=255)
                                                                                                                                                                                                                 gCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCg                                                                                         >  2:1047548/1‑161 (MQ=255)
                                                                                                                                                                                                                       tcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCa                                                                                                                                                                                                                <  2:70749/36‑1 (MQ=11)
                                                                                                                                                                                                                        cgcgGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGAAAGGTTGCTGCCCTGCGCCTGccc                                                               <  1:122334/180‑1 (MQ=255)
                                                                                                                                                                                                                         gcgGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAcccc                                                                                          >  1:250263/1‑152 (MQ=255)
                                                                                                                                                                                                                            ggCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGt                                                                                                                                                                                                           <  2:452278/36‑1 (MQ=255)
                                                                                                                                                                                                                            ggCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCg                                                                                                                                  <  2:151139/109‑1 (MQ=255)
                                                                                                                                                                                                                             gCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGTCAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGAAAGGTTGCTGCCCTGCGCCTGCCCCTGCTGTTGTTGCTGCCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAAcccc  <  1:737489/236‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
CACCGAAAGGCCGATCGCCAGCAACAGGCCGACGACAGCCAGCAGGATCAGCTTGAGTTTGCCTTTAGTGGCGGGGTCTTTGACTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCCCGGAAAGGTTGCTGCCCTGCGCCTGCCCCTGCTGTTGTTGCTGCCCGCGCGACTGGTCAGCCACGTTGACATCCGGCTGCGCCAACCCC  >  AE015451/4950749‑4951205

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