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.4% 608.5 / ‑4.2 188F676L (TTC→TTAnirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (87/99);  total (88/100)
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 188 total aligned reads displayed.

GTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATAT  >  CP009273/3489252‑3489540
                                                                                                                                                 |                                                                                                                                               
gTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAgg                                                                                                                                               <  2:2293843/148‑1 (MQ=255)
 tCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGTACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGgt                                                                                                                                              >  2:3544794/1‑148 (MQ=255)
   aTGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGgtgt                                                                                                                                            <  1:3052564/148‑1 (MQ=255)
    tGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGtc                                                                                                                                           >  2:1084418/1‑148 (MQ=255)
      ccTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGtctc                                                                                                                                         <  2:1871993/148‑1 (MQ=255)
      ccTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGAACGCCGCACACAATGAAGTTAGGTGtctc                                                                                                                                         >  1:2249030/1‑148 (MQ=255)
          tGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGc                                                                                                                                     <  2:1069746/148‑1 (MQ=255)
          tGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGc                                                                                                                                     <  2:720680/148‑1 (MQ=255)
           gCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCt                                                                                                                                    <  1:1147200/148‑1 (MQ=255)
             gAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGt                                                                                                                                  <  2:296273/148‑1 (MQ=255)
                 gCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTAccc                                                                                                                              >  2:2307321/1‑148 (MQ=255)
                 gCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTAccc                                                                                                                              >  2:3424298/1‑148 (MQ=255)
                   aCTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGt                                                                                                                            <  1:2824719/148‑1 (MQ=255)
                    cTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTg                                                                                                                           >  1:2630933/1‑148 (MQ=255)
                     tGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGa                                                                                                                          <  2:3038882/148‑1 (MQ=255)
                      gCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGaa                                                                                                                         >  2:2997755/1‑148 (MQ=255)
                      gCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGaa                                                                                                                         >  2:2997789/1‑148 (MQ=255)
                       cGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAAt                                                                                                                        <  1:2515388/148‑1 (MQ=255)
                         tATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGt                                                                                                                      <  2:3182648/148‑1 (MQ=255)
                             gCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAg                                                                                                                  <  2:1430811/148‑1 (MQ=255)
                                aaaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAg                                                                                                               >  1:2777617/1‑148 (MQ=255)
                                 aaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGc                                                                                                              <  1:2476122/148‑1 (MQ=255)
                                 aaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACAGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGc                                                                                                              <  2:794776/148‑1 (MQ=255)
                                  aaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCt                                                                                                             >  1:163289/1‑148 (MQ=255)
                                  aaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCt                                                                                                             <  1:2970007/148‑1 (MQ=255)
                                  aaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCt                                                                                                             >  1:419491/1‑148 (MQ=255)
                                     cTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAg                                                                                                          <  2:2668937/148‑1 (MQ=255)
                                      tGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAgg                                                                                                         <  1:3265511/148‑1 (MQ=255)
                                       gCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAggg                                                                                                        >  1:338834/1‑148 (MQ=255)
                                       gCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAggg                                                                                                        >  1:339256/1‑148 (MQ=255)
                                         gTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTa                                                                                                      <  2:2603622/148‑1 (MQ=255)
                                          tGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTaa                                                                                                     >  2:1804363/1‑148 (MQ=255)
                                          tGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTaa                                                                                                     <  1:1493714/148‑1 (MQ=255)
                                           gggTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTaaa                                                                                                    >  2:2694645/1‑148 (MQ=255)
                                            ggTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAg                                                                                                   <  1:3111283/148‑1 (MQ=255)
                                              tAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGACCTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGAc                                                                                                 >  2:2064801/1‑148 (MQ=255)
                                               aGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACg                                                                                                >  2:2466838/1‑148 (MQ=255)
                                               aGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACg                                                                                                >  2:197227/1‑148 (MQ=255)
                                                gCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGt                                                                                               >  1:1197552/1‑148 (MQ=255)
                                                 cACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACCAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTg                                                                                              >  2:1266683/1‑148 (MQ=255)
                                                 cACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTg                                                                                              >  1:1635546/1‑148 (MQ=255)
                                                      ggTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtat                                                                                         >  2:1500139/1‑148 (MQ=255)
                                                      ggTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtat                                                                                         <  1:939066/148‑1 (MQ=255)
                                                       gTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtatt                                                                                        >  2:1281505/1‑148 (MQ=255)
                                                             cTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGcc                                                                                  >  1:1991123/1‑148 (MQ=255)
                                                              tACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCa                                                                                 <  2:1673290/148‑1 (MQ=255)
                                                               aCGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCAc                                                                                >  1:796876/1‑148 (MQ=255)
                                                               aCGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCAc                                                                                >  1:453237/1‑148 (MQ=255)
                                                                cGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACt                                                                               >  1:1719726/1‑148 (MQ=255)
                                                                 ggCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTg                                                                              >  1:869656/1‑148 (MQ=255)
                                                                     ttGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGaaaa                                                                          >  2:2515379/1‑148 (MQ=255)
                                                                     ttGGCGACAGCGTCGGCCTCGGCGTGGAACTCGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGaaaa                                                                          <  2:1931798/148‑1 (MQ=255)
                                                                        gCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAgg                                                                       <  1:1726978/148‑1 (MQ=255)
                                                                         cGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGc                                                                      <  2:2741104/148‑1 (MQ=255)
                                                                          gACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCt                                                                     <  2:2881550/148‑1 (MQ=255)
                                                                           aCAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTg                                                                    <  2:926423/148‑1 (MQ=255)
                                                                             aGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGa                                                                  >  1:2149043/1‑148 (MQ=255)
                                                                                gTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAcc                                                               <  1:1829672/148‑1 (MQ=255)
                                                                                gTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAcc                                                               <  1:2840770/148‑1 (MQ=255)
                                                                                 tCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCt                                                              <  2:2795829/148‑1 (MQ=255)
                                                                                  cGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTg                                                             <  2:217265/148‑1 (MQ=255)
                                                                                   ggCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGt                                                            >  2:113755/1‑148 (MQ=255)
                                                                                    gCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTa                                                           <  2:1794180/148‑1 (MQ=255)
                                                                                    gCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTa                                                           >  2:1752208/1‑148 (MQ=255)
                                                                                    gCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTa                                                           <  1:1663914/148‑1 (MQ=255)
                                                                                      cTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATg                                                         <  2:3083698/148‑1 (MQ=255)
                                                                                       tCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGt                                                        <  2:477673/148‑1 (MQ=255)
                                                                                         ggagTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGttt                                                      <  1:2567821/145‑1 (MQ=255)
                                                                                             tGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCgg                                                  <  2:2236954/148‑1 (MQ=255)
                                                                                             tGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCgg                                                  <  2:190862/148‑1 (MQ=255)
                                                                                             tGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCgg                                                  >  2:497334/1‑148 (MQ=255)
                                                                                             tGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCgg                                                  <  1:2405100/148‑1 (MQ=255)
                                                                                              ggAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGt                                                 <  2:2969211/148‑1 (MQ=255)
                                                                                                     gAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAAcggcg                                          >  2:1429677/1‑148 (MQ=255)
                                                                                                     gAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAAcggcg                                          >  1:3320376/1‑148 (MQ=255)
                                                                                                       aaaCCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACggcggc                                        >  2:1415456/1‑148 (MQ=255)
                                                                                                           cGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGa                                    >  1:917736/1‑148 (MQ=255)
                                                                                                               aCAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAAcc                                >  2:3140317/1‑148 (MQ=255)
                                                                                                                cAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACcg                               >  2:1317545/1‑148 (MQ=255)
                                                                                                                  aaGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACcgcg                             <  2:1575980/148‑1 (MQ=255)
                                                                                                                  aaGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACcgcg                             <  2:579512/148‑1 (MQ=255)
                                                                                                                  aaGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACcgcg                             >  1:393933/1‑148 (MQ=255)
                                                                                                                  aaGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTAGGGTTGCGGTAACGGCGGCATGAAACcgcg                             >  2:2631730/1‑148 (MQ=255)
                                                                                                                   aGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGt                            >  1:1475736/1‑148 (MQ=255)
                                                                                                                      cATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCAt                         <  1:2936179/148‑1 (MQ=255)
                                                                                                                        tCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGc                       >  1:3470027/1‑148 (MQ=255)
                                                                                                                           gTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGa                    <  1:100544/148‑1 (MQ=255)
                                                                                                                            tACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGAt                   <  2:962632/148‑1 (MQ=255)
                                                                                                                            tACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGAt                   <  2:1991136/148‑1 (MQ=255)
                                                                                                                                ccgcACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATctgc               >  2:3354497/1‑148 (MQ=255)
                                                                                                                                   cacaAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTgg            >  1:2069682/1‑148 (MQ=255)
                                                                                                                                     caAACTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggcg          <  2:2287705/148‑1 (MQ=255)
                                                                                                                                      aaaaTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggcgg         >  2:1715417/1‑148 (MQ=255)
                                                                                                                                      aaaaTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggcgg         >  2:1714882/1‑148 (MQ=255)
                                                                                                                                         aTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGcggcgg      >  2:3368580/1‑148 (MQ=255)
                                                                                                                                         aTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGcggcgg      >  2:3368617/1‑148 (MQ=255)
                                                                                                                                         aTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGcggcgg      >  1:1680083/1‑148 (MQ=255)
                                                                                                                                            aaGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGata   >  1:821336/1‑148 (MQ=255)
                                                                                                                                             aGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGatat  <  1:652290/148‑1 (MQ=255)
                                                                                                                                             aGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGatat  <  1:1031144/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATAT  >  CP009273/3489252‑3489540

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: