Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 8,258 C→T 29.4% S7S (TCC→TCT)  talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,2580CT29.4% 117.1 / 80.1 93S7S (TCC→TCT) talBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base C (37/28);  new base T (12/15);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

GAAATCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTG  >  CP009273/8112‑8403
                                                                                                                                                  |                                                                                                                                                 
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                                                                                  >  1:954190/1‑150 (MQ=255)
                                                              cgtGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAAcc                                                                                  >  2:75391/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:789816/1‑149 (MQ=255)
                                                                               aGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCAcaacc                                                                  >  2:154419/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                         <  1:953207/143‑1 (MQ=255)
                                                                                                                       acacaaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTa                         <  1:859845/143‑1 (MQ=255)
                                                                                                                       acacaaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTa                         <  2:643746/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                     >  1:649239/2‑148 (MQ=255)
                                                                                                                             aaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                     >  2:658169/2‑148 (MQ=255)
                                                                                                                             aaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                     >  1:464207/2‑148 (MQ=255)
                                                                                                                             aaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                     >  1:562561/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:26879/150‑2 (MQ=255)
                                                                                                                               tGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAgc                 <  2:855228/150‑2 (MQ=255)
                                                                                                                               tGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAg                  <  2:761711/149‑1 (MQ=255)
                                                                                                                               tGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGt                 >  1:403580/1‑150 (MQ=255)
                                                                                                                                 aCGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAgct                >  2:75551/1‑147 (MQ=255)
                                                                                                                                 aCGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAgct                >  2:206074/1‑147 (MQ=255)
                                                                                                                                  cGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctga              >  2:986790/1‑150 (MQ=255)
                                                                                                                                  cGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctga              >  2:839585/1‑150 (MQ=255)
                                                                                                                                  cGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctga              >  1:1031049/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          <  2:317463/149‑1 (MQ=255)
                                                                                                                                      taaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgattga          >  1:7525/2‑150 (MQ=255)
                                                                                                                                      taaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgattga          >  1:208610/2‑150 (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)
                                                                                                                                                  |                                                                                                                                                 
GAAATCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTG  >  CP009273/8112‑8403

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: ‑