Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,936,836 C→A 100% A47A (GCC→GCA csdA → cysteine sulfinate desulfinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,936,8360CA99.6% 739.4 / ‑4.3 226A47A (GCC→GCAcsdAcysteine sulfinate desulfinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (105/119);  total (105/121)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-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 226 total aligned reads displayed.

GTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAAC  >  CP009273/2936691‑2936980
                                                                                                                                                 |                                                                                                                                                
gTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAgg                                                                                                                                                >  2:3289426/1‑148 (MQ=255)
  aCCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGaa                                                                                                                                              <  1:762326/148‑1 (MQ=255)
  aCCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGaa                                                                                                                                              <  1:1423021/148‑1 (MQ=255)
    cATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAAc                                                                                                                                            >  2:2347158/1‑148 (MQ=255)
      tGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGt                                                                                                                                          >  1:2377946/1‑148 (MQ=255)
       gAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTc                                                                                                                                         >  2:2224500/1‑148 (MQ=255)
        aaCGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTcc                                                                                                                                        <  2:316124/148‑1 (MQ=255)
            tttttAATCCCGCGCAGTTTCGCGCCAAGTTTCCCGCACTACAGGATGCTGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGACGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCg                                                                                                                                    <  1:1322856/148‑1 (MQ=255)
                   tCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAg                                                                                                                             <  1:733108/148‑1 (MQ=255)
                     ccGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGtt                                                                                                                           >  1:965508/1‑148 (MQ=255)
                          cAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCg                                                                                                                      >  2:3329507/1‑148 (MQ=255)
                             tttCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCCGTTTGCCGAAg                                                                                                                   >  1:711844/1‑148 (MQ=255)
                             tttCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAg                                                                                                                   >  2:2422775/1‑148 (MQ=255)
                             tttCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAg                                                                                                                   >  1:459722/1‑148 (MQ=255)
                             tttCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAg                                                                                                                   >  2:790094/1‑148 (MQ=255)
                              ttCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGc                                                                                                                  >  2:247136/1‑148 (MQ=255)
                              ttCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGc                                                                                                                  >  2:247175/1‑148 (MQ=255)
                              ttCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGc                                                                                                                  >  2:3224524/1‑148 (MQ=255)
                                     ccAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGc                                                                                                           >  2:2504337/1‑148 (MQ=255)
                                     ccAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGc                                                                                                           >  2:3761758/1‑148 (MQ=255)
                                      cAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACTTCCATCGCAGCCAGTTTGCCGAAGCCCAACGcc                                                                                                          <  2:3391299/148‑1 (MQ=255)
                                      cAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGcc                                                                                                          <  1:1918174/148‑1 (MQ=255)
                                         tttCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGa                                                                                                       >  2:468841/1‑148 (MQ=255)
                                          ttCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGAc                                                                                                      >  1:1472433/1‑148 (MQ=255)
                                              cGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACcgcg                                                                                                  <  2:1018414/148‑1 (MQ=255)
                                                 aCTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGt                                                                                               <  1:259131/148‑1 (MQ=255)
                                                  cTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGtt                                                                                              >  2:2854098/1‑148 (MQ=255)
                                                  cTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGtt                                                                                              >  2:3586506/1‑148 (MQ=255)
                                                   tACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTa                                                                                             >  1:209807/1‑148 (MQ=255)
                                                     cAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATg                                                                                           <  1:2685764/148‑1 (MQ=255)
                                                     cAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATg                                                                                           <  1:3063196/148‑1 (MQ=255)
                                                      aGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGa                                                                                          <  2:3400505/148‑1 (MQ=255)
                                                       ggATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGaa                                                                                         >  1:983691/1‑148 (MQ=255)
                                                        gATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTCTGAAg                                                                                        >  2:1906127/1‑148 (MQ=255)
                                                        gATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAg                                                                                        >  1:3490519/1‑148 (MQ=255)
                                                        gATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAg                                                                                        <  1:2882086/148‑1 (MQ=255)
                                                         aTGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGc                                                                                       <  1:2126133/148‑1 (MQ=255)
                                                         aTGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGc                                                                                       <  1:197970/148‑1 (MQ=255)
                                                               gCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACg                                                                                 <  1:3168331/148‑1 (MQ=255)
                                                                   cTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACgagag                                                                             >  1:1174570/1‑148 (MQ=255)
                                                                   cTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACgagag                                                                             >  2:2648149/1‑148 (MQ=255)
                                                                   cTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACgagag                                                                             >  1:381282/1‑148 (MQ=255)
                                                                   cTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACgagag                                                                             >  2:73504/1‑148 (MQ=255)
                                                                     aTCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGaa                                                                           <  2:1454051/148‑1 (MQ=255)
                                                                      tctcGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGaaa                                                                          <  1:3457181/148‑1 (MQ=255)
                                                                      tctcGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGaaa                                                                          <  1:1837667/148‑1 (MQ=255)
                                                                      tctcGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGaaa                                                                          <  1:264269/148‑1 (MQ=255)
                                                                      tctcGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGaaa                                                                          <  2:704224/148‑1 (MQ=255)
                                                                       ctcGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAg                                                                         <  1:3289426/148‑1 (MQ=255)
                                                                       ctcGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAg                                                                         >  2:3012133/1‑148 (MQ=255)
                                                                         cGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTg                                                                       <  2:242454/148‑1 (MQ=255)
                                                                          gACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTgg                                                                      <  2:2456940/148‑1 (MQ=255)
                                                                               cgccgcGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCaa                                                                 >  2:3515114/1‑148 (MQ=255)
                                                                               cgccgcGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCaa                                                                 >  1:3536644/1‑148 (MQ=255)
                                                                               cgccgcGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCaa                                                                 <  2:765785/148‑1 (MQ=255)
                                                                                gccgcGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAAt                                                                <  2:2381702/148‑1 (MQ=255)
                                                                                 ccgcGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGGAGCTGCACGAGAGAAAGTGGCGCAAtt                                                               <  1:2761248/148‑1 (MQ=255)
                                                                                 ccgcGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAAtt                                                               >  1:333825/1‑148 (MQ=255)
                                                                                 ccgcGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAAtt                                                               <  2:874136/148‑1 (MQ=255)
                                                                                 ccgcGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCCCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGGAGCTGCACGAGAGAAAGTGGCGCAAtt                                                               <  1:2466837/148‑1 (MQ=255)
                                                                                  cgcgACCGCGCTTACACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTa                                                              <  2:630323/148‑1 (MQ=255)
                                                                                  cgcgACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTa                                                              <  1:2347158/148‑1 (MQ=255)
                                                                                  cgcgACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCAAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTa                                                              <  1:424260/148‑1 (MQ=255)
                                                                                    cgACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACt                                                            >  1:1246597/1‑148 (MQ=255)
                                                                                    cgACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACt                                                            >  2:2182440/1‑148 (MQ=255)
                                                                                       ccGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGaa                                                         >  2:2139283/1‑148 (MQ=255)
                                                                                          cgcTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGACACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCCCGAGAGCAAGTGGCGCAATTACTGACTGc                                                      >  2:2870407/1‑148 (MQ=255)
                                                                                          cgcTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGc                                                      <  2:537133/148‑1 (MQ=255)
                                                                                          cgcTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGc                                                      >  1:3139458/1‑148 (MQ=255)
                                                                                          cgcTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGc                                                      >  1:3258100/1‑148 (MQ=255)
                                                                                                aaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGga                                                >  2:19306/1‑148 (MQ=255)
                                                                                                aaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGga                                                <  1:2719683/148‑1 (MQ=255)
                                                                                                aaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGga                                                <  2:3067871/148‑1 (MQ=255)
                                                                                                aaCCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGga                                                >  2:1046129/1‑148 (MQ=255)
                                                                                                     gAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATa                                           >  2:1801497/1‑148 (MQ=255)
                                                                                                       aGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATaaa                                         <  1:2456944/148‑1 (MQ=255)
                                                                                                       aGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATaaa                                         <  2:448195/148‑1 (MQ=255)
                                                                                                          cGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACt                                      <  1:1672550/148‑1 (MQ=255)
                                                                                                          cGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACt                                      <  1:2417372/148‑1 (MQ=255)
                                                                                                          cGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACt                                      <  1:2417680/148‑1 (MQ=255)
                                                                                                           gTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTa                                     <  2:277080/148‑1 (MQ=255)
                                                                                                            tGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTAt                                    >  2:2588805/1‑148 (MQ=255)
                                                                                                                   aGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTgg                             <  2:1128669/148‑1 (MQ=255)
                                                                                                                   aGCCACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTgg                             <  2:1967101/148‑1 (MQ=255)
                                                                                                                     ccACCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGAc                           >  2:2128261/1‑148 (MQ=255)
                                                                                                                       aCCCAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGAcgc                         >  2:1573022/1‑148 (MQ=255)
                                                                                                                        cccAACAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGAcgcg                        >  1:3762461/1‑148 (MQ=255)
                                                                                                                           aaCAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCgg                     <  2:5585/148‑1 (MQ=255)
                                                                                                                            aCAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGc                    >  1:1768843/1‑148 (MQ=255)
                                                                                                                             cAGTTTTACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGca                   >  1:666152/1‑148 (MQ=255)
                                                                                                                                ttttACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGcacca                >  1:898474/1‑148 (MQ=255)
                                                                                                                                 tttACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCACTTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGcaccac               >  1:3571997/1‑148 (MQ=255)
                                                                                                                                  ttACAGTCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACt              <  2:3641008/148‑1 (MQ=255)
                                                                                                                                        tCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATcc        >  2:1326818/1‑148 (MQ=255)
                                                                                                                                        tCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATcc        >  2:1547626/1‑148 (MQ=255)
                                                                                                                                        tCTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATcc        >  2:2621631/1‑148 (MQ=255)
                                                                                                                                         cTGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCa       <  2:752009/148‑1 (MQ=255)
                                                                                                                                          tGAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCAt      <  1:862356/148‑1 (MQ=255)
                                                                                                                                           gAGCGCAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATc     >  2:3106427/1‑148 (MQ=255)
                                                                                                                                              cgcAGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAAc  <  1:2504337/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCACCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAAC  >  CP009273/2936691‑2936980

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: