Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 536,758 C→T 91.7% P202L (CCA→CTA)  rplY → ribosomal protein L25

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730536,7580CT91.7% 167.2 / 4.0 68P202L (CCA→CTA) rplYribosomal protein L25
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base T (21/39);  total (25/43)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

TTGCTGATGTTAAAGTAACTGGCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACCAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATT  >  CP000730/536620‑536894
                                                                                                                                          |                                                                                                                                        
ttGCTGATGTTAAAGTAACTGGCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACGGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACGAg                                                                                                                                        >  2:198264/1‑141 (MQ=255)
 tGCTGATGTTAAAGTAACTGGCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGa                                                                                                                                       <  1:641283/141‑1 (MQ=255)
 tGCTGATGTTAAAGTAACTGGCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGa                                                                                                                                       <  1:570161/141‑1 (MQ=255)
 tGCTGATGTTAAAGTAACTGGCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGACGAACAACAAACTGAAGAACTAGa                                                                                                                                       <  1:485506/141‑1 (MQ=255)
  gCTGATGTTAAAGTAACTGGCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGaa                                                                                                                                      <  2:567901/141‑1 (MQ=255)
  gCTGATGTTAAAGTAACTGGCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGaa                                                                                                                                      >  1:207701/1‑141 (MQ=255)
  gCTGATGTTAAAGTAACTGGCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCGATGGAAGGCGAACAACAAACTGAAGAACCAGaa                                                                                                                                      >  2:74493/1‑141 (MQ=255)
  gCTGATGTTAAAGTAACTGGCGACTTCAAAATCGAAAACGATTCAGCGGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGaa                                                                                                                                      <  1:273329/141‑1 (MQ=255)
       tGTTAAAGTAACTGGCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTACTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAgttgt                                                                                                                                 >  2:624880/1‑141 (MQ=255)
           aaaGTAACTGGCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGc                                                                                                                             <  2:339363/141‑1 (MQ=255)
            aaGTAACTGGCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCg                                                                                                                            <  1:309202/141‑1 (MQ=255)
               tAACTGGCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGaaa                                                                                                                         >  2:560914/1‑141 (MQ=255)
                 aCTGGCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGc                                                                                                                       <  2:484024/141‑1 (MQ=255)
                    ggCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCaaa                                                                                                                    >  1:368994/1‑141 (MQ=255)
                                   gAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGaaaaa                                                                                                     >  1:9301/1‑141 (MQ=255)
                                   gAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCCGATTACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGaaaaa                                                                                                     >  1:453634/1‑141 (MQ=255)
                                       aCGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTg                                                                                                 <  1:653044/141‑1 (MQ=255)
                                       aCGATTCAGCTGAATCAGTAGCAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTg                                                                                                 <  1:14810/141‑1 (MQ=255)
                                        cGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGa                                                                                                >  1:505016/1‑141 (MQ=255)
                                            tCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAagag                                                                                            <  1:459467/141‑1 (MQ=255)
                                            tCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCCATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAagag                                                                                            <  2:202867/141‑1 (MQ=255)
                                             cAGCTGAATCAGTAGGAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGt                                                                                           >  1:334814/1‑141 (MQ=255)
                                               gCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACCAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTaa                                                                                         >  1:60126/1‑141 (MQ=255)
                                                cTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAAt                                                                                        <  1:51781/141‑1 (MQ=255)
                                                 tGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAAtt                                                                                       >  2:575320/1‑141 (MQ=255)
                                                   aaacaaaagaaaCAGAAATTGCTGCAACTGAAGAACCAACTGAAGAAGAGAACGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAAAAAAAAAAAAAAAAACTGAAGAGTAAtttt                                                                                     <  2:205965/131‑1 (MQ=255)
                                                    aTCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAGCTGAAGAGTAATTTTa                                                                                    <  1:331704/141‑1 (MQ=255)
                                                          agtaACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGt                                                                              <  1:412287/141‑1 (MQ=255)
                                                          agtaACAGTAGTTGCTCCAACTGAAGAAACAACAGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTAAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAAAAAAAAAAAAAACTGAAGAGTAATTTTAATCTGt                                                                              <  2:43397/141‑1 (MQ=255)
                                                           gtaACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGtt                                                                             >  1:42899/1‑141 (MQ=255)
                                                             aaCAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTAc                                                                           >  2:180270/1‑141 (MQ=255)
                                                                 aaagtAGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGACCAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAAAAAAAAAAAAAAAAACTGAAGAGTAATTTTAATCTGTTACATTa                                                                       <  2:8267/139‑1 (MQ=255)
                                                                   aggtGCTCTAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGATTAATTTTAATCTGTTACATTaaa                                                                     <  2:573997/138‑1 (MQ=255)
                                                                    gtTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAg                                                                    >  1:226149/1‑141 (MQ=255)
                                                                       gCTCCAACTGAAGAACGAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGttt                                                                 <  2:282361/141‑1 (MQ=255)
                                                                       gCTCCAACTGAAGAACCAACTGAAGAAGAAAACGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGttt                                                                 <  1:112519/141‑1 (MQ=255)
                                                                        cTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGtttt                                                                >  1:429238/1‑141 (MQ=255)
                                                                             aatataaaaccaactgaagaaGAAAGCGAAGCTATGGAAGGCCAACAACAAACAGAAGAACTAGAAGTTGTTGGCGAAAGCAAAAAAAAAAAAAAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATAc                                                           <  2:477019/134‑1 (MQ=255)
                                                                                  agaaccaactgaagaaGAAATCGTAGCTATGGAAGGCGAACAACAAACTGAAGAACTACAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTACACTTTGt                                                      <  1:418992/141‑1 (MQ=255)
                                                                                  agaaccaactgaagaaGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGt                                                      >  1:601687/1‑141 (MQ=255)
                                                                                    aaccaactgaagaaGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAGGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTGTAATTTGttt                                                    >  2:264540/1‑141 (MQ=255)
                                                                                    aaccaactgaagaaGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGttt                                                    >  2:502045/1‑141 (MQ=255)
                                                                                        agctgaagaaGAAATCGAAGCTATGGAAGGCGAACAACAAACGGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACa                                                <  2:282186/139‑1 (MQ=255)
                                                                                        aactgaagaaGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACa                                                <  1:337322/141‑1 (MQ=255)
                                                                                         actgaagaaGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACaa                                               <  1:120036/141‑1 (MQ=255)
                                                                                           tgaagaaGAAATCTAAGCGAAGGAAGGCGAACAACAATGTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGc                                             <  2:620768/141‑1 (MQ=255)
                                                                                            gaagaagaaATCGAAGCTCTGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGAAAATTTAATCAGTTACATCAACGCTTCAATACTTTGTTTAACAAGcg                                            >  2:136722/1‑140 (MQ=255)
                                                                                            gaagaagaaATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCa                                            <  1:529537/141‑1 (MQ=255)
                                                                                             aagcagaaATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCAc                                           <  1:62649/141‑1 (MQ=255)
                                                                                             aagaagaaATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGATGAAAAAACTGAAGAGTAATTCTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCAc                                           <  1:297956/141‑1 (MQ=255)
                                                                                                aagaaATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACtgt                                        <  2:334814/141‑1 (MQ=255)
                                                                                                 agaaATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACtgtg                                       >  1:233705/1‑141 (MQ=255)
                                                                                                  gaaATGGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGc                                      <  2:484032/141‑1 (MQ=255)
                                                                                                     aTCGAAGGTAGGGAAGGCGAACAACAAACTGAAGAACCAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTa                                   <  1:149293/141‑1 (MQ=255)
                                                                                                     aTCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACCAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTa                                   <  1:510018/141‑1 (MQ=255)
                                                                                                       cGAAGCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTAtt                                 <  1:59331/141‑1 (MQ=255)
                                                                                                           gCTATGGAAGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTTTGCTTATTTTaa                             >  2:408272/1‑141 (MQ=255)
                                                                                                                  aaGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGc                      >  2:218505/1‑141 (MQ=255)
                                                                                                                  aaGGCGAACAACAAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGc                      <  2:505016/141‑1 (MQ=255)
                                                                                                                     gCGAACAACAAACCGAAGAACTAGAAGCTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATCAAAGCaag                   <  2:576942/141‑3 (MQ=255)
                                                                                                                         accacaaACTGAAGTATGAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATCAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGc               <  1:541936/138‑1 (MQ=255)
                                                                                                                         acaacaAACTGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGc               <  2:320580/141‑1 (MQ=255)
                                                                                                                         acaacaAACTGAAGAACCAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGc               <  1:167806/141‑1 (MQ=255)
                                                                                                                            acaaACGGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATGTTAATACAAGCATGGTGCttt            <  1:346984/141‑1 (MQ=255)
                                                                                                                               aaCTGAAAAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAATCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtg         <  2:239892/141‑1 (MQ=255)
                                                                                                                                  tGAAGAACTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTAGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTgtgt      <  2:490794/141‑1 (MQ=255)
                                                                                                                                     agaaCTAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttat   >  2:352734/1‑141 (MQ=255)
                                                                                                                                      gaaCCAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAACGTTTTTATACATTGTTTAACAAGTACTGTGCTTATTTTATTATAAGGAAGTTGCTTTTTGTGttatt  >  2:130337/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                        
TTGCTGATGTTAAAGTAACTGGCGACTTCAAAATCGAAAACGATTCAGCTGAATCAGTAGTAACAGTAGTTGCTCCAACTGAAGAACCAACTGAAGAAGAAATCGAAGCTATGGAAGGCGAACAACAAACTGAAGAACCAGAAGTTGTTGGCGAAAGCAAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATT  >  CP000730/536620‑536894

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