Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,537,806 A→T 89.9% L171* (TTA→TAA)  narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,537,8060AT89.9% 914.1 / 74.8 298L171* (TTA→TAA) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base A (13/17);  new base T (128/140);  total (141/157)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 299 total aligned reads displayed.

CGGTACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTG  >  CP009273/1537661‑1537954
                                                                                                                                                 |                                                                                                                                                    
cGGTACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATc                                                                                                                                                   <  1:1669833/149‑1 (MQ=255)
 ggTACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATcc                                                                                                                                                  <  1:2679325/149‑1 (MQ=255)
    aCCCAACTCCATGCGGCATTCGCCAGCGACATCCCCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGcc                                                                                                                                               <  2:187113/149‑1 (MQ=255)
    aCCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGcc                                                                                                                                               >  1:37239/1‑149 (MQ=255)
      ccAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCAt                                                                                                                                             <  1:501206/149‑1 (MQ=255)
         aaTCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTaa                                                                                                                                          <  1:4400534/149‑1 (MQ=255)
          aTCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAAt                                                                                                                                         >  2:368234/1‑149 (MQ=255)
          aTCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAAt                                                                                                                                         >  2:3168482/1‑149 (MQ=255)
           tCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATc                                                                                                                                        <  1:680193/149‑1 (MQ=255)
              aTGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCACCCAGCTGCATTACACGTACACCTAAGTTTCCTTATCCGCCAGTAATCCCa                                                                                                                                     >  2:4789970/1‑149 (MQ=255)
                     aTTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAgcgc                                                                                                                              <  2:4285149/149‑1 (MQ=255)
                      ttCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCt                                                                                                                             <  2:3895143/149‑1 (MQ=255)
                        cGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTcc                                                                                                                           <  1:981331/149‑1 (MQ=255)
                         gCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTccc                                                                                                                          >  1:1749945/1‑149 (MQ=255)
                          ccAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCt                                                                                                                         <  1:2322001/149‑1 (MQ=255)
                           cAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCtt                                                                                                                        >  2:2557681/1‑149 (MQ=255)
                             gCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGt                                                                                                                      <  2:593698/149‑1 (MQ=255)
                              cGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGtt                                                                                                                     <  1:2213414/149‑1 (MQ=255)
                              cGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGtt                                                                                                                     <  1:1434961/149‑1 (MQ=255)
                               gACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGttt                                                                                                                    <  1:806732/149‑1 (MQ=255)
                                aCATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTg                                                                                                                   >  1:3329845/1‑149 (MQ=255)
                                  aTCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGc                                                                                                                 <  1:2145480/149‑1 (MQ=255)
                                    cACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCtt                                                                                                               <  1:410083/149‑1 (MQ=255)
                                     aCCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCttt                                                                                                              >  1:2000505/1‑149 (MQ=255)
                                     aCCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCttt                                                                                                              >  2:230874/1‑149 (MQ=255)
                                        gAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCACGAGCGCTCCCTTGTTTGGCTTTTgg                                                                                                           >  2:4131955/1‑149 (MQ=255)
                                         aaCCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGa                                                                                                          >  1:1888912/1‑149 (MQ=255)
                                         aaCCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTGTTGGa                                                                                                          >  1:447863/1‑149 (MQ=255)
                                         aaCCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGa                                                                                                          >  1:4504128/1‑149 (MQ=255)
                                           ccGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGaaa                                                                                                        <  1:1610547/149‑1 (MQ=255)
                                            cGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAaag                                                                                                       <  1:1586246/149‑1 (MQ=255)
                                             gTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAaaga                                                                                                      >  2:4226219/1‑149 (MQ=255)
                                              tCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAaagaa                                                                                                     <  1:399751/149‑1 (MQ=255)
                                                 gCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAgaagaa                                                                                                  <  2:4145494/149‑1 (MQ=255)
                                                 gCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAgaagaa                                                                                                  <  2:2462815/149‑1 (MQ=255)
                                                  ccTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGaaa                                                                                                 >  1:3497244/1‑149 (MQ=255)
                                                      cGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGa                                                                                             <  2:798532/149‑1 (MQ=255)
                                                       ggTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGat                                                                                            >  2:695647/1‑149 (MQ=255)
                                                       ggTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGat                                                                                            >  1:1538334/1‑149 (MQ=255)
                                                        gTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGata                                                                                           <  1:2708552/149‑1 (MQ=255)
                                                        gTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGata                                                                                           <  2:4011810/149‑1 (MQ=255)
                                                          aCGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATAtt                                                                                         <  1:3296972/149‑1 (MQ=255)
                                                          aCGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATAtt                                                                                         <  1:2556568/149‑1 (MQ=255)
                                                               aTTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCa                                                                                    <  2:2093831/149‑1 (MQ=255)
                                                                ttGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCAt                                                                                   >  2:3962170/1‑149 (MQ=255)
                                                                ttGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCAt                                                                                   >  1:2280040/1‑149 (MQ=255)
                                                                ttGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCAt                                                                                   <  1:1928262/149‑1 (MQ=255)
                                                                 tGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATg                                                                                  >  1:282034/1‑149 (MQ=255)
                                                                  gACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGc                                                                                 >  1:516510/1‑149 (MQ=255)
                                                                      ccGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGa                                                                             >  2:1096235/1‑149 (MQ=255)
                                                                      ccGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGa                                                                             >  2:2769863/1‑149 (MQ=255)
                                                                            aaGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCaaa                                                                       >  1:3901548/1‑149 (MQ=255)
                                                                                 aaaTACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTg                                                                  >  2:4682756/1‑149 (MQ=255)
                                                                                   aTACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCa                                                                <  2:4757574/149‑1 (MQ=255)
                                                                                    tACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCAc                                                               >  2:1732763/1‑149 (MQ=255)
                                                                                     aCAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCAcc                                                              >  2:2484440/1‑149 (MQ=255)
                                                                                     aCAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCAcc                                                              <  1:842000/149‑1 (MQ=255)
                                                                                        ggTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGc                                                           >  1:2436599/1‑149 (MQ=255)
                                                                                        ggTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGc                                                           >  2:1207338/1‑149 (MQ=255)
                                                                                           aCAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCaaa                                                        >  2:1392361/1‑149 (MQ=255)
                                                                                           aCAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCaaa                                                        >  1:4537342/1‑149 (MQ=255)
                                                                                             aaaaaTGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGACGCAAAGTTTGCACCTGCAAAAc                                                      >  2:4293804/1‑149 (MQ=255)
                                                                                             aaaaaTGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAAc                                                      <  2:3008866/149‑1 (MQ=255)
                                                                                             aaaaaTGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAAc                                                      <  2:3223568/149‑1 (MQ=255)
                                                                                               aaaTGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCg                                                    >  2:4393372/1‑149 (MQ=255)
                                                                                               aaaTGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTGTGCACCTGCAAAACCg                                                    >  2:2002559/1‑149 (MQ=255)
                                                                                                aaTGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGc                                                   <  2:256187/149‑1 (MQ=255)
                                                                                                 aTGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCa                                                  >  2:3369755/1‑149 (MQ=255)
                                                                                                  tGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCAt                                                 >  1:2093972/1‑149 (MQ=255)
                                                                                                   gACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATa                                                <  1:3159819/149‑1 (MQ=255)
                                                                                                   gACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATa                                                <  1:3656735/149‑1 (MQ=255)
                                                                                                       aGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCaa                                            >  2:656017/1‑149 (MQ=255)
                                                                                                        gCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCaaa                                           <  1:4373904/149‑1 (MQ=255)
                                                                                                            tGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcga                                       <  2:3605062/149‑1 (MQ=255)
                                                                                                              cAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgata                                     <  1:1265434/149‑1 (MQ=255)
                                                                                                                aCCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataac                                   >  2:4447531/1‑149 (MQ=255)
                                                                                                                aCCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataac                                   >  1:2502982/1‑149 (MQ=255)
                                                                                                                aCCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataac                                   >  1:2856002/1‑149 (MQ=255)
                                                                                                                  cAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataacga                                 <  1:1197489/149‑1 (MQ=255)
                                                                                                                   aGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataacgat                                >  2:647826/1‑149 (MQ=255)
                                                                                                                   aGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataacgat                                <  1:2978256/149‑1 (MQ=255)
                                                                                                                    gCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataacgata                               >  1:2454685/1‑149 (MQ=255)
                                                                                                                          ttacacttacacCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATc                         >  2:341129/1‑149 (MQ=255)
                                                                                                                          ttacacttacacCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATc                         >  1:3978364/1‑149 (MQ=255)
                                                                                                                           tacacttacacCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATcc                        <  1:3622392/149‑1 (MQ=255)
                                                                                                                                 tacacCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAg                  <  1:1204403/149‑1 (MQ=255)
                                                                                                                                  acacCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAgg                 <  1:4689544/149‑1 (MQ=255)
                                                                                                                                      cTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTa             <  1:953174/149‑1 (MQ=255)
                                                                                                                                       tAAGTTTCCTTATCCTCCATTAATCCCAAGAGCGATCCATTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTTCAAAACCTCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTAt            <  1:1720049/149‑1 (MQ=255)
                                                                                                                                       tAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTAt            >  1:4632364/1‑149 (MQ=255)
                                                                                                                                       tAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTAt            <  1:1033630/149‑1 (MQ=255)
                                                                                                                                        aaGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTAtt           <  1:2616196/149‑1 (MQ=255)
                                                                                                                                         aGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTc          >  2:1578894/1‑149 (MQ=255)
                                                                                                                                         aGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTc          <  1:4230331/149‑1 (MQ=255)
                                                                                                                                         aGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTc          <  1:3025161/149‑1 (MQ=255)
                                                                                                                                         aGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTc          <  2:976645/149‑1 (MQ=255)
                                                                                                                                          gTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCg         <  2:2528246/149‑1 (MQ=255)
                                                                                                                                            ttCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGa       >  1:1338726/1‑149 (MQ=255)
                                                                                                                                             tCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGAt      >  2:3823993/1‑149 (MQ=255)
                                                                                                                                             tCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGAt      <  1:3368408/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGGTACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTG  >  CP009273/1537661‑1537954

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: