Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,927,123 C→A 16.6% T138K (ACG→AAG)  rbsD → putative cytoplasmic sugar‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,927,1230CA16.6% 788.3 / 101.9 265T138K (ACG→AAG) rbsDputative cytoplasmic sugar‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (115/105);  new base A (31/14);  total (146/119)
Fisher's exact test for biased strand distribution p-value = 4.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 267 total aligned reads displayed.

TGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCA  >  CP009273/3926974‑3927270
                                                                                                                                                     |                                                                                                                                                   
 gAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACg                                                                                                                                                    <  1:5099394/150‑1 (MQ=255)
 gAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACg                                                                                                                                                    <  1:3881521/150‑1 (MQ=255)
  agcagcTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGt                                                                                                                                                   <  2:4635325/150‑1 (MQ=255)
   gcagcTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGtt                                                                                                                                                  >  1:1155367/1‑150 (MQ=255)
     agcTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCt                                                                                                                                                >  1:2235009/1‑150 (MQ=255)
       cTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGa                                                                                                                                              <  1:1540124/150‑1 (MQ=255)
           aaaaaCACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGAAGTTCTGAGGcc                                                                                                                                          <  2:3282129/150‑1 (MQ=255)
            aaaaCACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCg                                                                                                                                         >  1:1979673/1‑150 (MQ=255)
             aaaCACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGt                                                                                                                                        <  2:1556873/150‑1 (MQ=255)
              aaCACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGAAGTTCTGAGGCCGTc                                                                                                                                       >  2:4740217/1‑150 (MQ=255)
                cacCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCAt                                                                                                                                     <  2:1504268/150‑1 (MQ=255)
                cacCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCAt                                                                                                                                     >  2:2266952/1‑150 (MQ=255)
                 acCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGAAGTTCTGAGGCCGTCATg                                                                                                                                    >  2:5133000/1‑150 (MQ=255)
                     gggAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAg                                                                                                                                >  2:1371451/1‑150 (MQ=255)
                      ggAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGc                                                                                                                               >  2:1036542/1‑150 (MQ=255)
                               aTTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCa                                                                                                                      <  1:4217508/150‑1 (MQ=255)
                               aTTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCa                                                                                                                      >  1:778276/1‑150 (MQ=255)
                                ttGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAg                                                                                                                     <  2:2111975/150‑1 (MQ=255)
                                ttGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAg                                                                                                                     <  1:3767242/150‑1 (MQ=255)
                                ttGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAg                                                                                                                     <  2:2868507/150‑1 (MQ=255)
                                  gAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGAAGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCt                                                                                                                   >  2:699107/1‑150 (MQ=255)
                                   aaaTTCGTTACACCACGCATGAACAATTCAAACAAAAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATACTCTGTGCTGGCGTGAAGTTCTGAGGCCGTCATGGAAGCATTACTTCAGcgt                                                                                                                  >  1:3411556/1‑148 (MQ=255)
                                   aaaTTAGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGAAGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCtt                                                                                                                  <  1:1010989/150‑1 (MQ=255)
                                           tACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCAt                                                                                                          >  1:1331409/1‑150 (MQ=255)
                                            acacCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATc                                                                                                         >  1:2610411/1‑150 (MQ=255)
                                            acacCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGGATGCGAATATCATTCTCTGTGCTGGCGTGACGGTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATc                                                                                                         >  2:5107265/1‑150 (MQ=255)
                                            acacCACGCATGAACAATTCAAACAACAAACCGCAGAAAGACAGGAGGTAATTCGCAGAGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTTCTGGCGTGAAGTTCTGAGGCCTTCATGGAAGCATTACTTCAGCTTAAAGGCATc                                                                                                         >  1:32249/1‑150 (MQ=255)
                                              accacGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCAGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGa                                                                                                       <  1:720909/150‑1 (MQ=255)
                                                     aTGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGcc                                                                                                >  2:2766483/1‑150 (MQ=255)
                                                       gAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCtt                                                                                              >  1:3101523/1‑150 (MQ=255)
                                                       gAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCtt                                                                                              >  1:2492631/1‑150 (MQ=255)
                                                       gAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCtt                                                                                              <  1:4357581/150‑1 (MQ=255)
                                                         aCAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTcc                                                                                            >  2:626897/1‑150 (MQ=255)
                                                         aCAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTcc                                                                                            <  1:2806489/150‑1 (MQ=255)
                                                             ttCAAACAACAAACCTCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCggg                                                                                        <  1:4172628/150‑1 (MQ=255)
                                                             ttCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCggg                                                                                        >  2:1914609/1‑150 (MQ=255)
                                                             ttCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCggg                                                                                        <  2:4588405/150‑1 (MQ=255)
                                                              tCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGc                                                                                       >  2:66184/1‑150 (MQ=255)
                                                                 aacaacAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTa                                                                                    <  2:2993171/150‑1 (MQ=255)
                                                                   caacaaACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTaaa                                                                                  >  1:2989395/1‑150 (MQ=255)
                                                                   caacaaACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTaaa                                                                                  <  1:4669668/150‑1 (MQ=255)
                                                                    aacaaACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTaaaa                                                                                 >  2:1981509/1‑150 (MQ=255)
                                                                    aacaaACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTaaaa                                                                                 >  1:5172551/1‑150 (MQ=255)
                                                                     acaaACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAg                                                                                >  1:829961/1‑150 (MQ=255)
                                                                       aaaCCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGcc                                                                              >  2:4373018/1‑150 (MQ=255)
                                                                       aaaCCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGcc                                                                              >  2:4543749/1‑150 (MQ=255)
                                                                       aaaCCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGcc                                                                              >  2:1490065/1‑150 (MQ=255)
                                                                        aaCCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGccc                                                                             >  2:2067364/1‑150 (MQ=255)
                                                                        aaCCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGccc                                                                             >  2:2066940/1‑150 (MQ=255)
                                                                        aaCCGCAGAAAGTCAGGCGGTAATTCGCAACGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGccc                                                                             <  1:595051/150‑1 (MQ=255)
                                                                               gAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCggg                                                                      >  1:1528050/1‑150 (MQ=255)
                                                                                 aaGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGgcg                                                                    >  2:1067100/1‑150 (MQ=255)
                                                                                 aaGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGgcg                                                                    <  2:378631/150‑1 (MQ=255)
                                                                                 aaGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGgcg                                                                    <  2:3022805/150‑1 (MQ=255)
                                                                                 aaGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGAAGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGgcg                                                                    >  1:3040362/1‑150 (MQ=255)
                                                                                   gTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCa                                                                  >  2:1155492/1‑150 (MQ=255)
                                                                                      aGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCg                                                               <  2:3580835/150‑1 (MQ=255)
                                                                                      aGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCg                                                               >  1:5004258/1‑150 (MQ=255)
                                                                                        gCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGtt                                                             <  1:2332962/150‑1 (MQ=255)
                                                                                           gTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTaaa                                                          >  1:2342047/1‑150 (MQ=255)
                                                                                           gTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTaaa                                                          <  1:4956029/150‑1 (MQ=255)
                                                                                            tAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAAt                                                         <  2:3427450/150‑1 (MQ=255)
                                                                                            tAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAAt                                                         <  1:2266952/150‑1 (MQ=255)
                                                                                               ttCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTc                                                      >  1:2948244/1‑150 (MQ=255)
                                                                                                tCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCt                                                     >  1:4546417/1‑150 (MQ=255)
                                                                                                tCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCt                                                     <  1:3543390/150‑1 (MQ=255)
                                                                                                  gCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTAt                                                   >  2:4188223/1‑150 (MQ=255)
                                                                                                        gagaATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGc                                             <  2:78410/150‑1 (MQ=255)
                                                                                                          gaATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCcg                                           >  2:2342112/1‑150 (MQ=255)
                                                                                                           aaTGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCcgc                                          <  1:1048863/150‑1 (MQ=255)
                                                                                                           aaTGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCcgc                                          >  2:5047120/1‑150 (MQ=255)
                                                                                                           aaTGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCcgc                                          <  1:4270612/150‑1 (MQ=255)
                                                                                                                   ccGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGgc                                  >  2:2034345/1‑150 (MQ=255)
                                                                                                                     gTATGCGAATATCATTCTCTGTGCTGGCGTGAAGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGgcgc                                >  2:4922020/1‑150 (MQ=255)
                                                                                                                      tATGCGAATATCATTCTCTGTGCTGGCGTGAAGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCt                               <  1:2054753/150‑1 (MQ=255)
                                                                                                                        tGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCtgg                             >  2:1130833/1‑150 (MQ=255)
                                                                                                                        tGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCtgg                             >  2:4975059/1‑150 (MQ=255)
                                                                                                                         gCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCtggt                            <  1:881843/150‑1 (MQ=255)
                                                                                                                          cGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCtggtg                           >  1:93778/1‑150 (MQ=255)
                                                                                                                          cGAATATCATTCTCTGTGCTGGCGTGAAGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCtggtg                           <  1:624375/150‑1 (MQ=255)
                                                                                                                             atatCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGc                        >  2:2506120/1‑150 (MQ=255)
                                                                                                                             atatCATTCTCTGTGCTGGCGTGAAGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGc                        >  2:3823706/1‑150 (MQ=255)
                                                                                                                               atCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGa                      >  2:3921177/1‑150 (MQ=255)
                                                                                                                                tCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGaa                     >  2:4418368/1‑150 (MQ=255)
                                                                                                                                 cATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGaaa                    <  2:2607219/150‑1 (MQ=255)
                                                                                                                                  aTTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGaaaa                   >  1:2185656/1‑150 (MQ=255)
                                                                                                                                  aTTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGaaaa                   >  1:343433/1‑150 (MQ=255)
                                                                                                                                  aTTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGaaaa                   >  1:4781552/1‑150 (MQ=255)
                                                                                                                                   ttCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAAc                  <  1:1914609/150‑1 (MQ=255)
                                                                                                                                     ctctGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGAGGGCGAAAACgg                >  1:1759929/1‑150 (MQ=255)
                                                                                                                                        tgtgCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGgcgc             <  1:1658433/150‑1 (MQ=255)
                                                                                                                                        tgtgCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGgcgc             <  1:1658032/150‑1 (MQ=255)
                                                                                                                                        tgtgCTGGCGTGAAGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGgcgc             >  1:1435655/1‑150 (MQ=255)
                                                                                                                                            cTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGt         >  1:2768637/1‑150 (MQ=255)
                                                                                                                                              ggCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTaa       >  1:2657533/1‑150 (MQ=255)
                                                                                                                                              ggCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTaa       >  1:127920/1‑150 (MQ=255)
                                                                                                                                              ggCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTaa       >  1:1977969/1‑150 (MQ=255)
                                                                                                                                              ggCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTaa       >  2:1282186/1‑150 (MQ=255)
                                                                                                                                                 gTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATc    <  2:574433/150‑1 (MQ=255)
                                                                                                                                                  tGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATcc   >  1:999511/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCA  >  CP009273/3926974‑3927270

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: