Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,679,847 G→A G111S (GGT→AGT)  cbpD → possible secretory antigen

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,679,8470GA98.6% 209.1 / NA 70G111S (GGT→AGT) cbpDpossible secretory antigen
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/38);  minor base C (0/1);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTAACTATATTTACTAATCAAAAAACTTCTATCACGAACGCTTCAA  >  CP000730/2679718‑2679977
                                                                                                                                  |                                                                                                                                  
aaaaaTTGGTATTACAATGAACGAAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACagtagt                                                                                                                               >  1:491221/1‑135 (MQ=255)
 aaaaTTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCTACAACCTATATACAAGTCGCAGTTGTTGAACAAGTATATAAAAACAgtagta                                                                                                                              >  1:123540/1‑135 (MQ=255)
 aaaaTGTGTATCACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAgtagta                                                                                                                              <  2:625077/135‑1 (MQ=255)
 aaaaATGGTATTACAATGCACGTAAATCAAAATATTCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAgtagta                                                                                                                              <  2:622851/135‑1 (MQ=255)
     ttGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAAAAGTAGTATTaa                                                                                                                          >  1:227389/1‑135 (MQ=255)
                 tGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATa                                                                                                              <  1:636866/135‑1 (MQ=255)
                  gCACGTAAATCAAAATATGCGACTGGTCGTACTCCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATAc                                                                                                             <  1:407663/135‑1 (MQ=255)
                  gCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATAc                                                                                                             <  2:190195/135‑1 (MQ=255)
                    aCGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACaa                                                                                                           >  2:552834/1‑135 (MQ=255)
                    aCGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTGGAACAAGTATATAAAAACAGTAGTATTAAAGTTGCAAAATACaa                                                                                                           >  2:542127/1‑135 (MQ=255)
                    aCGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTACGGACATGTCCCAGTTGTTGAACAAGTAAAAAAAAACAGAAGTCTTAAGGTTTCAGCATACaa                                                                                                           >  2:39706/1‑135 (MQ=255)
                     cGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCCTGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTAAGAATACAAc                                                                                                          >  1:400032/1‑135 (MQ=255)
                     cGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGAATACAAAAACAGTAGTATTAAAGTTTCAGAATACAAc                                                                                                          >  2:534798/1‑135 (MQ=255)
                      gTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACt                                                                                                         >  1:454929/1‑135 (MQ=255)
                        aaaTCAAAATCTGCGACTGTTCATACACCAAAAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAAAttt                                                                                                       <  1:163131/135‑1 (MQ=255)
                           tgaaaaTATGCTACAGCTCGTACAGTAAGAAAATTCGCGGTCATGCAATCA‑ACTGCAGGATATTACGGACCTGTCGCAGTTGTTCAATAAGTATATTAAAACCGTAGTTTTAAATTTTCAGAATACAACTTTTAt                                                                                                    <  2:527867/133‑1 (MQ=255)
                           tCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTAt                                                                                                    >  1:619383/1‑135 (MQ=255)
                           tCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTAt                                                                                                    <  2:192070/135‑1 (MQ=255)
                           tCAAAAGATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTAt                                                                                                    <  2:527591/135‑1 (MQ=255)
                            cAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGACCAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATc                                                                                                   <  1:394983/135‑1 (MQ=255)
                            cAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATG‑ACTAACACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATc                                                                                                   <  2:636349/135‑1 (MQ=255)
                            cAAAAGATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATc                                                                                                   <  2:166511/135‑1 (MQ=255)
                              aaaaATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTCTTAAAGTTTCAAAATACAACCTTTAtacc                                                                                                 >  2:533046/1‑132 (MQ=255)
                                atatGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAAACA‑AATGCAGGATATTATGGACATGTAGCAATTGTTGAACAATCATATAAAAACAATAGTATTAAAGTTTCAGAATACAACTTTTATCGccc                                                                                               >  2:383727/1‑135 (MQ=255)
                                atatGCGACTGGTCGTACAACAAGAAAATTAGCTGTAATACAAACA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAATATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCTccc                                                                                               >  1:185867/1‑135 (MQ=255)
                                  atGCGACTGGTGGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCAt                                                                                             >  1:289116/1‑135 (MQ=255)
                                          ggTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCTCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACCACTTTTATCGCCCATTAAaatac                                                                                     <  1:374842/135‑1 (MQ=255)
                                              gTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAaatacaata                                                                                 <  2:427918/135‑1 (MQ=255)
                                                acacCAAGAAAACTCTCTGTCATGCAATCA‑ACTGCAGGATAATATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACATTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAatacaataca                                                                               <  1:175919/135‑1 (MQ=255)
                                                  acCAATAAAACTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACg                                                                             <  2:169214/135‑1 (MQ=255)
                                                  acCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTACAATACAATACACg                                                                             >  1:526996/1‑135 (MQ=255)
                                                  acCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACg                                                                             >  1:483803/1‑135 (MQ=255)
                                                   ccAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACgt                                                                            <  2:6663/135‑1 (MQ=255)
                                                   ccAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACgt                                                                            >  1:373555/1‑135 (MQ=255)
                                                   ccAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGGATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATAAACgt                                                                            >  2:536680/1‑135 (MQ=255)
                                                   ccAAGAAAATTTGCTGTCATACAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACgt                                                                            >  1:292755/1‑135 (MQ=255)
                                                      aggacaTTGGCTGCAATGCAATTA‑TGTGCAGGATATTGTGGACATGTCGCTGTTGTTGAACAAGTATATAAAAACAGTAGTATTTAAGTTTCAGAATACAACTTTCATCGCCCATTAAAATACAATTCGCGTGTa                                                                         <  2:74511/130‑1 (MQ=255)
                                                      aGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTa                                                                         <  2:639777/135‑1 (MQ=255)
                                                         aaaTTTGCTGTCATGCAATCA‑ACTGCAGAATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTa                                                                      >  1:581347/1‑135 (MQ=255)
                                                            tttGCTGTCATGCAATCA‑ACTACAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTCGTATTAAAGTTTCAGAATACAACTTTTATGGCCAATTAAAATACAATAAACGTGTACTaacc                                                                   >  1:136617/1‑133 (MQ=255)
                                                            ttggCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATTGTCCATTAATATATTATACTCGTGTAATAAGc                                                                   <  1:8420/132‑1 (MQ=255)
                                                            attGCTGTCATGATTACA‑ACTTCAGGCTATTATGGTCATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAGAGTTTCATAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGc                                                                   >  1:462197/2‑135 (MQ=255)
                                                             ttGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAGACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCa                                                                  <  1:534798/135‑1 (MQ=255)
                                                             ttGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCa                                                                  <  1:382243/135‑1 (MQ=255)
                                                                ttGTCACGTAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATATCCGTGCAGTAAGCaaaa                                                               <  1:624504/134‑1 (MQ=255)
                                                                cTGTCCTGCAATCA‑ACTGCAGGATATTATGGACATGTAGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCGTTAAAATACAATACACGTGTACTAAGAaaaa                                                               >  1:7602/1‑135 (MQ=255)
                                                                   tCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAgg                                                            >  1:675110/1‑135 (MQ=255)
                                                                      tGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACTCGTGTACTAAGCAAAAAGgcag                                                         <  1:553999/135‑1 (MQ=255)
                                                                       gtaaTCA‑ACGGCATCATATGAGGGAAATGTCGCAGTTGTGGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGgcagc                                                        <  2:123540/133‑1 (MQ=255)
                                                                       gCAATCA‑ACTGCAGGATATTATGGACATGTCCCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTATTAAGCAAAAAGgcagc                                                        <  1:291382/135‑1 (MQ=255)
                                                                        cAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGgtagca                                                       <  2:373374/135‑1 (MQ=255)
                                                                            cA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTa                                                   <  2:573027/135‑1 (MQ=255)
                                                                                 ggtaGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACtt                                               <  2:137284/132‑1 (MQ=255)
                                                                                    aGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTaa                                            >  2:511580/1‑135 (MQ=255)
                                                                                    aGGATATTATGGACATGTCGAAAATAAATAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTaa                                            >  1:668597/1‑135 (MQ=255)
                                                                                     ggATATTATGGAACTGTCGCAGTTGTTGAACAAGCAGATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTATGTAAAATGGTAGCACGTAACTTTAAc                                           <  1:291364/135‑1 (MQ=255)
                                                                                      ggtattatGGTCTTGTCGCAGTTGTTGTACAACTATATAGAACCAGTAGTATTATAGTTTCAGGATACAACTTTTATTGGCTGTTATAATTCTTTGCACGAATACAAAGCAAAAAGGCAGCACGTAACTTTAACt                                          <  2:270351/133‑1 (MQ=255)
                                                                                       atatTATGGACATGTCACAGTTATAGAACAAGTATATAAAAACAGTAGTATTAAAATTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAACGGCAGCACGTAACTTTAAcga                                         >  1:48606/1‑133 (MQ=255)
                                                                                        tattatGTACATGTCGCAGTTGTTTAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTAACtat                                        <  2:377890/135‑1 (MQ=255)
                                                                                                    gTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCCTTAAAATACAATACACGTGTACTAAGAAAAAAGGCAGCACGTAACTTTAACTATATTTACTAATCa                            >  2:293200/1‑135 (MQ=255)
                                                                                                    gTCGCAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTAACTATATTTACTACTCa                            >  2:46868/1‑135 (MQ=255)
                                                                                                    gTCGCAGTTGTTGAACAAGAATATAAAAACAGTAGTATTAAAGTTTCAGAATAAAACTTTTATAGACAATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTAACTATATTTACTAATCa                            >  1:300817/1‑135 (MQ=255)
                                                                                                    cgctcAGTTGTTGAACAAGTATATAAAAACAGTACTATTAAAGTTTCATAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCACCACGTAACTTTAACTATATTTACTAATCa                            <  2:468416/131‑1 (MQ=255)
                                                                                                        cAGTTGTTGAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAATTTTAACTATATTTCCTAATCaaaaa                        <  1:282760/135‑1 (MQ=255)
                                                                                                          gttgttGAACAAGTATATAAAAACAGTAGTCTTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCCAAAAGGCAGCACGTCACTTTACCTATATTTACTAATAAaacaac                      >  2:72203/1‑135 (MQ=255)
                                                                                                            tgttgAACAAGTATATAAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTAACTATATTTACTAATCAAAAAACtt                    >  1:383689/1‑135 (MQ=255)
                                                                                                                       tatataAAAACAGTAGTATTAATGTTTCAGCATACAACTTTTATTGCCCATTAAAATACAATACACGTGTACTAAGCAATCAGGCAGCACGTAACTTTAACTATATTTTCTATTCTAAACACTTCTATCACGAAc         <  2:38326/135‑1 (MQ=255)
                                                                                                                       tatagaAAAACAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTAACTATATTTACTAATCAAAAAACTTCTATCACGAAc         <  2:44985/135‑1 (MQ=255)
                                                                                                                              aaaCAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTAACTATATTTACTAATTAAAAATCTTTTATCACGAACGCTTCaa  <  1:551202/135‑1 (MQ=255)
                                                                                                                              aaaCAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTAACTATATTTACTAATCAAAAAACTTCTATCACGAACGCTTCaa  <  2:454929/135‑1 (MQ=255)
                                                                                                                              aaaCAGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTAACTATATTTACTAATCAAAAAACTTCTATCACGAACGCTTCaa  <  1:560616/135‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                  
AAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCA‑ACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTAACTATATTTACTAATCAAAAAACTTCTATCACGAACGCTTCAA  >  CP000730/2679718‑2679977

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