Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,009,187 G→T 49.8% 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,1870GT49.8% 57.6 / 307.9 201intergenic (+90/‑232)mocA/ygfKCTP:molybdopterin cytidylyltransferase/putative Fe‑S subunit oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base G (61/40);  new base T (47/53);  total (108/93)
Fisher's exact test for biased strand distribution p-value = 6.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 223 total aligned reads displayed.

AGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTG  >  CP009273/3009040‑3009335
                                                                                                                                                   |                                                                                                                                                    
aGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaaaagataa                                                                                                                                                    >  2:519659/1‑150 (MQ=255)
aGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaaaagagaa                                                                                                                                                    <  1:2809630/150‑1 (MQ=255)
   tatTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGaaaaa                                                                                                                                                 >  2:84648/1‑150 (MQ=255)
    atTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAAt                                                                                                                                                <  2:2827054/150‑1 (MQ=255)
       gATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGa                                                                                                                                             <  1:1250041/150‑1 (MQ=255)
        aTACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGAt                                                                                                                                            >  2:1636885/1‑150 (MQ=255)
        aTACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGAt                                                                                                                                            >  1:2883256/1‑150 (MQ=255)
         tACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATa                                                                                                                                           >  1:1045842/1‑150 (MQ=255)
         tACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATa                                                                                                                                           >  1:921388/1‑150 (MQ=255)
          aCTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATaa                                                                                                                                          <  2:3967362/150‑1 (MQ=255)
            tCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAAt                                                                                                                                        >  1:4022858/1‑150 (MQ=255)
              cGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATgg                                                                                                                                      >  1:1717718/1‑150 (MQ=255)
               ggACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCATCCGGCACTCATTATCAAAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATgtg                                                                                                                                     >  1:837944/1‑148 (MQ=255)
                  cgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTg                                                                                                                                  >  2:4133615/1‑150 (MQ=255)
                  cgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTg                                                                                                                                  <  1:1927399/150‑1 (MQ=255)
                   gacTTTATTACAGCTAAGGAATGTTATACTTATATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTgg                                                                                                                                 <  1:2550234/150‑1 (MQ=255)
                   gacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTgg                                                                                                                                 >  1:749302/1‑150 (MQ=255)
                        tattaCAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCaata                                                                                                                            >  2:25350/1‑150 (MQ=255)
                           taCAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCaataata                                                                                                                         >  1:2133367/1‑150 (MQ=255)
                                cGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAAtttt                                                                                                                    >  1:260225/1‑150 (MQ=255)
                                cGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAAtttt                                                                                                                    >  1:4014752/1‑150 (MQ=255)
                                 gAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTc                                                                                                                   >  1:4079261/1‑150 (MQ=255)
                                 gAACGAAAGGTATACTGAAATATAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTc                                                                                                                   <  2:3290481/150‑1 (MQ=255)
                                  aaGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCa                                                                                                                  >  1:4275462/1‑150 (MQ=255)
                                    ggAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAgg                                                                                                                >  2:2903598/1‑150 (MQ=255)
                                    ggAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAgg                                                                                                                <  1:1591550/150‑1 (MQ=255)
                                     gAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGa                                                                                                               <  1:698507/150‑1 (MQ=255)
                                     gAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGa                                                                                                               <  1:1828580/150‑1 (MQ=255)
                                      aaaGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAATTGGGCTGGCAATAATAATTTTCAGGAg                                                                                                              >  2:2676924/1‑150 (MQ=255)
                                      aaaGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAg                                                                                                              <  2:189645/150‑1 (MQ=255)
                                       aaGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGt                                                                                                             >  1:3532321/1‑150 (MQ=255)
                                          gTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAt                                                                                                          <  2:2625153/150‑1 (MQ=255)
                                          gTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAt                                                                                                          <  1:1839875/150‑1 (MQ=255)
                                           tataCTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAtt                                                                                                         <  2:231806/150‑1 (MQ=255)
                                           tataCTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAtt                                                                                                         <  2:2407847/150‑1 (MQ=255)
                                            ataCTGAAATGTAGAAAACGTAGTTAAACGATTGAGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAttt                                                                                                        <  2:2977247/150‑1 (MQ=255)
                                              aCTTAAATTTAGAAAACGAAGTTAAACCATTGATTAAAAATATTTAATCCTTCCTGCACTCATTATCATAATGATAAGGAATCAATGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTa                                                                                                      <  2:3762142/150‑1 (MQ=255)
                                                   aaTTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATtct                                                                                                 <  1:1538093/150‑1 (MQ=255)
                                                     tttAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCa                                                                                               >  1:1975354/1‑150 (MQ=255)
                                                        aGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCAttt                                                                                            <  1:3235677/150‑1 (MQ=255)
                                                        aGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCAttt                                                                                            >  2:3949851/1‑150 (MQ=255)
                                                        aGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCAttt                                                                                            <  1:3120179/150‑1 (MQ=255)
                                                         gAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCCTTTa                                                                                           >  2:439695/1‑150 (MQ=255)
                                                            aaCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCa                                                                                        <  2:2424328/150‑1 (MQ=255)
                                                            aaCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCa                                                                                        <  1:20753/150‑1 (MQ=255)
                                                             aCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCAc                                                                                       <  2:4275462/150‑1 (MQ=255)
                                                               gTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCt                                                                                     <  1:3032181/150‑1 (MQ=255)
                                                                tAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCtg                                                                                    <  2:846043/150‑1 (MQ=255)
                                                                     aaaCGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGt                                                                               <  1:2073903/150‑1 (MQ=255)
                                                                           tttcGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttcattc                                                                         <  2:372822/147‑1 (MQ=255)
                                                                           ttGCGTTCAAATAATTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttcattc                                                                         <  1:1291893/150‑1 (MQ=255)
                                                                            tGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttcattca                                                                        >  1:2071569/1‑150 (MQ=255)
                                                                                ttCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACaaa                                                                    >  1:3186186/1‑150 (MQ=255)
                                                                                  cAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACaaaaa                                                                  >  2:3046066/1‑150 (MQ=255)
                                                                                   aaaTATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAg                                                                 <  1:1049717/150‑1 (MQ=255)
                                                                                     atatTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATTGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGcc                                                               >  1:2879921/1‑150 (MQ=255)
                                                                                             tCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGtt                                                       >  2:3180883/1‑150 (MQ=255)
                                                                                              ccTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGttt                                                      <  1:2596339/150‑1 (MQ=255)
                                                                                              ccTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGttt                                                      >  2:3103506/1‑150 (MQ=255)
                                                                                               cTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTc                                                     >  1:1016678/1‑150 (MQ=255)
                                                                                                ttCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCg                                                    >  1:2753254/1‑150 (MQ=255)
                                                                                                 tCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGa                                                   >  1:3956983/1‑150 (MQ=255)
                                                                                                   cGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAgc                                                 <  2:1962549/150‑1 (MQ=255)
                                                                                                    ggCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAgcc                                                >  2:404117/1‑150 (MQ=255)
                                                                                                        cTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCt                                            <  2:114216/150‑1 (MQ=255)
                                                                                                         tCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTc                                           >  1:242246/1‑150 (MQ=255)
                                                                                                           aTTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTcac                                         >  1:492436/1‑150 (MQ=255)
                                                                                                            ttATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTcaca                                        <  1:1254880/150‑1 (MQ=255)
                                                                                                            ttATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTcaca                                        >  2:1438180/1‑150 (MQ=255)
                                                                                                             tATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAt                                       <  1:10482/150‑1 (MQ=255)
                                                                                                              aTCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAtt                                      <  2:4237079/150‑1 (MQ=255)
                                                                                                              aTCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGAAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAtt                                      <  1:3754076/150‑1 (MQ=255)
                                                                                                                cATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAtttt                                    >  2:1082584/1‑150 (MQ=255)
                                                                                                                cATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAtttt                                    <  2:3600289/150‑1 (MQ=255)
                                                                                                                     tGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCACTAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTtat                               >  2:4220090/1‑150 (MQ=255)
                                                                                                                     tGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTtat                               <  1:3362180/150‑1 (MQ=255)
                                                                                                                     tGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTtat                               >  1:2105701/1‑150 (MQ=255)
                                                                                                                      gATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTtata                              <  1:3524016/150‑1 (MQ=255)
                                                                                                                       aTAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTatat                             >  1:3388013/1‑150 (MQ=255)
                                                                                                                        tAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATAtt                            >  2:3623728/1‑150 (MQ=255)
                                                                                                                        tAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATAtt                            >  2:3396757/1‑150 (MQ=255)
                                                                                                                           ggAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTcc                         <  2:749302/150‑1 (MQ=255)
                                                                                                                            gAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTccc                        <  1:2239714/150‑1 (MQ=255)
                                                                                                                               tCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTgccgcc                     >  1:3606507/1‑150 (MQ=255)
                                                                                                                                 aCTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCaa                   >  1:748863/1‑150 (MQ=255)
                                                                                                                                  cTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCaaa                  <  1:519994/150‑1 (MQ=255)
                                                                                                                                     aTTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCt               >  2:1003927/1‑150 (MQ=255)
                                                                                                                                        ttGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGc            <  2:1045842/150‑1 (MQ=255)
                                                                                                                                        ttGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGc            <  1:1041720/150‑1 (MQ=255)
                                                                                                                                        ttGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGc            <  1:1108587/150‑1 (MQ=255)
                                                                                                                                         tGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCa           >  2:2895727/1‑150 (MQ=255)
                                                                                                                                          gagaaaagagaaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCACGCCAAACCTGGCaa          >  2:3090879/1‑150 (MQ=255)
                                                                                                                                            gaaaagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAaga        >  2:392996/1‑150 (MQ=255)
                                                                                                                                            gaaaagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAaga        >  2:1688486/1‑150 (MQ=255)
                                                                                                                                            gaaaagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAaga        >  1:3405903/1‑150 (MQ=255)
                                                                                                                                             aaaagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAagag       >  1:878810/1‑150 (MQ=255)
                                                                                                                                                agataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAgtgg    <  1:2472268/150‑1 (MQ=255)
                                                                                                                                                agagaaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAgtgg    <  2:2288038/150‑1 (MQ=255)
                                                                                                                                                agagaaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAgtgg    >  1:3735395/1‑150 (MQ=255)
                                                                                                                                                  agaaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAgtggtg  >  2:1698801/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTG  >  CP009273/3009040‑3009335

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: