Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,489,397 C→A 100% F676L (TTC→TTA nirB → nitrite reductase, large subunit, NAD(P)H‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,489,3970CA99.3% 988.6 / ‑1.9 292F676L (TTC→TTAnirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (137/152);  total (138/154)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 292 total aligned reads displayed.

GGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATATCG  >  CP009273/3489251‑3489542
                                                                                                                                                  |                                                                                                                                                 
ggTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAgg                                                                                                                                                 <  2:2561680/149‑1 (MQ=255)
 gTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGgt                                                                                                                                                >  1:3077996/1‑149 (MQ=255)
 gTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGgt                                                                                                                                                >  1:3321147/1‑149 (MQ=255)
 gTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGgt                                                                                                                                                >  1:1267617/1‑149 (MQ=255)
 gTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGgt                                                                                                                                                <  2:1179923/149‑1 (MQ=255)
   cATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGgtgt                                                                                                                                              >  1:3442243/1‑149 (MQ=255)
   cATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGgtgt                                                                                                                                              >  1:2561246/1‑149 (MQ=255)
    aTGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGtc                                                                                                                                             >  1:2707517/1‑149 (MQ=255)
        cTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCg                                                                                                                                         >  2:2957036/1‑149 (MQ=255)
          aTGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGc                                                                                                                                       >  1:1012004/1‑149 (MQ=255)
          aTGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGc                                                                                                                                       <  1:2615852/149‑1 (MQ=255)
                aaGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTAcc                                                                                                                                 >  1:1417657/1‑149 (MQ=255)
                 aGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTAccc                                                                                                                                <  1:3721558/149‑1 (MQ=255)
                  gCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGACCTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCg                                                                                                                               >  1:3190720/1‑149 (MQ=255)
                    aCTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAACGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTg                                                                                                                             <  1:3560533/149‑1 (MQ=255)
                    aCTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTg                                                                                                                             <  1:1595960/149‑1 (MQ=255)
                    aCTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTg                                                                                                                             >  2:826261/1‑149 (MQ=255)
                     cTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGa                                                                                                                            <  1:155319/149‑1 (MQ=255)
                     cTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGa                                                                                                                            >  2:2746899/1‑149 (MQ=255)
                      tGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGaa                                                                                                                           <  2:4055080/149‑1 (MQ=255)
                      tGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGaa                                                                                                                           >  2:4703619/1‑149 (MQ=255)
                       gCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAAt                                                                                                                          <  2:3952888/149‑1 (MQ=255)
                         gTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGt                                                                                                                        <  1:2390204/149‑1 (MQ=255)
                            tGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCa                                                                                                                     <  1:2743307/149‑1 (MQ=255)
                             ggCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAg                                                                                                                    >  1:765443/1‑149 (MQ=255)
                              gCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGa                                                                                                                   <  2:3720997/149‑1 (MQ=255)
                               cGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGaa                                                                                                                  <  1:3398375/149‑1 (MQ=255)
                               cGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGaa                                                                                                                  <  2:4347832/149‑1 (MQ=255)
                                 aaaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGc                                                                                                                >  1:3828710/1‑149 (MQ=255)
                                  aaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCt                                                                                                               <  1:1759092/149‑1 (MQ=255)
                                  aaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCt                                                                                                               <  2:1408132/149‑1 (MQ=255)
                                  aaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCt                                                                                                               <  2:3923214/149‑1 (MQ=255)
                                  aaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCt                                                                                                               <  1:4148713/149‑1 (MQ=255)
                                    aCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCa                                                                                                             >  2:847873/1‑149 (MQ=255)
                                      cTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAgg                                                                                                           <  2:2081579/149‑1 (MQ=255)
                                       tGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAggg                                                                                                          <  2:194997/149‑1 (MQ=255)
                                       tGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAggg                                                                                                          >  2:4288550/1‑149 (MQ=255)
                                             ggTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGa                                                                                                    <  1:2719726/149‑1 (MQ=255)
                                             ggTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGa                                                                                                    <  2:3607645/149‑1 (MQ=255)
                                              gTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGAc                                                                                                   <  1:727818/149‑1 (MQ=255)
                                              gTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGAc                                                                                                   <  2:2806448/149‑1 (MQ=255)
                                                 gCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTg                                                                                                >  1:4515663/1‑149 (MQ=255)
                                                 gCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTg                                                                                                <  1:1051187/149‑1 (MQ=255)
                                                  cACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTgg                                                                                               <  2:4253536/149‑1 (MQ=255)
                                                   aCCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTggg                                                                                              <  2:2242130/149‑1 (MQ=255)
                                                   aCCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTggg                                                                                              <  2:3174121/149‑1 (MQ=255)
                                                    ccTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGt                                                                                             <  2:2176254/149‑1 (MQ=255)
                                                    ccTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGt                                                                                             <  2:2313666/149‑1 (MQ=255)
                                                    ccTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGt                                                                                             <  1:437780/149‑1 (MQ=255)
                                                       ggTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtatt                                                                                          >  1:2200215/1‑149 (MQ=255)
                                                               tACGGCGTTGGCGACAGCGTCGGTCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCAc                                                                                  >  2:1375318/1‑149 (MQ=255)
                                                               tACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAATTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCAc                                                                                  <  1:4436981/149‑1 (MQ=255)
                                                               tACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCAc                                                                                  <  2:1952555/149‑1 (MQ=255)
                                                               tACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACCGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCAc                                                                                  <  2:333246/149‑1 (MQ=255)
                                                                   gCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGaa                                                                              >  2:3764979/1‑149 (MQ=255)
                                                                    cGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGaaa                                                                             <  2:1771583/149‑1 (MQ=255)
                                                                           gACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCACACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTg                                                                      >  2:2663216/1‑149 (MQ=255)
                                                                            aCAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTgg                                                                     <  1:2878678/149‑1 (MQ=255)
                                                                             cAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGa                                                                    <  1:1772975/149‑1 (MQ=255)
                                                                              aGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGACAAAGGCTGGaa                                                                   >  1:1145313/1‑149 (MQ=255)
                                                                                cGTCGGCCTCGGCGTGGAACTGGAAAACCGCGACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAcc                                                                 <  2:3814755/149‑1 (MQ=255)
                                                                                 gTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCt                                                                <  1:3230971/149‑1 (MQ=255)
                                                                                  tCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTg                                                               <  2:158519/149‑1 (MQ=255)
                                                                                  tCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTg                                                               <  2:3006762/149‑1 (MQ=255)
                                                                                      ccTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATg                                                           <  1:3513893/149‑1 (MQ=255)
                                                                                        tCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGtt                                                         <  2:4473997/149‑1 (MQ=255)
                                                                                          ggCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTg                                                       <  1:1987266/149‑1 (MQ=255)
                                                                                           gCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGc                                                      >  1:848546/1‑149 (MQ=255)
                                                                                             gTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCgg                                                    <  2:4203839/149‑1 (MQ=255)
                                                                                             gTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCgg                                                    <  2:4405382/149‑1 (MQ=255)
                                                                                                gAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTaa                                                 <  2:563130/149‑1 (MQ=255)
                                                                                                  aCTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAAcg                                               >  2:1656599/1‑149 (MQ=255)
                                                                                                   cTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAAcgg                                              <  1:1445087/149‑1 (MQ=255)
                                                                                                        aaaCCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCa                                         >  1:588491/1‑149 (MQ=255)
                                                                                                         aaCCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGTCGTGGGTATTATCGCCACTGAAAAAGGCTGGAAACTGTATGTTTGCGGTAACGGCGGCAt                                        <  2:1095450/149‑1 (MQ=255)
                                                                                                         aaCCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCAt                                        >  1:4286116/1‑149 (MQ=255)
                                                                                                         aaCCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCAt                                        <  2:721475/149‑1 (MQ=255)
                                                                                                          aCCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATg                                       <  2:921466/149‑1 (MQ=255)
                                                                                                           ccGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGa                                      >  2:1311881/1‑149 (MQ=255)
                                                                                                             gCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGaaa                                    <  2:1012004/149‑1 (MQ=255)
                                                                                                                 cAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACcgc                                <  1:4565363/149‑1 (MQ=255)
                                                                                                                 cAAAGGCATCCGTAAGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCACAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGACCCTCTATGTATGCGGTAACGGCGGCAGGAAACcgc                                >  2:2872559/1‑149 (MQ=255)
                                                                                                                   aaGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGt                              <  2:3568917/149‑1 (MQ=255)
                                                                                                                   aaGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGt                              <  2:3007589/149‑1 (MQ=255)
                                                                                                                      gCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCAt                           >  1:1035683/1‑149 (MQ=255)
                                                                                                                       cATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATg                          <  2:569315/149‑1 (MQ=255)
                                                                                                                         tCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAAcc                                  <  2:3402934/139‑1 (MQ=255)
                                                                                                                          ccGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCgg                       <  1:4054109/149‑1 (MQ=255)
                                                                                                                          ccGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCgg                       >  1:316941/1‑149 (MQ=255)
                                                                                                                           cGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGa                      <  1:702233/149‑1 (MQ=255)
                                                                                                                             tACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATc                    <  1:3486427/149‑1 (MQ=255)
                                                                                                                                gccgcACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATctgc                 <  2:4425529/149‑1 (MQ=255)
                                                                                                                                     acaAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggcg            >  2:617529/1‑149 (MQ=255)
                                                                                                                                      caAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCGGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggcgg           >  2:2358511/1‑149 (MQ=255)
                                                                                                                                        aaaTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGgcggcg         >  2:4186184/1‑149 (MQ=255)
                                                                                                                                          aTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGa       >  2:2591794/1‑149 (MQ=255)
                                                                                                                                           tGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGat      >  2:3199855/1‑149 (MQ=255)
                                                                                                                                            gAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCGTGCGGATCTGCTGGCGGCGGata     <  1:995905/149‑1 (MQ=255)
                                                                                                                                             aaGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCATTAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGatat    <  1:4193407/149‑1 (MQ=255)
                                                                                                                                               gTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATatcg  >  1:4748263/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATATCG  >  CP009273/3489251‑3489542

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: