Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,377,438 A→T 8.4% intergenic (+35/+36) 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,4380AT 8.4% 196.0 / 7.2 68intergenic (+35/+36)elaD/yfbKprotease, capable of cleaving an AMC‑ubiquitin model substrate/Von Willebrand factor domain putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (26/36);  new base T (4/2);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

AGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGC  >  CP009273/2377290‑2377585
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                              >  2:217663/1‑149 (MQ=255)
                                       aatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCac                                                                                                              >  1:1042900/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                        >  2:600891/1‑149 (MQ=255)
                                                                                                                             attaattCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCctg                        <  1:798611/149‑1 (MQ=255)
                                                                                                                              ttaattCGACCAGTCAGCAAATATGCTGACTGCTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCctgc                       <  1:845146/149‑1 (MQ=255)
                                                                                                                                    cGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGcc                 <  2:992158/149‑1 (MQ=255)
                                                                                                                                     gACCAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGccc                >  1:129299/1‑149 (MQ=255)
                                                                                                                                     gACCAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGccc                >  1:1229674/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)
                                                                                                                                                   tatGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCGGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATctgc  >  2:687407/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGC  >  CP009273/2377290‑2377585

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: