Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,173,864 G→T 4.8% E29* (GAG→TAG) ‡ ygiC → ATP‑Grasp family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,173,8640GT 4.8% 276.2 / 3.4 95E29* (GAG→TAG) ygiCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base G (44/46);  new base T (2/3);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

AGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCC  >  CP009273/3173721‑3174008
                                                                                                                                               |                                                                                                                                                 
aGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCg                                                                                                                                              >  1:609132/1‑149 (MQ=255)
  cACTATGCAGCGTAGTGCAACCGGTGCCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTAAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGACCACGAATACGGTTTCAATTTTCACACCATGTAGGGCgtgtcgag                                                                                                                                            >  1:486838/1‑142 (MQ=255)
        gCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTACATGgtg                                                                                                                                     <  1:1125187/149‑1 (MQ=255)
           gggTAGTGCAACCGGTCCCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTCCCGAGCGCCCGGCCTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑Tg                                                                                                                                        <  2:264960/141‑1 (MQ=255)
           gCGTAGTGCAACCGGTACCTCTTCCCGTTCAGTGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTAcatg                                                                                                                                        <  2:1137273/143‑4 (MQ=255)
            cGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGGGCCGTAC‑Tgg                                                                                                                                       >  2:121788/1‑143 (MQ=255)
            cGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑Tgg                                                                                                                                       >  2:833672/1‑143 (MQ=255)
            agtggTGGAACCGGTACCTGTTCTGTTTCAATGGGTCGCTGATACCGAGGGAAAGAGTCAGTATGGCCGAGCGCCCGGACTGGCGTGAGACAGCCCGCGAATACGGTTTCAATTTTCGCACCATGTACGGCGAGCCGTAC‑Tgg                                                                                                                                       <  2:373357/139‑1 (MQ=255)
                gTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGGTTCAATTTTCACACCATGTACGGCGAGCCGGAC‑TGGtgtg                                                                                                                                   >  2:1011279/1‑143 (MQ=255)
                gTGCAACCGGTACCTCTTCGCGTTCAATGGGTGGCTGATAGCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGATAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGtgtg                                                                                                                                   <  2:135733/143‑1 (MQ=255)
                  gCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGaa                                                                                                                                 <  2:986707/143‑1 (MQ=255)
                  ccAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCCGGATTACCGGGCGCCGGGACTGGCGTGAGAACGCCCACGAATACGGTTTCAATTTTCCCACCATGTACGGCGAGCCGTAC‑TGGTGTGaa                                                                                                                                 <  2:1054522/142‑1 (MQ=255)
                          tACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTta                                                                                                                         <  2:1043202/143‑1 (MQ=255)
                          tACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACaa                                                                                                                   <  1:738249/149‑1 (MQ=255)
                           cccTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAg                                                                                                                  <  1:771978/148‑1 (MQ=255)
                           acgttctctcGTCCAACGGGGGGCTGCTACCGATGGAAAGAGCCAGAATGACCGGGCGCCCGGACTGGCGTGGGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTtac                                                                                                                        <  2:654952/137‑1 (MQ=255)
                             ctctTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTCCCGAGCGCCCGGACTGGCGTGAGCAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGACGCGTAC‑TGGTGTGAAGGTGATtacac                                                                                                                      >  2:258581/1‑141 (MQ=255)
                                ttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCGTACTACAAGTTGAc                                                                                                             >  1:1264424/1‑149 (MQ=255)
                                ttctCGTGCAGTGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACaa                                                                                                                   <  2:1051286/143‑1 (MQ=255)
                                ttcgCGTTAAATGGGTGGGTGATACCGATGGAAAGAGTCAGGATTACCGAGCGCCCGGCCGGGCGTGAGAGAGCCCACGAATACGGTTCCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACaa                                                                                                                   <  2:131019/143‑1 (MQ=255)
                                    cGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGGTACTTACTACAAGTTg                                                                                                               >  2:778826/1‑143 (MQ=255)
                                    cGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTc                                                                                                         <  1:704100/149‑1 (MQ=255)
                                     gTTCCATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGTGCTCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCg                                                                                                        <  1:91729/149‑1 (MQ=255)
                                     gTTCAATGGGTGGCTGATACCGATGGAAAGGGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCg                                                                                                        <  1:468809/149‑1 (MQ=255)
                                      ttCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGCGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGc                                                                                                       >  1:487037/1‑149 (MQ=255)
                                        cAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGAccc                                                                                                           >  2:432910/1‑143 (MQ=255)
                                          agggttGGCGGATACCGAGGGAAAGAGTCAGTATTACCGCGCGCCCGGACTGGCGGGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTc                                                                                                         <  2:394009/138‑1 (MQ=255)
                                              gTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGATg                                                                                               >  1:1140149/1‑149 (MQ=255)
                                              gTGGCGGAGACCGATGGAAAGAGACAGTATTACCGGGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGccc                                                                                                     <  2:441530/143‑1 (MQ=255)
                                                 gCTGATACCGATGGCAAGAGCCGGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGaaa                                                                                            <  1:352774/149‑1 (MQ=255)
                                                   ggATACCGGGGGAAGGAGTCCGTATTACTGAGCGCCCGGCCGGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGtt                                                                                                <  2:586890/142‑1 (MQ=255)
                                                      tACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTg                                                                                       <  1:118894/149‑1 (MQ=255)
                                                       aCCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCTTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTgg                                                                                      <  1:597145/149‑1 (MQ=255)
                                                       aCCGATGGAAAGAGTCAGTATGACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTgg                                                                                      <  1:672633/149‑1 (MQ=255)
                                                         cGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaa                                                                                    <  1:282053/149‑1 (MQ=255)
                                                          gATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaag                                                                                   <  1:374147/149‑1 (MQ=255)
                                                          gATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTAAAATTTCACACCATGTACGGCGAGACGTAC‑TGGTGTGAAGATGCTTACTACAAGTGGACCCTCGCCCAGGTTGAAAAGc                                                                                         >  2:1083666/1‑143 (MQ=255)
                                                             ggaagggggcAGTATTACCGGGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTgg                                                                                      <  2:410095/134‑1 (MQ=255)
                                                              gAAAGGGTCAGTATCCCCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGga                                                                                     <  2:700662/143‑1 (MQ=255)
                                                               aaaGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGGTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGCGAGAAGCTGCTTACTACAAGGTGACCCCCGCCCAGGTTGAAAACCTGgaa                                                                                    >  2:84040/1‑143 (MQ=255)
                                                                     tCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTAAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgcc                                                                        >  1:1016066/1‑149 (MQ=255)
                                                                       aGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCAACGAATACGGTGTCAATTTTCACACCAGGTACGGCGAGCCGTAC‑GGAGGTGCAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGGAGTCAc                                                                            >  2:1119770/1‑143 (MQ=255)
                                                                          aTTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGGGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgc                                                                         <  2:1118224/143‑1 (MQ=255)
                                                                           ttACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTAACACCATGTACGGCGAGCCGTAC‑TGGTGTAAAGAGGGTTACTAGAAGTCGCCCTTCGACCAGTGTGAAACGCTGAAAAGAGCCAccaaccacccg                                                                  >  1:1144091/1‑140 (MQ=255)
                                                                            tACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGCAGTCACCGCCGAACTGc                                                                 >  1:1159328/1‑149 (MQ=255)
                                                                               cGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAcc                                                              <  1:1278243/149‑1 (MQ=255)
                                                                               cGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCCCCGACCGGCAcc                                                              >  1:876503/1‑149 (MQ=255)
                                                                               cGAGCCCCCGGCCTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAcc                                                              <  1:423729/149‑1 (MQ=255)
                                                                                     ccccGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAcc                                                              >  2:1053464/1‑143 (MQ=255)
                                                                                     cccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTCGACCCTCACCCAGACTGAAAAGCTGGACGAAGCCACCGCCAACCTCCAcc                                                              >  2:397683/1‑143 (MQ=255)
                                                                                     cccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACAACATGTACGGCGAGCCGTCC‑GGGTGTGCAGATGCTTACTACAAGTTGACACTCGCCCAGTTTGAAAAGCTGGAAGAAGTCACAGGCGAACTGCAcc                                                              >  2:119890/1‑143 (MQ=255)
                                                                                        ggACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGcc                                                     >  1:1276904/1‑149 (MQ=255)
                                                                                           cTGGGGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGa                                                  <  1:1071519/149‑1 (MQ=255)
                                                                                                gTGAGAAAGCCCACGAATACGGTTTTAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCTCCGAACTGCACCAGATGTGCCTGAAAgtg                                             <  1:795111/149‑1 (MQ=255)
                                                                                                gTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAGGTGCCTg                                                   >  2:534365/1‑143 (MQ=255)
                                                                                                gTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCGGGAAGAAGTCACCGCCGAACTGCACCAGATGGGCCTg                                                   >  2:605814/1‑143 (MQ=255)
                                                                                                 ggagaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGa                                                  <  2:174064/142‑1 (MQ=255)
                                                                                                  cagaAAGCCCACGGATACGGTTTCAATGTTCACACCATGTACGGCGAGCCGTAC‑GGGTGTGCAGATGCTTACTACAAGTTGACCCTCGCCCAGTTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGaa                                                 <  2:374878/142‑1 (MQ=255)
                                                                                                    gaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGtggtgg                                         <  1:503695/149‑1 (MQ=255)
                                                                                                     aaaGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGCGACGTAC‑AGGTTTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAACAGCTGGAAGAAGTCACCGCCGAACTGCACCAGCTGTGCCTGAAAgt                                              >  2:188430/1‑143 (MQ=255)
                                                                                                     aaaGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACCAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGa                                        >  1:1240162/1‑149 (MQ=255)
                                                                                                     aaaGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGa                                        >  1:470299/1‑149 (MQ=255)
                                                                                                      aaGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGGGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtg                                             >  2:274180/1‑143 (MQ=255)
                                                                                                       aGCCCACGAATACGGTTTCAATTTTCGCACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGCCCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                            <  2:709945/143‑1 (MQ=255)
                                                                                                        gCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaa                                     >  1:945148/1‑149 (MQ=255)
                                                                                                         cccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAATgtggtg                                          >  2:818602/1‑143 (MQ=255)
                                                                                                            aCGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaa                                       >  2:878610/1‑143 (MQ=255)
                                                                                                            aCGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaa                                       >  2:825484/1‑143 (MQ=255)
                                                                                                            aCGAATACGGTTTCAATTTTCACAACATGTATGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaa                                       >  2:129526/1‑143 (MQ=255)
                                                                                                              gAAGACGGTTTCCATTTTCCCACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGAt                               <  1:490927/149‑1 (MQ=255)
                                                                                                               gaTACGGGTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaaa                                    <  2:972656/142‑1 (MQ=255)
                                                                                                               aaTACGGGTTCTATGTTCACAACATGGACGGCGAGCCGGAC‑TGGTGTGAAGATGCGTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATc                              <  1:1195087/149‑1 (MQ=255)
                                                                                                                  aCGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAACGTg                                 >  2:655856/1‑143 (MQ=255)
                                                                                                                   tagtaagcatcttCACACCA‑GTACGGCTCGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAg                         >  1:912889/12‑149 (MQ=38)
                                                                                                                   cGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCTGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGa                                >  2:476792/1‑143 (MQ=255)
                                                                                                                    ggTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCac                         >  1:904546/1‑148 (MQ=255)
                                                                                                                    ggTTTCAATTTCCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAg                         <  1:1128412/149‑1 (MQ=255)
                                                                                                                     gTTTCAATTTTCACCCCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGagc                              >  2:411858/1‑141 (MQ=255)
                                                                                                                      tttCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGCACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCg                             >  2:614777/1‑143 (MQ=255)
                                                                                                                       ttCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGGGGCAAAGGTGacagc                            >  2:452190/1‑139 (MQ=255)
                                                                                                                           aTTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTCCAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAg                  <  1:101287/149‑1 (MQ=255)
                                                                                                                            ttttCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGc                 >  1:24972/1‑149 (MQ=255)
                                                                                                                             cttCACACCA‑GTACGGCTCGCCGTACATGGTGTGAAAATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt                <  1:1182819/148‑1 (MQ=38)
                                                                                                                              ttCACACCA‑GTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGACTGAAAGGGGTGGAAAAAGTGATAGCCAGCGAt                     >  2:796155/1‑143 (MQ=38)
                                                                                                                               tCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTAAAAGTTGGCCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtga              <  1:619809/149‑1 (MQ=255)
                                                                                                                                 acacCATGTACTGGGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAg                  <  2:162227/143‑1 (MQ=255)
                                                                                                                                 acacCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgatg            <  1:140832/149‑1 (MQ=255)
                                                                                                                                  cacCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgatga           >  1:784158/1‑149 (MQ=255)
                                                                                                                                   acCATGTACGGCGAGCAGTAC‑TGGTGTGAAGATGATTACTACAAGTTGACCCTCGCCCAGGTCGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCGGAAAGTGGTGGAAAAAGTGGTCGCCAGCGATGGGCTGATGAc          >  1:573616/1‑149 (MQ=255)
                                                                                                                                     cATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCa        <  1:1019585/149‑1 (MQ=255)
                                                                                                                                     cATGTACGGCGAGCCGTAC‑TGGTGTGAACATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCa        >  1:784290/1‑149 (MQ=255)
                                                                                                                                         tACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAc          >  2:294806/1‑143 (MQ=255)
                                                                                                                                          aCGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTc   <  1:198131/149‑1 (MQ=255)
                                                                                                                                           cGGCGAGCCGTAC‑TGGTGTGAAGATGCTCACGACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGGTGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTcc  <  1:1237793/149‑1 (MQ=255)
                                                                                                                                            ggCGGGCCGTAC‑TGGTGTGAGGGGGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCaa       <  2:456956/143‑1 (MQ=255)
                                                                                                                                            ggCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCac       >  2:1031379/1‑142 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
AGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCC  >  CP009273/3173721‑3174008

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑