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,3970CA98.8% 805.5 / 0.2 251F676L (TTC→TTAnirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base A (115/132);  total (115/136)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01
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 251 total aligned reads displayed.

GGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATATCG  >  CP009273/3489251‑3489542
                                                                                                                                                  |                                                                                                                                                 
 gTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGgt                                                                                                                                                >  1:2457997/1‑149 (MQ=255)
 gTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGgt                                                                                                                                                >  1:3741399/1‑149 (MQ=255)
    aTGCCTATGCGAAAGCACTGCGGATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGtc                                                                                                                                             <  2:1832391/149‑1 (MQ=255)
     tGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGtct                                                                                                                                            >  2:4025944/1‑149 (MQ=255)
       ccTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTcc                                                                                                                                          >  2:2546059/1‑149 (MQ=255)
                aaGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTAcc                                                                                                                                 <  2:3229705/149‑1 (MQ=255)
                 aGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTAccc                                                                                                                                >  1:1279231/1‑149 (MQ=255)
                     cTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGa                                                                                                                            >  1:3433907/1‑149 (MQ=255)
                      tGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGaa                                                                                                                           >  1:87357/1‑149 (MQ=255)
                      tGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGaa                                                                                                                           >  2:3629535/1‑149 (MQ=255)
                       gCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAAt                                                                                                                          >  1:1955316/1‑149 (MQ=255)
                         gTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGt                                                                                                                        >  2:2161358/1‑149 (MQ=255)
                          tATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGtt                                                                                                                       <  2:810637/149‑1 (MQ=255)
                            tGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCa                                                                                                                     >  2:1904723/1‑149 (MQ=255)
                            tGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCa                                                                                                                     <  2:867422/149‑1 (MQ=255)
                            tGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCa                                                                                                                     <  1:3475319/149‑1 (MQ=255)
                             ggCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAg                                                                                                                    >  2:1536984/1‑149 (MQ=255)
                               cGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGaa                                                                                                                  <  1:216615/149‑1 (MQ=255)
                               cGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGaa                                                                                                                  >  2:3571641/1‑149 (MQ=255)
                                gAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGACCAGGAAAACCGCTACAAAGGCAACCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAg                                                                                                                 >  1:2609587/1‑149 (MQ=255)
                                gAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAg                                                                                                                 <  2:290860/149‑1 (MQ=255)
                                 aaaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGc                                                                                                                <  2:880906/149‑1 (MQ=255)
                                  aaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCt                                                                                                               <  1:1129204/149‑1 (MQ=255)
                                   aaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTc                                                                                                              <  1:3276701/149‑1 (MQ=255)
                                     ccTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAg                                                                                                            <  2:3713160/149‑1 (MQ=255)
                                      cTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAgg                                                                                                           >  1:4069348/1‑149 (MQ=255)
                                        gCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGt                                                                                                         >  2:2759880/1‑149 (MQ=255)
                                        gCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGt                                                                                                         <  2:3938492/149‑1 (MQ=255)
                                           tGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTaaa                                                                                                      <  1:1536984/149‑1 (MQ=255)
                                           tGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTaaa                                                                                                      <  2:3985429/149‑1 (MQ=255)
                                           tGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTaaa                                                                                                      <  1:3541873/149‑1 (MQ=255)
                                                aGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGt                                                                                                 <  1:3330546/149‑1 (MQ=255)
                                                 gCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTg                                                                                                <  1:972074/149‑1 (MQ=255)
                                                 gCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTg                                                                                                >  1:768792/1‑149 (MQ=255)
                                                 gCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTg                                                                                                <  2:3080369/149‑1 (MQ=255)
                                                 gCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTg                                                                                                >  1:332033/1‑149 (MQ=255)
                                                  cACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTgg                                                                                               >  1:1501518/1‑149 (MQ=255)
                                                     cTGGTGCCGCTACGGAGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGta                                                                                            <  1:1830152/149‑1 (MQ=255)
                                                       ggTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtatt                                                                                          >  1:1732019/1‑149 (MQ=255)
                                                        gTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtatta                                                                                         >  1:1304030/1‑149 (MQ=255)
                                                         tGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCNCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtattat                                                                                        <  1:2082918/149‑1 (MQ=255)
                                                         tGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtattat                                                                                        <  2:2393458/149‑1 (MQ=255)
                                                         tGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtattat                                                                                        <  1:1296391/149‑1 (MQ=255)
                                                           ccGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCg                                                                                      >  2:1961100/1‑149 (MQ=255)
                                                            cGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGc                                                                                     <  2:2813551/149‑1 (MQ=255)
                                                               tACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCAc                                                                                  >  2:3137152/1‑149 (MQ=255)
                                                               tACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCAc                                                                                  >  2:4114467/1‑149 (MQ=255)
                                                                 cGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTg                                                                                >  1:2930021/1‑149 (MQ=255)
                                                                       tGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAg                                                                          <  1:142032/149‑1 (MQ=255)
                                                                       tGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAg                                                                          <  1:3090320/149‑1 (MQ=255)
                                                                       tGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAg                                                                          <  2:1524359/149‑1 (MQ=255)
                                                                         gCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGc                                                                        >  2:3932599/1‑149 (MQ=255)
                                                                          cGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCt                                                                       <  2:2379646/149‑1 (MQ=255)
                                                                              aGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGaa                                                                   <  2:927133/149‑1 (MQ=255)
                                                                               gCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAc                                                                  >  1:602/1‑149 (MQ=255)
                                                                               gCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAc                                                                  <  2:2358240/149‑1 (MQ=255)
                                                                                cGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAcc                                                                 <  1:566522/149‑1 (MQ=255)
                                                                                cGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAcc                                                                 <  1:3124396/149‑1 (MQ=255)
                                                                                cGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCCCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAcc                                                                 <  1:1069604/149‑1 (MQ=255)
                                                                                   cGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGt                                                              <  1:3571641/149‑1 (MQ=255)
                                                                                   cGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGt                                                              <  1:670734/149‑1 (MQ=255)
                                                                                    ggCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTa                                                             <  2:1305931/149‑1 (MQ=255)
                                                                                      ccTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATg                                                           <  2:2436791/149‑1 (MQ=255)
                                                                                       cTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGt                                                          <  2:332033/149‑1 (MQ=255)
                                                                                       cTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGt                                                          <  1:3882703/149‑1 (MQ=255)
                                                                                         cGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGttt                                                        <  2:3551774/149‑1 (MQ=255)
                                                                                           gCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGc                                                      <  1:2224542/149‑1 (MQ=255)
                                                                                           gCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGc                                                      >  1:1582783/1‑149 (MQ=255)
                                                                                                gAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTaa                                                 >  2:4085712/1‑149 (MQ=255)
                                                                                                 aaCTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAAc                                                <  1:2545768/149‑1 (MQ=255)
                                                                                                 aaCTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAAc                                                >  2:1234348/1‑149 (MQ=255)
                                                                                                        aaaCCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCa                                         >  2:3288593/1‑149 (MQ=255)
                                                                                                            cGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGaa                                     >  2:3374356/1‑149 (MQ=255)
                                                                                                            cGCTACAAAGGCATCCGTACGACGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCGACTGAAAAAGGCTGGAACCTGTCTGTGTGCGGTAACGGCGGCATGaa                                     >  2:2221078/1‑149 (MQ=255)
                                                                                                             gCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGaaa                                    >  1:2324472/1‑149 (MQ=255)
                                                                                                             gCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGaaa                                    <  2:3486646/149‑1 (MQ=255)
                                                                                                             gCTACAAAGGCATCCGCACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGaaa                                    <  2:1070014/149‑1 (MQ=255)
                                                                                                              cTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAAc                                   <  1:4088919/149‑1 (MQ=255)
                                                                                                              cTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAAc                                   <  1:1227088/149‑1 (MQ=255)
                                                                                                              cTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAAc                                   <  2:3651075/149‑1 (MQ=255)
                                                                                                                 cAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACcgc                                >  1:513579/1‑149 (MQ=255)
                                                                                                                  aaaGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACcgcg                               <  2:756949/149‑1 (MQ=255)
                                                                                                                  aaaGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACcgcg                               >  1:2146270/1‑149 (MQ=255)
                                                                                                                  aaaGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACcgcg                               >  2:1749935/1‑149 (MQ=255)
                                                                                                                     ggCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCa                            >  2:74340/1‑149 (MQ=255)
                                                                                                                        aTCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGc                         >  2:267553/1‑149 (MQ=255)
                                                                                                                                gccgcACAAAATGAAGTGAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATctgc                 <  2:2966787/149‑1 (MQ=255)
                                                                                                                                     acaAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggcg            >  2:841651/1‑149 (MQ=255)
                                                                                                                                       aaaaTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggcggc          >  1:1898465/1‑149 (MQ=255)
                                                                                                                                       aaaaTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggcggc          <  2:2265249/149‑1 (MQ=255)
                                                                                                                                        taaTGAAGTTAGGTGTCTCCGGGTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGTCGTGGGTATTTTCGCCACTGAAATAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGgcggcg         <  2:2168174/148‑1 (MQ=255)
                                                                                                                                          aTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGa       <  1:2665119/149‑1 (MQ=255)
                                                                                                                                          aTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGa       >  1:3795070/1‑149 (MQ=255)
                                                                                                                                           tGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGat      >  2:3791302/1‑149 (MQ=255)
                                                                                                                                           tGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGat      <  2:1641625/149‑1 (MQ=255)
                                                                                                                                           tGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGat      >  1:208787/1‑149 (MQ=255)
                                                                                                                                           tGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGat      >  1:966026/1‑149 (MQ=255)
                                                                                                                                           tGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGat      >  2:3250595/1‑149 (MQ=255)
                                                                                                                                              aGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATatc   <  2:1955316/149‑1 (MQ=255)
                                                                                                                                               gTTTGGTGTCTACGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGTGTATTATCGCCACTGAAAAAGGCTGGAACATGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATatcg  <  1:3948459/149‑1 (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: