Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,173,856 Δ1 bp 5.3% coding (77/1161 nt) ygiC → ATP‑Grasp family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,173,8560T. 5.3% 281.3 / 9.9 94coding (77/1161 nt)ygiCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base T (44/44);  new base . (2/3);  total (47/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

TGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCA  >  CP009273/3173708‑3174002
                                                                                                                                                    |                                                                                                                                                   
tGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAACGGGCGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCCCGAATACGGTTTCAATTTTGACACCAt                                                                                                                                                     >  1:186540/1‑149 (MQ=255)
     ccAAACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACg                                                                                                                                                <  1:1096392/149‑1 (MQ=255)
        aacaaaGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGGGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGt                                                                                                                                                   >  2:712053/1‑143 (MQ=255)
           aaaGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGGGAGAAAGCCCACGAATACGTGTTGAAGTTTCACACCATGGACg                                                                                                                                                >  2:365463/1‑143 (MQ=255)
             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:825484/1‑143 (MQ=255)
                                                                                                                         aCGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaa                                 >  2:878610/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)
                                                                                                                                                    |                                                                                                                                                   
TGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCA  >  CP009273/3173708‑3174002

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