Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,679,803 G→T G96V (GGA→GTA)  cbpD → possible secretory antigen

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,679,8030GT98.7% 258.4 / ‑3.0 76G96V (GGA→GTA) cbpDpossible secretory antigen
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (33/42);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

AAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTAAC  >  CP000730/2679667‑2679936
                                                                                                                                        |                                                                                                                                     
aaaCGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTatat                                                                                                                                   >  1:706119/1‑141 (MQ=255)
 aaCGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTAtatt                                                                                                                                  >  1:437077/1‑141 (MQ=255)
 aaCGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTAtatt                                                                                                                                  >  2:824430/1‑141 (MQ=255)
   cGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTAtattat                                                                                                                                <  2:408013/141‑1 (MQ=255)
    gTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATg                                                                                                                               <  2:114175/141‑1 (MQ=255)
         gcACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACAt                                                                                                                          >  2:639467/1‑141 (MQ=255)
          cacaAGTCGGTAAACCTGTGTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACAGCAGTATATTATGGACATg                                                                                                                         >  2:111813/1‑141 (MQ=255)
            caAGTCGGTAAACCTGTTTCTAAGAGATGGGGAAATGCTAAAAAGTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAAGCAACTGCAGTATATTATGGACATGTc                                                                                                                       <  1:618383/141‑1 (MQ=255)
                   gTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAgttg                                                                                                                >  2:270195/1‑141 (MQ=255)
                    tAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAgttgt                                                                                                               <  1:864520/141‑1 (MQ=255)
                      aaCCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGttgttg                                                                                                             >  1:414626/1‑141 (MQ=255)
                        ccTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGaa                                                                                                           <  2:437077/141‑1 (MQ=255)
                        ccTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGaa                                                                                                           >  1:979596/1‑141 (MQ=255)
                        ccTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGaa                                                                                                           >  1:513421/1‑141 (MQ=255)
                         cTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAAc                                                                                                          <  1:431968/141‑1 (MQ=255)
                             ttCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGt                                                                                                      <  1:505954/141‑1 (MQ=255)
                                  aTAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGtatata                                                                                                 <  1:564812/141‑1 (MQ=255)
                                      aTGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATaaaaa                                                                                             >  1:703856/1‑141 (MQ=255)
                                        ggggCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACg                                                                                           >  1:371930/1‑141 (MQ=255)
                                         gggCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACgg                                                                                          <  2:399919/141‑1 (MQ=255)
                                         gggCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACgg                                                                                          <  2:399946/141‑1 (MQ=255)
                                          ggCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGgt                                                                                         >  1:985870/1‑141 (MQ=255)
                                              aTGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGgtagta                                                                                     >  2:983723/1‑141 (MQ=255)
                                                gCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTAtt                                                                                   <  1:54933/141‑1 (MQ=255)
                                                gCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTAtt                                                                                   <  2:74410/141‑1 (MQ=255)
                                                gCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTAtt                                                                                   <  2:717757/141‑1 (MQ=255)
                                                gCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTAtt                                                                                   <  1:621816/141‑1 (MQ=255)
                                                 cTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTa                                                                                  <  1:701962/141‑1 (MQ=255)
                                                 cTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTa                                                                                  <  1:612080/141‑1 (MQ=255)
                                                   aaaaaTTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAgttgt                                                                                                               <  2:889880/110‑1 (MQ=255)
                                                   aaaaaTTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAgttgt                                                                                                               >  1:889880/1‑110 (MQ=255)
                                                    aaaaTTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGGCATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAg                                                                               <  2:931037/141‑1 (MQ=255)
                                                     aaaTTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGt                                                                              >  1:707919/1‑141 (MQ=255)
                                                     aaaTTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGt                                                                              <  1:314493/141‑1 (MQ=255)
                                                                 cAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACaa                                                                  <  1:474351/141‑1 (MQ=255)
                                                                 cAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACaa                                                                  >  1:466241/1‑141 (MQ=255)
                                                                  aaTGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAAc                                                                 <  1:788178/141‑1 (MQ=255)
                                                                    tGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGtatat                                                                                                  <  2:273750/106‑1 (MQ=255)
                                                                    tGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGtatat                                                                                                  >  1:273750/1‑106 (MQ=255)
                                                                       aCGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTa                                                            >  2:49851/1‑141 (MQ=255)
                                                                          tAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCg                                                         <  1:767660/141‑1 (MQ=255)
                                                                          tAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCg                                                         <  2:643269/141‑1 (MQ=255)
                                                                           aaaTCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGc                                                        <  2:979596/141‑1 (MQ=255)
                                                                           aaaTCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGc                                                        <  2:371930/141‑1 (MQ=255)
                                                                           aaaTCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCGGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGc                                                        >  1:803453/1‑141 (MQ=255)
                                                                            aaTCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGcc                                                       >  2:962066/1‑141 (MQ=255)
                                                                             aTCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGccc                                                      >  1:513533/1‑141 (MQ=255)
                                                                             aTCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGccc                                                      <  2:243476/141‑1 (MQ=255)
                                                                              tCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCa                                                     <  1:177566/141‑1 (MQ=255)
                                                                               cAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCAt                                                    <  1:991745/141‑1 (MQ=255)
                                                                               cAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCAt                                                    <  1:914577/141‑1 (MQ=255)
                                                                                aaaaTATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCAtt                                                   >  2:177425/1‑141 (MQ=255)
                                                                                  aaTATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTaa                                                 <  2:92515/141‑1 (MQ=255)
                                                                                       gcgcCTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAaatac                                            <  1:445030/141‑1 (MQ=255)
                                                                                       gCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAaatac                                            <  2:781435/141‑1 (MQ=255)
                                                                                       gCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAaatac                                            <  2:941380/141‑1 (MQ=255)
                                                                                                     acCAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACt                              >  1:316220/1‑141 (MQ=255)
                                                                                                       cAAGAAAATTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAGTACAATACACGTGTACTaa                            >  2:433845/1‑141 (MQ=255)
                                                                                                              aTTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAGGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCaaaaa                     >  1:471061/1‑141 (MQ=255)
                                                                                                              aTTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGcc                                                       >  2:13385/1‑107 (MQ=255)
                                                                                                              aTTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGcc                                                       <  1:13385/107‑1 (MQ=255)
                                                                                                              aTTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCaaaaa                     >  2:905839/1‑141 (MQ=255)
                                                                                                              aTTTGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGAACTAAGCaaaaa                     >  2:184750/1‑141 (MQ=255)
                                                                                                                ttGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAgg                   >  2:189426/1‑141 (MQ=255)
                                                                                                                 tGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGgc                  >  1:1025349/1‑141 (MQ=255)
                                                                                                                 tGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGgc                  >  2:942538/1‑141 (MQ=255)
                                                                                                                 tGCTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGgc                  >  1:52287/1‑141 (MQ=255)
                                                                                                                   cTGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACTGTAGTATTAAAGTTTCAGAATACAACTTTTGTCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGgcag                >  1:617088/1‑141 (MQ=255)
                                                                                                                    tGTCATGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGgcagc               <  1:527737/141‑1 (MQ=255)
                                                                                                                        aTGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGt           <  1:321783/141‑1 (MQ=255)
                                                                                                                        aTGCAATCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGt           <  2:755538/141‑1 (MQ=255)
                                                                                                                            aaTCAACNGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACt       <  2:771557/141‑1 (MQ=255)
                                                                                                                              tCAACTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACttt     <  1:707136/141‑1 (MQ=255)
                                                                                                                                 aCTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTAAc  >  1:567348/1‑141 (MQ=255)
                                                                                                                                 aCTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTAAc  <  1:933509/141‑1 (MQ=255)
                                                                                                                                 aCTGCAGTATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAACAAAAAGGCAGCACGTAACTTTAAc  <  2:129834/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
AAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTTTCAGAATACAACTTTTATCGCCCATTAAAATACAATACACGTGTACTAAGCAAAAAGGCAGCACGTAACTTTAAC  >  CP000730/2679667‑2679936

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