Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,009,187 G→T 93.4% intergenic (+90/‑232) mocA → / → ygfK CTP:molybdopterin cytidylyltransferase/putative Fe‑S subunit oxidoreductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,009,1870GT93.4% 773.0 / 39.8 251intergenic (+90/‑232)mocA/ygfKCTP:molybdopterin cytidylyltransferase/putative Fe‑S subunit oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base G (10/7);  new base T (121/113);  total (131/120)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 271 total aligned reads displayed.

AGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGG  >  CP009273/3009040‑3009333
                                                                                                                                                   |                                                                                                                                                  
aGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaaaagataa                                                                                                                                                  <  2:61843/150‑1 (MQ=255)
 gATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaaaagagaaa                                                                                                                                                 <  2:3724618/150‑1 (MQ=255)
   tatTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATaaaaa                                                                                                                                               >  1:1215157/1‑150 (MQ=255)
      tGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTg                                                                                                                                            >  1:3001600/1‑150 (MQ=255)
       gATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGACTCACTGATTTTGAGAAAAGATAAAAATTGa                                                                                                                                           >  1:3816927/1‑150 (MQ=255)
       gATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGa                                                                                                                                           >  1:2524953/1‑150 (MQ=255)
       gATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGa                                                                                                                                           >  2:2675813/1‑150 (MQ=255)
        aTACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGAt                                                                                                                                          <  2:1213553/150‑1 (MQ=255)
          aCTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATaa                                                                                                                                        >  2:3138295/1‑150 (MQ=255)
             ccGTACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATg                                                                                                                                     >  2:721/1‑150 (MQ=255)
                 acgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCt                                                                                                                                 >  2:170189/1‑150 (MQ=255)
                 acgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCt                                                                                                                                 <  1:1675732/150‑1 (MQ=255)
                   gacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTgg                                                                                                                               <  2:2861938/150‑1 (MQ=255)
                   gacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTgg                                                                                                                               >  1:228431/1‑150 (MQ=255)
                   gacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTgg                                                                                                                               >  1:1242238/1‑150 (MQ=255)
                    acTTTATTAAAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCTTTCAAATATTTAATCCTTCCTGCACTCATTATCATAATGATAAGTAATCAATGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGc                                                                                                                              >  2:2223618/1‑150 (MQ=255)
                           taCAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCaataata                                                                                                                       <  1:3610130/150‑1 (MQ=255)
                            aCAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCaataataa                                                                                                                      >  2:2828961/1‑150 (MQ=255)
                             cAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCaataataat                                                                                                                     >  1:623718/1‑150 (MQ=255)
                                 gAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTc                                                                                                                 <  1:1405081/150‑1 (MQ=255)
                                 gAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATATTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTc                                                                                                                 >  2:3387257/1‑150 (MQ=255)
                                  aaGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCa                                                                                                                >  2:1030462/1‑150 (MQ=255)
                                  aaGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCa                                                                                                                >  2:2924893/1‑150 (MQ=255)
                                  aaGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCa                                                                                                                >  2:1148503/1‑150 (MQ=255)
                                   aGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAg                                                                                                               <  2:1769162/150‑1 (MQ=255)
                                      aaaGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAg                                                                                                            <  1:2119928/150‑1 (MQ=255)
                                      aaaGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAg                                                                                                            <  1:4110639/150‑1 (MQ=255)
                                        aGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTg                                                                                                          >  1:3800148/1‑150 (MQ=255)
                                        aGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGAACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATATTCAGGAGTg                                                                                                          >  1:1190890/1‑150 (MQ=255)
                                          gTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAt                                                                                                        >  2:653892/1‑150 (MQ=255)
                                             taCTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAtttt                                                                                                     <  2:1525650/150‑1 (MQ=255)
                                                tGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAAt                                                                                                  >  1:2313851/1‑150 (MQ=255)
                                                 gAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAAtt                                                                                                 >  2:4097756/1‑150 (MQ=255)
                                                    aTTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATtctc                                                                                              >  1:4448257/1‑150 (MQ=255)
                                                      ttAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCAt                                                                                            <  1:3731433/150‑1 (MQ=255)
                                                          aaaaCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAg                                                                                        <  1:3527787/150‑1 (MQ=255)
                                                           aaaCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGc                                                                                       >  2:1638114/1‑150 (MQ=255)
                                                            aaCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCa                                                                                      >  2:2512555/1‑150 (MQ=255)
                                                            aaCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCa                                                                                      >  1:3881726/1‑150 (MQ=255)
                                                                tAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCtg                                                                                  >  1:2603080/1‑150 (MQ=255)
                                                                   ttAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATGATCATTTTCAGGAGTGATTTTAATTCTCATGTAGCACCTGTGa                                                                               >  2:273294/1‑150 (MQ=255)
                                                                   ttAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGa                                                                               <  1:1378039/150‑1 (MQ=255)
                                                                   ttAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGa                                                                               <  1:290674/150‑1 (MQ=255)
                                                                   ttAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGa                                                                               >  1:4221365/1‑150 (MQ=255)
                                                                   ttAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGa                                                                               <  1:3431905/150‑1 (MQ=255)
                                                                     aaaCGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGt                                                                             <  2:1925520/150‑1 (MQ=255)
                                                                       aCGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttc                                                                           >  1:4430289/1‑150 (MQ=255)
                                                                       aCGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttc                                                                           >  2:4528981/1‑150 (MQ=255)
                                                                       aCGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttc                                                                           <  2:3422324/150‑1 (MQ=255)
                                                                        cGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttca                                                                          >  2:826227/1‑150 (MQ=255)
                                                                         gATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttcat                                                                         <  1:552961/150‑1 (MQ=255)
                                                                         gATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttcat                                                                         >  2:3809797/1‑150 (MQ=255)
                                                                          aTTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttcatt                                                                        >  2:3910286/1‑150 (MQ=255)
                                                                              cGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTTAGTTCATTcaca                                                                    <  1:550556/150‑1 (MQ=255)
                                                                              cGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTcaca                                                                    >  1:653030/1‑150 (MQ=255)
                                                                               gTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACaa                                                                   >  1:3071670/1‑150 (MQ=255)
                                                                               gTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACaa                                                                   >  1:2508206/1‑150 (MQ=255)
                                                                                ttCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACaaa                                                                  >  2:2266760/1‑150 (MQ=255)
                                                                                 tCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGTCTGGCAATAATACTTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCAGTCACaaaa                                                                 >  2:3383822/1‑150 (MQ=255)
                                                                                   aaaTATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAg                                                               >  2:2690378/1‑150 (MQ=255)
                                                                                   aaaTATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAg                                                               >  1:4368175/1‑150 (MQ=255)
                                                                                       atTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCt                                                           >  1:2055572/1‑150 (MQ=255)
                                                                                           aaTCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAg                                                       <  2:2687074/150‑1 (MQ=255)
                                                                                            aTCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGt                                                      >  1:846416/1‑150 (MQ=255)
                                                                                            aTCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGt                                                      >  1:485375/1‑150 (MQ=255)
                                                                                            aTCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGt                                                      >  2:2562041/1‑150 (MQ=255)
                                                                                              ccTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGttt                                                    >  1:423239/1‑150 (MQ=255)
                                                                                               cTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTc                                                   <  1:4134998/150‑1 (MQ=255)
                                                                                               cTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTc                                                   >  2:4604951/1‑150 (MQ=255)
                                                                                                ttCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCg                                                  <  1:1614929/150‑1 (MQ=255)
                                                                                                   cGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAgc                                               >  1:143538/1‑150 (MQ=255)
                                                                                                   cGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAgc                                               >  1:143570/1‑150 (MQ=255)
                                                                                                        cTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCt                                          <  1:1460792/150‑1 (MQ=255)
                                                                                                         tCATTATCATAATGATACGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAACTGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTc                                         <  2:3658823/150‑1 (MQ=255)
                                                                                                             tATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAt                                     <  1:3144085/150‑1 (MQ=255)
                                                                                                              aTCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAtt                                    <  1:1894401/150‑1 (MQ=255)
                                                                                                                cATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAtttt                                  >  2:1301465/1‑150 (MQ=255)
                                                                                                                 aTAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAttttt                                 >  2:103267/1‑150 (MQ=255)
                                                                                                                   aaTGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTCAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAttttttt                               <  2:429187/150‑1 (MQ=255)
                                                                                                                    aTGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTta                              >  1:239027/1‑150 (MQ=255)
                                                                                                                     tGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTtat                             >  1:2320398/1‑150 (MQ=255)
                                                                                                                      gATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTtata                            <  1:1099348/150‑1 (MQ=255)
                                                                                                                        tAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATAtt                          >  2:3478843/1‑150 (MQ=255)
                                                                                                                         aaGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATAttt                         <  2:2095022/150‑1 (MQ=255)
                                                                                                                           ggAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCGCACATTTTTTTATATTTcc                       >  1:1072663/1‑150 (MQ=255)
                                                                                                                             aaTCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCg                     <  2:241877/150‑1 (MQ=255)
                                                                                                                                cACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCa                  <  1:2817557/150‑1 (MQ=255)
                                                                                                                                cACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCa                  <  1:462560/150‑1 (MQ=255)
                                                                                                                                cACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCa                  <  2:1331245/150‑1 (MQ=255)
                                                                                                                                 aCTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCaa                 >  2:4135742/1‑150 (MQ=255)
                                                                                                                                   tGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAAc               >  1:3152966/1‑150 (MQ=255)
                                                                                                                                   tGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAAc               <  1:3339216/150‑1 (MQ=255)
                                                                                                                                   tGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAAc               >  2:1446401/1‑150 (MQ=255)
                                                                                                                                      ttttGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTg            <  2:4145453/150‑1 (MQ=255)
                                                                                                                                       tttGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTgg           <  2:340422/150‑1 (MQ=255)
                                                                                                                                       tttGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTgg           >  1:2211958/1‑150 (MQ=255)
                                                                                                                                           agaaaagagaaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAag       <  2:1274190/150‑1 (MQ=255)
                                                                                                                                            gaaaagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAaga      >  2:4403091/1‑150 (MQ=255)
                                                                                                                                            gaaaagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAaga      >  2:1533485/1‑150 (MQ=255)
                                                                                                                                             aaaagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAagag     <  2:4574094/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGG  >  CP009273/3009040‑3009333

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: