Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 8,261 T→G 29.7% L8L (CTT→CTG talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,2610TG29.7% 140.8 / 90.3 91L8L (CTT→CTGtalBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base T (36/28);  new base G (12/15);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

GAAATCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCT  >  CP009273/8112‑8409
                                                                                                                                                     |                                                                                                                                                    
gAAATCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCtt                                                                                                                                                      <  1:61454/150‑1 (MQ=255)
 aaaTCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGTGCTGACGTTGCGTCCTGATATAATTAGGGCCGACCGGTTTTTTCCCCCTCACAAGCTGTTTTAAGAGGAATACTATCATGACGGACAAATTGACCTCCCttt                                                                                                                                                     >  2:975754/1‑149 (MQ=255)
   aTCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGGTATCATCAGGGCAGACCGGGTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGt                                                                                                                                                   >  2:667214/1‑150 (MQ=255)
        aaaCAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTa                                                                                                                                              <  1:660055/150‑1 (MQ=255)
         aaCAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTa                                                                                                                                              >  1:519727/1‑149 (MQ=255)
             gTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTAcacca                                                                                                                                         <  1:389418/150‑1 (MQ=255)
                 cGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAAATGACCTCCCTTCGTCAGTACACCACCg                                                                                                                                      >  1:893918/1‑149 (MQ=255)
                  gCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACaccacc                                                                                                                                       <  1:467866/147‑1 (MQ=255)
                   ccTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAATAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAg                                                                                                                                   >  1:1000076/1‑150 (MQ=255)
                    cTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTTCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCCTACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCCTAGt                                                                                                                                  >  1:708670/1‑150 (MQ=255)
                    cTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGt                                                                                                                                  <  2:793423/150‑1 (MQ=255)
                      ggCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTgg                                                                                                                                >  1:396905/1‑150 (MQ=255)
                      ggCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTgg                                                                                                                                >  1:90980/1‑150 (MQ=255)
                       gCGATAACCGTCTTGTCGGCGTTTGCGCTGACGTTACGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAGCAAGCTGTTTAAAGAGAAATACTATCATGACGTACAAATTGACCAAGCTTCGTCAGTACACGGCCGTAGTg                                                                                                                                 <  2:163907/148‑1 (MQ=255)
                         gATAACCGTCTTGTCGGCGGTTTCGCTGACGTTGCGTCGTGTTATTTTTTTTGCTTTCCGGTTTCTTTCCCCTTACAAGCTTTTTTAATAGAAATTCTTTCATTTCGGTCAAATTTACCTCCCTTCTTCATTACACCACCGTAGTGGCCg                                                                                                                             >  1:206007/1‑150 (MQ=255)
                          aTAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGGTATCATCCGGGCAGACCGGTTACATCCCCCTAACAAACTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGa                                                                                                                            >  2:538235/1‑150 (MQ=255)
                           tAACCGTCTTGTCGGCGGCTCTGCTGCCGCTGCGTCCTGATATCACCAGGGGAGACCGGTTTCCTCCCCCTAAAAAGCTGTTTTAAGAGAAATACTATCATGACGGACAAATTGACCTTCCCTCGCCAGGACACCACCGAAGTGGCCgc                                                                                                                            >  2:987542/1‑148 (MQ=255)
                             aCCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGacac                                                                                                                         <  2:780705/150‑1 (MQ=255)
                             aCCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGaca                                                                                                                          <  2:147352/149‑1 (MQ=255)
                                gTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTgg                                                                                                                      <  2:350647/150‑1 (MQ=255)
                                           ggTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCgg                                                                                                           <  2:582330/150‑1 (MQ=255)
                                                     aCGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAg                                                                                                   <  1:363010/148‑1 (MQ=255)
                                                           cgtcgtGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATc                                                                                            >  2:132322/1‑149 (MQ=255)
                                                           cgtcgtGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATc                                                                                            >  1:62772/1‑149 (MQ=255)
                                                              cgtGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAAcc                                                                                        >  2:75391/1‑150 (MQ=255)
                                                              cgtGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAAcc                                                                                        >  1:954190/1‑150 (MQ=255)
                                                                tGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCg                                                                                       >  2:378573/1‑149 (MQ=255)
                                                                tGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGc                                                                                      <  2:917768/150‑1 (MQ=255)
                                                                  atatCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAg                                                                                    <  1:795639/150‑1 (MQ=255)
                                                                        tcaGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGcc                                                                              <  1:83526/150‑1 (MQ=255)
                                                                              cAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCAcaacc                                                                        <  1:23789/150‑1 (MQ=255)
                                                                               aGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCAcaacc                                                                        >  2:154419/1‑149 (MQ=255)
                                                                               aGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCAcaacc                                                                        >  2:789816/1‑149 (MQ=255)
                                                                                 aCCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCAcaaccaac                                                                     >  1:596127/1‑150 (MQ=255)
                                                                                  ccGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACaaccaacc                                                                    <  1:862905/150‑1 (MQ=255)
                                                                                     gTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCATAACCAACCCtt                                                                 >  1:68190/1‑150 (MQ=255)
                                                                                         cATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGACGCCACAACCAACCCTtctct                                                             >  2:108836/1‑150 (MQ=255)
                                                                                          aTCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTtctct                                                             >  2:287614/1‑149 (MQ=255)
                                                                                           tCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCa                                                           >  2:311304/1‑150 (MQ=255)
                                                                                            cccccTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCAt                                                          <  2:206007/150‑1 (MQ=255)
                                                                                             ccccTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCAtt                                                         >  1:264520/1‑150 (MQ=255)
                                                                                               ccTAACAAGCTGTTTAAATAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCt                                                       <  2:1000076/150‑1 (MQ=255)
                                                                                                  aaCAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACa                                                                                                                 >  1:559161/1‑89 (MQ=255)
                                                                                                   aCAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACAt                                                                                                                <  2:559161/89‑1 (MQ=255)
                                                                                                           tttAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAAcgcagcg                                             >  2:612642/1‑148 (MQ=255)
                                                                                                            ttAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAAcgcagcg                                             >  2:758/1‑147 (MQ=255)
                                                                                                                agagAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTc                                      >  2:979953/1‑150 (MQ=255)
                                                                                                                 gagaAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTcc                                     <  2:417583/150‑1 (MQ=255)
                                                                                                                 gagaAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTcc                                     <  1:287614/150‑1 (MQ=255)
                                                                                                                   gaAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCg                                    >  1:250908/1‑149 (MQ=255)
                                                                                                                    aaaTACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGa                                  <  2:304990/150‑1 (MQ=255)
                                                                                                                       acacaaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTa                               <  2:643746/143‑1 (MQ=255)
                                                                                                                       acacaaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTa                               <  1:953207/143‑1 (MQ=255)
                                                                                                                       acacaaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTa                               <  1:859845/143‑1 (MQ=255)
                                                                                                                         cTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATAcc                             <  1:975754/150‑1 (MQ=255)
                                                                                                                          caaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTAcca                            <  1:166744/146‑2 (MQ=255)
                                                                                                                          caaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCg                            <  2:674975/146‑1 (MQ=255)
                                                                                                                          caaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCg                            <  2:813193/146‑1 (MQ=255)
                                                                                                                           aaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                           <  2:679568/147‑1 (MQ=255)
                                                                                                                           aTCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCg                            >  2:305898/1‑149 (MQ=255)
                                                                                                                             aaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                           >  2:658169/2‑148 (MQ=255)
                                                                                                                             aaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                           >  1:562561/2‑148 (MQ=255)
                                                                                                                             aaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                           >  1:649239/2‑148 (MQ=255)
                                                                                                                             aaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                           >  1:464207/2‑148 (MQ=255)
                                                                                                                              aTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                           >  1:761711/1‑147 (MQ=255)
                                                                                                                              aTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAg                        <  1:266261/150‑1 (MQ=255)
                                                                                                                              aTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAg                        <  1:43749/150‑1 (MQ=255)
                                                                                                                              aTGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAg                        <  1:667214/150‑1 (MQ=255)
                                                                                                                               tGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAgc                       <  2:855228/150‑2 (MQ=255)
                                                                                                                               tGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAgc                       <  2:26879/150‑2 (MQ=255)
                                                                                                                               tGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAg                        <  2:761711/149‑1 (MQ=255)
                                                                                                                               tGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGt                       >  1:403580/1‑150 (MQ=255)
                                                                                                                                 aCGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAgct                      >  2:206074/1‑147 (MQ=255)
                                                                                                                                 aCGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAgct                      >  2:75551/1‑147 (MQ=255)
                                                                                                                                  cGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctga                    >  2:839585/1‑150 (MQ=255)
                                                                                                                                  cGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctga                    >  1:1031049/1‑150 (MQ=255)
                                                                                                                                  cGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctga                    >  2:986790/1‑150 (MQ=255)
                                                                                                                                  cGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGttga                    >  1:455146/1‑150 (MQ=255)
                                                                                                                                   ggATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgat                   <  1:477985/150‑1 (MQ=255)
                                                                                                                                   ggACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGttgat                   >  2:679368/1‑150 (MQ=255)
                                                                                                                                    gACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGttgatt                  <  2:578449/150‑1 (MQ=255)
                                                                                                                                     aCAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGttgatt                  >  1:307169/1‑149 (MQ=255)
                                                                                                                                      taaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgattga                >  1:7525/2‑150 (MQ=255)
                                                                                                                                      taaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgattga                >  1:208610/2‑150 (MQ=255)
                                                                                                                                      taaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgattga                <  2:317463/149‑1 (MQ=255)
                                                                                                                                        aaTTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTtgattgat               >  2:975469/1‑149 (MQ=255)
                                                                                                                                        aaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgattgac               <  2:994661/149‑2 (MQ=255)
                                                                                                                                          ttGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTgatgat            <  2:396905/150‑1 (MQ=255)
                                                                                                                                          ttGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTgatgat            <  1:75391/150‑1 (MQ=255)
                                                                                                                                           tGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGatgatg           >  2:392195/1‑150 (MQ=255)
                                                                                                                                             aCCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGc          <  2:447864/149‑1 (MQ=255)
                                                                                                                                             aCCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCt         <  2:556892/150‑1 (MQ=255)
                                                                                                                                              ccTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTg        >  1:290390/1‑150 (MQ=255)
                                                                                                                                                    ttCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCt  >  1:1022081/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAAATCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCT  >  CP009273/8112‑8409

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: