Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,377,435 A→T 9.6% intergenic (+32/+39) elaD → / ← yfbK protease, capable of cleaving an AMC‑ubiquitin model substrate/Von Willebrand factor domain putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,377,4350AT 9.6% 203.5 / 8.3 68intergenic (+32/+39)elaD/yfbKprotease, capable of cleaving an AMC‑ubiquitin model substrate/Von Willebrand factor domain putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (27/34);  new base T (4/3);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCT  >  CP009273/2377287‑2377583
                                                                                                                                                    |                                                                                                                                                    
aaCAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACTTCATTGTTAATTAATTCGACCAGTCAGCaa                                                                                                                                                      <  1:650480/149‑1 (MQ=255)
 aCAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCaaa                                                                                                                                                     >  2:164425/1‑149 (MQ=255)
   aGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAata                                                                                                                                                   <  1:37775/149‑1 (MQ=255)
   aGAAATTATTATTTAGATAAATACGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAata                                                                                                                                                   >  2:845146/1‑149 (MQ=255)
          gTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGAc                                                                                                                                            <  1:986549/149‑1 (MQ=255)
             aTTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTgg                                                                                                                                         <  1:1026483/149‑1 (MQ=255)
              tttAGAAAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGt                                                                                                                                        <  2:454074/149‑1 (MQ=255)
               ttAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTc                                                                                                                                       <  1:45240/149‑1 (MQ=255)
               ttAGATAAAAATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAGAACACTCACAACCCAAAAATAGAGCATCCAGCAAAAGAGTGAGTTCAAGTCTTTGTTAATTAATCCGACCAGTCAGCAAATATGCTGACTGGTc                                                                                                                                       <  2:952825/149‑1 (MQ=255)
                   aTAAATATGGCGATGCAAATATCAACCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAAc                                                                                                                                   <  2:688317/149‑1 (MQ=255)
                      aaTATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATa                                                                                                                                <  1:179580/149‑1 (MQ=255)
                      aaTATGGCGATGCAAATAACAATCATTATTATAGAAACCTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGCGTGAGTTAACATCATTGTTAATTAATCCGACCAGTCAGCAAATATGCTGACTGGTCGAACATa                                                                                                                                <  1:704174/149‑1 (MQ=255)
                         atGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCTCCCTATCACGTTTGGTCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGtt                                                                                                                             <  1:596081/149‑1 (MQ=255)
                            gCGATTCAAATATCACTCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAACTAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATTCTGACTGGTCGAACATAGTTTGc                                                                                                                          <  2:392752/149‑1 (MQ=255)
                                  cAAATATCAATCATTATGATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCAATGTTAATTAATTCGACCAGTCAGCTAATATGCTGACTGGTCGAACATAGTTTGCTAAACa                                                                                                                    <  2:880600/149‑1 (MQ=255)
                                   aaaTATCAATCATTATTATAGAAACTTAGAAAGAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAg                                                                                                                   <  2:465083/149‑1 (MQ=255)
                                    aaTATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGtc                                                                                                                 <  1:1170374/149‑1 (MQ=255)
                                          aatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCac                                                                                                            >  1:1042900/1‑149 (MQ=255)
                                          aatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCac                                                                                                            >  2:217663/1‑149 (MQ=255)
                                          aatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTAATCac                                                                                                            >  2:798611/1‑149 (MQ=255)
                                          aatcaTTATTATAGAAACTTAGAAACAAAACACTCACCACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCac                                                                                                            <  1:769668/149‑1 (MQ=255)
                                             caTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactga                                                                                                         <  1:217663/149‑1 (MQ=255)
                                                tattatAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgactg                                                                                                      >  1:1079421/1‑149 (MQ=255)
                                                     taGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGGTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGTTAAACAGTCATCactgacctgtgtc                                                                                                 >  1:710225/1‑142 (MQ=255)
                                                        agaCTTAGAAATAAAACACTCACAACCCAAACATAGAGCATCGGGCAAAAGAGTGTGGTAATAACGTTGTTAATTAAGTCGACCGGTCCGAATATATGCTGACTGGTCGATCTTTGTTTGCTAAACAGTCATCACTGACTGTTgtcagt                                                                                              <  2:317097/147‑1 (MQ=255)
                                                        aaaCTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtcagt                                                                                              <  1:792122/149‑1 (MQ=255)
                                                        aaaCTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGATTGATAAACAATCATCACTGACTGATgtcaat                                                                                              >  1:739281/1‑147 (MQ=255)
                                                          aCTCAGAAAGAAAACACACACAACCAAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtcagtca                                                                                            <  2:84417/149‑1 (MQ=255)
                                                            ttAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTCACCAGTCGTCACTGACTGATGTCAGTcaca                                                                                          <  2:578214/149‑1 (MQ=255)
                                                            ctAGAAATAAAACACTCACTACACAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTAGACCAGTCAACAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTcaca                                                                                          <  1:962718/148‑1 (MQ=255)
                                                              aGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACAcc                                                                                        <  2:899215/149‑1 (MQ=255)
                                                                   tAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGGCAGTCCAACCATCCg                                                                                   >  2:3750/1‑149 (MQ=255)
                                                                     aaaCACTCACAACCCAAAAATAGAGCATCCGGCAAAAAAGTGAGTTAACATCATAGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAGCAGTCATCACTAACGGGTGTCAAGCACGCCATCCGcc                                                                                 >  2:1158190/1‑149 (MQ=255)
                                                                      aaCACTCACAACCCAAAAATAGAGCATCCGCAAAAAGGGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCa                                                                                <  1:538566/149‑1 (MQ=255)
                                                                           tCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCa                                                                           >  1:388701/1‑149 (MQ=255)
                                                                             acaACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCTGTCACACCATCCGCCAGTTCAAt                                                                         >  2:257376/1‑149 (MQ=255)
                                                                             acaACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAAt                                                                         >  1:399785/1‑149 (MQ=255)
                                                                             acaACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGACAGTACAAt                                                                         >  2:703358/1‑149 (MQ=255)
                                                                             acaACCCAAAAATAGAGCAGCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAAt                                                                         >  2:839405/1‑149 (MQ=255)
                                                                               aaCCCAAAAATAGGGTAGCCGGCAAAAGAGTGAGGTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCa                                                                       <  1:108939/149‑1 (MQ=255)
                                                                                  ccAAAAAGAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGc                                                                    <  1:1086856/149‑1 (MQ=255)
                                                                                     aaaaTAGCGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGaa                                                                 <  1:422108/149‑1 (MQ=255)
                                                                                       aaTAGACCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATa                                                               <  2:1142176/149‑1 (MQ=255)
                                                                                          agagCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACACTCATCACGGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAAt                                                            >  1:916513/1‑149 (MQ=255)
                                                                                              cATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCg                                                        >  1:1128247/1‑149 (MQ=255)
                                                                                              cATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCg                                                        >  1:992158/1‑149 (MQ=255)
                                                                                              cATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGCCAGTCACACCATCCGCCAGACCAATCAGGCGAATAAATTCCg                                                        >  2:3018/1‑149 (MQ=255)
                                                                                               aTCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGc                                                       >  2:499714/1‑149 (MQ=255)
                                                                                                   ggCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCCGGCGCATAAATTCCGCCCTg                                                   >  1:319750/1‑149 (MQ=255)
                                                                                                   ggCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTg                                                   >  2:1109072/1‑149 (MQ=255)
                                                                                                             tgagtTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCAATAACGGGTTGCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCtgtg                                         >  1:1215025/1‑149 (MQ=255)
                                                                                                                gtTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGAt                                      >  1:1199508/1‑149 (MQ=255)
                                                                                                                  taattttatTGTTAACTAATTTGACCAGTCAGCATAGTTGCTGACTTGTCGATCATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCt                                    <  2:853689/142‑1 (MQ=255)
                                                                                                                  tAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCt                                    <  2:506208/149‑1 (MQ=255)
                                                                                                                    aCATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCGGGCTAATAAATTCCGCCCTGTAACCCTGTGGATCTTc                                  <  1:788232/149‑1 (MQ=255)
                                                                                                                    aCATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTc                                  >  1:70354/1‑149 (MQ=255)
                                                                                                                             ttaattaaTTCGACCAGTCAGCAAATATGCTGACTGGTCGAACCTAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGc                         <  1:164425/149‑1 (MQ=255)
                                                                                                                                attaattCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCctg                      <  1:798611/149‑1 (MQ=255)
                                                                                                                                attaattCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCctg                      >  2:600891/1‑149 (MQ=255)
                                                                                                                                 ttaattCGACCAGTCAGCAAATATGCTGACTGCTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCctgc                     <  1:845146/149‑1 (MQ=255)
                                                                                                                                       cGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGcc               <  2:992158/149‑1 (MQ=255)
                                                                                                                                        gACCAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGccc              >  1:1229674/1‑149 (MQ=255)
                                                                                                                                        gACCAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGccc              >  1:129299/1‑149 (MQ=255)
                                                                                                                                         aCCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACGGTACCCATTGTCGGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTTCCTGCTGAGCCCa             <  1:847445/149‑1 (MQ=255)
                                                                                                                                            aGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCAcgg          >  2:537121/1‑147 (MQ=255)
                                                                                                                                             gTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGt         <  2:390534/149‑1 (MQ=255)
                                                                                                                                                aGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGCGCCCCCTGttgg      >  2:176627/1‑147 (MQ=255)
                                                                                                                                                aGCAAATATGCTGACTGGTCGAACACAGTTCGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTg      >  1:642870/1‑149 (MQ=255)
                                                                                                                                                 gCAAATATGCTGACTGGTCGAACCTGGTTTGCTATACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGa     <  2:777076/149‑1 (MQ=255)
                                                                                                                                                 gCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCAAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGa     >  1:837509/1‑149 (MQ=255)
                                                                                                                                                   aaaTATGCTGACTGGTCGAACATAGTTTGCTAAACCGTCATCACTGCCTGATGTCAGTCACACCCTCCGCCAGTTCAATCAGGCGAACAAACTCCGCCCTGTACCCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATc   <  2:604409/149‑1 (MQ=255)
                                                                                                                                                    aaTATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATct  <  2:388701/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCT  >  CP009273/2377287‑2377583

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