Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,272,132 G→C 5.0% intergenic (+28/‑38) clpQ → / → clpY T01 family HslV component of HsIUV peptidase/T01 family HslU component of HsIUV peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,272,1320GC 5.0% 245.4 / 3.5 75intergenic (+28/‑38)clpQ/clpYT01 family HslV component of HsIUV peptidase/T01 family HslU component of HsIUV peptidase
Reads supporting (aligned to +/- strand):  ref base G (37/33);  new base C (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AACGCCATGCATCGCATTTGTCTGCTGAAGAAATGGCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAATGATGCTAAACGTAAAGTGGCAAT  >  CP000730/1271995‑1272267
                                                                                                                                         |                                                                                                                                       
aaCGCCATGCATCGCATTTGTCTGCTGAAGAAATGGCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCaa                                                                                                                                      <  1:909454/141‑1 (MQ=255)
        gCATCGCATTTGTCTGCTGAAGAAATGGCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAAtttt                                                                                                                              >  1:722901/1‑141 (MQ=255)
         cATCGCATTTGTCTGCTGAAGAAATGGCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTTGAAACACTATAATAATCAGAGCACGATAAATAATTACGACCAATTAATTTTa                                                                                                                             <  1:309334/141‑1 (MQ=255)
         cATCGCATTTGTCTGCTGAAGAAATGGCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTa                                                                                                                             <  1:459463/141‑1 (MQ=255)
              cATTTGTCTGCTGAAGAAATGGCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTaa                                                                                                                        >  2:838529/1‑141 (MQ=255)
                           aagaaATGGCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAAt                                                                                                           <  1:252696/141‑1 (MQ=255)
                           aagaaATGGCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCCTTACCAACGATAATATTGTTGTCGAAACACTATAATACTCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAAt                                                                                                           <  1:222606/141‑1 (MQ=255)
                            agaaATGGCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATg                                                                                                          <  1:420421/141‑1 (MQ=255)
                              aaaTGGCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGat                                                                                                        <  1:1031418/141‑2 (MQ=255)
                                aTGGCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAAt                                                                                                      >  1:862117/1‑141 (MQ=255)
                                  ggCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTa                                                                                                    >  1:150180/1‑141 (MQ=255)
                                  ggCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTa                                                                                                    >  1:75318/1‑141 (MQ=255)
                                   gCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTaa                                                                                                   >  1:315089/1‑141 (MQ=255)
                                     atatGAGAGCTTGAAAGTAGCGGCTTATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATg                                                                                                 <  1:670354/141‑1 (MQ=255)
                                       atGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGa                                                                                               <  2:476670/141‑1 (MQ=255)
                                       atGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGa                                                                                               <  1:986112/141‑1 (MQ=255)
                                            agCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAg                                                                                          <  2:131220/141‑1 (MQ=255)
                                             gCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGc                                                                                         >  1:1076423/1‑141 (MQ=255)
                                              cTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCt                                                                                        <  2:357886/141‑1 (MQ=255)
                                                tGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTgg                                                                                      <  1:971818/141‑1 (MQ=255)
                                                tGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTgg                                                                                      <  1:994110/141‑1 (MQ=255)
                                                   aaGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTTGAAACACTATAATAATCAGAGCACGATAAATAATTACGACCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAAt                                                                                   >  1:475761/1‑141 (MQ=255)
                                                   aaGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAAt                                                                                   >  2:602573/1‑141 (MQ=255)
                                                    aGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATa                                                                                  >  1:345210/1‑141 (MQ=255)
                                                         cGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGAtt                                                                             >  1:23549/1‑141 (MQ=255)
                                                                  tttGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCaaa                                                                    <  1:625168/141‑1 (MQ=255)
                                                                    tgtgTCTTGACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATCAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAg                                                                  <  1:942492/141‑1 (MQ=255)
                                                                     gtgtCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGTTACAGCTGGAATAAGATTAACTCCAAAAGa                                                                 <  2:421774/141‑1 (MQ=255)
                                                                     gtgtCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGa                                                                 >  2:296057/1‑141 (MQ=255)
                                                                     gtgtCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGa                                                                 >  2:381594/1‑141 (MQ=255)
                                                                     gtgtCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGa                                                                 >  2:624112/1‑141 (MQ=255)
                                                                       gtCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGaaa                                                               >  1:580570/1‑141 (MQ=255)
                                                                         cTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATc                                                             <  2:163130/141‑1 (MQ=255)
                                                                             aCCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTAt                                                         >  2:612340/1‑141 (MQ=255)
                                                                              ccAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATc                                                        <  2:990919/141‑1 (MQ=255)
                                                                              ccAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATc                                                        <  1:642419/141‑1 (MQ=255)
                                                                              ccAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATc                                                        <  2:826526/141‑1 (MQ=255)
                                                                               catcgagcgtatTGCTGTCGAAACGCCATAATAATGAGCGCAAGATTGAGAACTACGAGCAATGAATTTTAGCTACCAGACGGAGGAATGACCTTAATGGATACAGCTGGAATAGGATTATCGCCAAAGGAACGCGTACCt                                                       <  2:374658/132‑1 (MQ=255)
                                                                                   gATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCtaaat                                                   >  2:717235/1‑141 (MQ=255)
                                                                                   gATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAAGTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCACATGTGTTCGTATCtaaat                                                   <  1:1071897/141‑1 (MQ=255)
                                                                                   gATAATATTGTTGTCGAAACACTATAATAACCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCtaaat                                                   <  1:585703/141‑1 (MQ=255)
                                                                                    ataataTTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCtaaatt                                                  >  1:1274/1‑141 (MQ=255)
                                                                                      aataTTGTTGGCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCtaaattaa                                                >  2:768156/1‑141 (MQ=255)
                                                                                        tatTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCtaaattaaat                                              >  2:470764/1‑141 (MQ=255)
                                                                                            gttgtCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAAt                                          >  1:54183/1‑141 (MQ=255)
                                                                                               gtCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACa                                       >  2:740252/1‑141 (MQ=255)
                                                                                                   aaaCACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGt                                   >  2:41449/1‑141 (MQ=255)
                                                                                                    aaCACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGtt                                  >  1:856317/1‑141 (MQ=255)
                                                                                                    aaCACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGtt                                  >  1:345143/1‑141 (MQ=255)
                                                                                                     acacTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTg                                 <  1:78209/141‑1 (MQ=255)
                                                                                                      cacTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTgg                                >  2:276430/1‑141 (MQ=255)
                                                                                                      cacTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTgg                                <  1:921456/141‑1 (MQ=255)
                                                                                                      cacTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAACGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTgg                                >  2:415418/1‑141 (MQ=255)
                                                                                                      cacTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTGAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTgg                                <  2:801277/141‑1 (MQ=255)
                                                                                                         tataATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACa                             >  1:974094/1‑141 (MQ=255)
                                                                                                         tataATAATCAGAGCACGATAAATAATTACGACCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATATATCGTTGGACa                             >  1:996687/1‑141 (MQ=255)
                                                                                                          ataataATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACaa                            >  1:452453/1‑141 (MQ=255)
                                                                                                           taataaTCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACaaa                           >  2:445328/1‑141 (MQ=255)
                                                                                                             ataatCAGAGCACGATAAATAATTACGATCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATcattttgtccaac                         <  1:233282/141‑13 (MQ=255)
                                                                                                             ataatCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATcattttgtccaac                         <  1:86595/141‑13 (MQ=255)
                                                                                                             ataatCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAat                         <  2:100809/141‑1 (MQ=255)
                                                                                                                  cAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAATGATGc                    >  1:570358/1‑141 (MQ=255)
                                                                                                                  aagagCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTATCTTCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAATGATGc                    <  1:420337/140‑1 (MQ=255)
                                                                                                                   agagCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAATGATGCt                   >  2:279457/1‑141 (MQ=255)
                                                                                                                      gCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAATGATGCTaaa                >  2:476656/1‑141 (MQ=255)
                                                                                                                        aCGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAATGATGCTAAACg              <  1:137645/141‑1 (MQ=255)
                                                                                                                         cGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTCACTCAAAAAGAAATCGTATCCAAATAAAATGAATACATCGCTGGCCACAATGCAGGTAAACGt             >  1:322252/1‑141 (MQ=255)
                                                                                                                         cGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAATGATGCTAAACGt             >  1:593019/1‑141 (MQ=255)
                                                                                                                         cGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAATGATGCTAAACGt             <  2:809283/141‑1 (MQ=255)
                                                                                                                               ataatTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAATGATGCTAAACGTAAAGTg       <  1:803742/141‑1 (MQ=255)
                                                                                                                                taatTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAATGATGCTAAACGTAAAGTgg      <  1:85637/141‑1 (MQ=255)
                                                                                                                                 aatTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAATGATGCTAAACGTAAAGTGGc     >  2:1036799/1‑141 (MQ=255)
                                                                                                                                  atTACGACCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATATATCGTTGGACAAAATGATGCTAAACGTAAAGTGGCa    <  2:487244/141‑1 (MQ=255)
                                                                                                                                   ttACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAATGATGCTAAACGTAAAGTGGCaa   >  2:925863/1‑141 (MQ=255)
                                                                                                                                    tACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAATGATGCTAAACGTAAAGTGGCAAt  >  2:975975/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                       
AACGCCATGCATCGCATTTGTCTGCTGAAGAAATGGCATATGAGAGCTTGAAAGTAGCGGCTGATATTTGTGTCTTTACCAACGATAATATTGTTGTCGAAACACTATAATAATCAGAGCACGATAAATAATTACGAGCAATTAATTTTAGTTAAAAGACGGAGGAATGAAATTAATGGATACAGCTGGAATAAGATTAACTCCAAAAGAAATCGTATCTAAATTAAATGAATACATCGTTGGACAAAATGATGCTAAACGTAAAGTGGCAAT  >  CP000730/1271995‑1272267

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