Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,544,339 C→T 100% V189V (GTC→GTT hlgC → gamma hemolysin component C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,544,3390CT95.6% 175.6 / 1.2 70V189V (GTC→GTThlgCgamma hemolysin component C
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (41/24);  total (43/27)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.41e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTCAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGATTATTTCGTTCCAGACAGTGAGTTACCAC  >  CP000730/2544202‑2544463
                                                                                                                                          |                                                                                                                            
atCGGTGGTAGTTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTaa                                                                                                                            >  2:514077/1‑140 (MQ=255)
atCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTGTCTGGGGCGTTaa                                                                                                                            >  2:744559/1‑140 (MQ=255)
atCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTATAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTaa                                                                                                                            >  2:66587/1‑140 (MQ=255)
 tCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTAAGTAAGTGAAGTAGAACAACAACACTCAAAAAGTGTTTTATGGGACGTTaaa                                                                                                                           >  2:31044/1‑140 (MQ=255)
          aTTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAattca                                                                                                                  >  2:451060/1‑140 (MQ=255)
          aTTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAAAAACAAAACTCAAAAAGGGTTTTATGGGGCGTTAAAGCGAattaa                                                                                                                  >  2:20823/1‑138 (MQ=255)
            ttCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACCAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAattcatt                                                                                                                >  2:296749/1‑140 (MQ=255)
            ttCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATGAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGCTTTATGGGGCGTTACAGCGAattcatt                                                                                                                >  2:442444/1‑140 (MQ=255)
             tCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAattcattc                                                                                                               >  2:715276/1‑140 (MQ=255)
             tCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAattcattc                                                                                                               >  2:729093/1‑140 (MQ=255)
                  tCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCAc                                                                                                          >  2:752121/1‑140 (MQ=255)
                  tCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCAc                                                                                                          >  1:194888/1‑140 (MQ=255)
                  tCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAAAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCAc                                                                                                          >  1:194937/1‑140 (MQ=255)
                    aGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAATTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTg                                                                                                        >  2:169903/1‑140 (MQ=255)
                    aGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTg                                                                                                        <  2:663571/140‑1 (MQ=255)
                       cccATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAAt                                                                                                     >  2:362631/1‑140 (MQ=255)
                         cATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCa                                                                                                   >  1:759357/1‑140 (MQ=255)
                            tacccGGTAGTAATGGAACAATTAACTATTCAAAATCAATCACCAAAACAAAACCTAAATATATAAAACAAATAAAAAAATAAAAATAAAAAAGTATTTTATGGCGCGTTAAAGCTAATTCATTCGCCACTGAATCAGGt                                                                                                <  2:283898/136‑1 (MQ=255)
                                 ggtggtAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAATGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCaaaa                                                                                           >  1:102073/1‑140 (MQ=255)
                                    ggtAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAATTCCTTTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATc                                                                                        <  1:590032/140‑1 (MQ=255)
                                        aTGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGcc                                                                                    <  1:567271/140‑1 (MQ=255)
                                        aTGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGcc                                                                                    >  1:96729/1‑140 (MQ=255)
                                          ggATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGAGAATTAATACGCCCCTGAATCAGAACAAAAagcagcata                                                                                  >  2:716535/1‑137 (MQ=255)
                                                  tAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCg                                                                          >  1:454680/1‑140 (MQ=255)
                                                   aaCTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGa                                                                         <  1:598544/140‑1 (MQ=255)
                                                   aaCTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGAATAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGa                                                                         >  1:589925/1‑140 (MQ=255)
                                                       aatCTAAATCGATTAACAAAACACAACAAAATTATGTAAGTGAAGTAGAAAAAAAAAAAAAAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGattta                                                                     <  2:446256/138‑1 (MQ=255)
                                                       aTTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAATAAGTGTTTTATTGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGattta                                                                     <  2:428890/140‑1 (MQ=255)
                                                        ttCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATACGCCACTGAATCAGGTCAAAAATCAGCCTTTGATGGCGatttat                                                                    >  2:171786/1‑140 (MQ=255)
                                                         tCTAAATCGATTATCTTTCCACAACAAAATTATGTAAGTGAAGTATAATTACTTCACTCTTAAAGTGTTTTATGGGGCGTTAAAGTGAATTTATTCGTCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGatttatt                                                                   <  1:670580/140‑1 (MQ=255)
                                                                cGATTAGCTACACAGAACCAAGTTATACAAGAGAAGNAGAACAACAACCCTCAAGAAGCATCTAACGGCGGTGTAAAGCTAATTCATTGGACAAAAAATCAAAACAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGc                                                            <  2:1933/140‑1 (MQ=255)
                                                                 gATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCt                                                           <  1:677040/140‑1 (MQ=255)
                                                                  aTTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTa                                                          <  1:339599/140‑1 (MQ=255)
                                                                  aTTAGCTAGACACAACTAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTa                                                          <  2:392039/140‑1 (MQ=255)
                                                                    tAGCTATACAAAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTCAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACa                                                        <  1:427427/140‑1 (MQ=255)
                                                                      gCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGCATTACTTCGCCACTGCATCCGGACACAACTCAGACTTGGATAGGGATTTAATTGTAGGCTACaaa                                                      >  2:522635/1‑140 (MQ=255)
                                                                      gCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCAATGCATCAGGTCAGCAATCGGCCTTTGAGAGCGAGTTATCTGTAGGCTACaag                                                      >  2:522654/1‑139 (MQ=255)
                                                                        tataCACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTCAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAAcc                                                    >  2:289957/1‑140 (MQ=255)
                                                                         ataCACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAACGAATTCATTCGCCACTAAATAAGGTCAAAAAGCAGCCATTGATAGCGGTTTATTTGTAGCCTACAAACCt                                                   >  1:235177/1‑140 (MQ=255)
                                                                         ataCACAACAAAATTATGTAAGTGAACTAGAACAACATAAAGCAAACAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGACACTGAATCAGGTCACAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCt                                                   <  2:728629/140‑1 (MQ=255)
                                                                            cacaACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCAt                                                >  1:597636/1‑140 (MQ=255)
                                                                            cacaACAAAATTATGTAAGAGAAATAGAACAACAAAACTCAAAAATTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCAt                                                <  1:476785/140‑1 (MQ=255)
                                                                              caacaaAATTATGTAAGTGAAGTAGACCAAGTCGACGCTGTAGGTGTTTAATGGGGCGTCAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGATACacagctactag                                              <  2:343443/140‑11 (MQ=255)
                                                                               aacaaAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACATCATAGt                                             >  2:119152/1‑140 (MQ=255)
                                                                                 caaAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGATAAAGCGAATTCATTCGCCACTGACTCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTAATAGTaa                                           >  1:392025/1‑140 (MQ=255)
                                                                                    aaTTATGTAAGTGAAGTAGTACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCAAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGAGTTATTTGTAGGCTACAAACCGCATAGTAAAGa                                        >  2:425781/1‑140 (MQ=255)
                                                                                    aaTTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGa                                        >  1:679810/1‑140 (MQ=255)
                                                                                          ggaatttaaggaaaacaacTAAACTCAAAAAGTGTTTTATGGGGAGATAAAGCGAATTAATTCGCCACTGAATCAGGTCAAAAAACAGCCTAAGATCGTCATTTCTTTTAAGGCTACAAACCACAAAGTACAGACCCTag                                  >  2:728232/14‑140 (MQ=255)
                                                                                          gTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCAag                                  >  1:538534/1‑140 (MQ=255)
                                                                                          gTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCACATAGCAAAGATCCTag                                  >  2:2801/1‑140 (MQ=255)
                                                                                               tGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGatt                             <  1:264425/140‑1 (MQ=255)
                                                                                                gAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAATCCTCATAGTAAAGATCCTAGAGatta                            <  2:454680/140‑1 (MQ=255)
                                                                                                 aaGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCGGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATACTAGAGcttat                           >  2:382596/1‑140 (MQ=255)
                                                                                                 aaGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGattat                           <  1:581570/140‑1 (MQ=255)
                                                                                                 aaGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATACTAGAGattat                           >  2:458866/1‑140 (MQ=255)
                                                                                                   gTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGATTAttt                         <  1:451060/140‑1 (MQ=255)
                                                                                                     aGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGATTAtttgg                       >  1:166463/1‑138 (MQ=255)
                                                                                                     aGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGATTATTTCg                       >  1:276228/1‑140 (MQ=255)
                                                                                                      gAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATACTAGAGATTATTTCGt                      >  2:529140/1‑140 (MQ=255)
                                                                                                      gAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTAGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGATTATTTCGt                      >  2:560162/1‑140 (MQ=255)
                                                                                                       aacaacAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGATTATTTCGtt                     >  2:628253/1‑140 (MQ=255)
                                                                                                       aacaacAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGATTATTTCGtt                     >  2:70081/1‑140 (MQ=255)
                                                                                                       aacaacAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGATTATTTCGtt                     >  1:724876/1‑140 (MQ=255)
                                                                                                               aaCTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGATTATTTCGTTCCAGACAg             <  1:433985/140‑1 (MQ=255)
                                                                                                                    aaaaaGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGATTATTTCGTTCCAGACAGTGAGt        <  1:611545/140‑1 (MQ=255)
                                                                                                                     aaaaGTGTTTGATGGGTCGTTAACGCGAAGGCATCCGCCACTGCTTGCTTACAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGATTATTTCGTTCCAGACAGTGAGtt       <  1:171786/140‑1 (MQ=255)
                                                                                                                       aaGTGTTTTATGGGGCTTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTAAAAACCTCATAGTAAAGATCCTAGAGATTATTTCGTTCCAGACAGTGAGTTac     <  2:60172/140‑1 (MQ=255)
                                                                                                                       aaGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGATTATTTCGTTCCAGACCGTGAGTTac     <  1:557378/140‑1 (MQ=255)
                                                                                                                       aaGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGATTATTTCGTTCCAGACAGTGAGTTac     <  1:331575/140‑1 (MQ=255)
                                                                                                                          tgtTTTATGGGGCGTTAATGCCAATTCATTCGTCACTGAATCAGGTCTAATATCAGTCTTTGATAGCTATTTATTTGTAGGCTACAAACCTTATAGTAAAGATCCTAGAGATTATTTCGTTCCAGACAGTGAGTTaccac  <  2:423670/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                            
ATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTCAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCTCATAGTAAAGATCCTAGAGATTATTTCGTTCCAGACAGTGAGTTACCAC  >  CP000730/2544202‑2544463

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: