Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 821,301 T→C 5.5% P346P (CCA→CCG)  pepT ← M20 family peptidase T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730821,3010TC 5.5% 205.4 / 2.9 73P346P (CCA→CCG) pepTM20 family peptidase T
Reads supporting (aligned to +/- strand):  ref base T (31/37);  new base C (2/2);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

TTTCCATAACATCAATGGATGCATATTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATCAACTTTAACTGGATAATTTTCAAAATGGGCATTGATATCGTC  >  CP000730/821165‑821441
                                                                                                                                        |                                                                                                                                            
tttCCATAACATCAATGGATGCATATTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGaa                                                                                                                                          <  2:934677/141‑1 (MQ=255)
   ccATAACATCAATGGATGCATATTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAGTGACACTTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAatat                                                                                                                                       <  2:811649/141‑1 (MQ=255)
    cATAACATCAATGGATGCATATTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGCTGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATGGGTTCAGGGTTTGCTGGAATATc                                                                                                                                      >  1:801417/1‑141 (MQ=255)
    cATAACATCAATGGATGCATATTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATc                                                                                                                                      <  1:626068/141‑1 (MQ=255)
      tAACATCAATGGATGCATATTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTa                                                                                                                                    >  1:106031/1‑141 (MQ=255)
         cATCAATGGATGCATATTCATATGGGCCGTGGATGTTACTGCAACCTGTAAAGTGGTGTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAAtt                                                                                                                                 >  1:931885/1‑141 (MQ=255)
                ggATGCATTGTCATATGGGCCGCGGCAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCaaa                                                                                                                          <  1:170026/141‑1 (MQ=255)
                 gATGCATATTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAg                                                                                                                         <  2:672256/141‑1 (MQ=255)
                 gATGCATATTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAg                                                                                                                         <  1:411058/141‑1 (MQ=255)
                  aTGCATATTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGa                                                                                                                        <  1:673135/141‑1 (MQ=255)
                   tGCATATTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGac                                                                                                                       >  1:945966/1‑141 (MQ=255)
                   tGCATATTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGac                                                                                                                       >  1:753287/1‑141 (MQ=255)
                      attttCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACg                                                                                                                    <  2:51836/138‑1 (MQ=255)
                      atatTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACg                                                                                                                    <  2:470240/141‑1 (MQ=255)
                      atatTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACg                                                                                                                    <  2:383143/141‑1 (MQ=255)
                      atatTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAACTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCCAAgccccc                                                                                                                    >  2:470841/1‑136 (MQ=255)
                         ttCCTATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACATGttt                                                                                                                 <  1:898828/141‑1 (MQ=255)
                         ttCATATGGGCCGTGGAAGTTACCGCTACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCGAATTTGGCAAAGACACGttt                                                                                                                 <  1:570888/141‑1 (MQ=255)
                         ttCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACATGttt                                                                                                                 <  1:898838/141‑1 (MQ=255)
                         ttCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGttt                                                                                                                 <  1:199717/141‑1 (MQ=255)
                         ttCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGttt                                                                                                                 <  2:608955/141‑1 (MQ=255)
                               tGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACTCGTTTAGGtat                                                                                                           >  1:94573/1‑141 (MQ=255)
                                  gCCGTCGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATc                                                                                                        <  2:508507/141‑1 (MQ=255)
                                    cGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTATGCTGGAATATCTAATTTGGCAACGACACGTTTAGGTATATCaa                                                                                                      <  1:574655/141‑1 (MQ=255)
                                     tttgAAGTTACTGAAACCTGTAAAGAGGTTTGGAGTTGGTAACCCTATAAATGACAATTGTGAGCCAACTGTACCACCGTGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCaat                                                                                                     <  2:840537/138‑1 (MQ=255)
                                       ggAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCCGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCaataa                                                                                                   >  1:475588/1‑141 (MQ=255)
                                       ggAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCCGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCaataa                                                                                                   >  1:719378/1‑141 (MQ=255)
                                        gAAGTTAGCGCAACCTGTAAAGATGTTTGGAGATGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTGGCTGGAATATCTAATTTGGCAAAGACACGTTGAGGTATTTCaataat                                                                                                  >  1:590262/1‑141 (MQ=255)
                                         aaGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACCATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACAAGATTAGGAATATCatacgta                                                                                                 >  2:830046/1‑135 (MQ=255)
                                           gTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAGTAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCACTAATATg                                                                                               >  1:802936/1‑141 (MQ=255)
                                           gTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATg                                                                                               >  1:88270/1‑141 (MQ=255)
                                            ttACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGa                                                                                              <  2:39994/141‑1 (MQ=255)
                                            ttACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGa                                                                                              >  2:482880/1‑141 (MQ=255)
                                             tACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAg                                                                                             <  2:942869/141‑1 (MQ=255)
                                              aCCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGATGGAATATCTAATTTGGCAAAGACGCGTTTAGTTATATGAATAATATGAgg                                                                                            <  2:789023/141‑1 (MQ=255)
                                               ccGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCTCTAATATGATGc                                                                                           >  1:97484/1‑141 (MQ=255)
                                                    aCCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATgg                                                                                      <  1:105017/141‑1 (MQ=255)
                                                        gTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATTTGAGGCAATGGTAat                                                                                  >  1:805961/1‑141 (MQ=255)
                                                        gTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAat                                                                                  <  1:43998/141‑1 (MQ=255)
                                                        gTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATAGCAATAATATGAGGCACTGGGAat                                                                                  >  1:719330/1‑141 (MQ=255)
                                                           aaGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATAtt                                                                               >  2:753421/1‑141 (MQ=255)
                                                                gTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATTTGAGGCAATGGTAAttttttttt                                                                          >  1:915815/1‑140 (MQ=255)
                                                                gTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATAttttttt                                                                          >  1:108286/1‑140 (MQ=255)
                                                                    ggAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTACAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAAAGGTAAAATGTTTTCTGcc                                                                      >  2:393651/1‑141 (MQ=255)
                                                                     gAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCa                                                                     >  1:161322/1‑141 (MQ=255)
                                                                      aGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCat                                                                    >  1:622022/1‑141 (MQ=255)
                                                                        ttgtgaagcccATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCCGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCatat                                                                  <  2:719378/133‑1 (MQ=255)
                                                                        ttGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCCGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCatat                                                                  <  2:475588/141‑1 (MQ=255)
                                                                          ggTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTACCATAtta                                                                >  2:78970/1‑140 (MQ=255)
                                                                             aaCCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATAtttgca                                                             >  2:508235/1‑137 (MQ=255)
                                                                             aaCCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGaaa                                                             >  1:769484/1‑141 (MQ=255)
                                                                                 ccATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTTGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATAtt                                                         <  2:931885/141‑1 (MQ=255)
                                                                                   aTAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGa                                                       <  1:753421/141‑1 (MQ=255)
                                                                                       aTGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTTCCATCTTGAAATATTGATCCg                                                   <  1:482880/141‑1 (MQ=255)
                                                                                          aCAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGata                                                <  1:663570/141‑1 (MQ=255)
                                                                                          aCAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGata                                                <  1:850097/141‑1 (MQ=255)
                                                                                             aTTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGatatat                                             <  2:92125/141‑1 (MQ=255)
                                                                                              ttGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATc                                            <  1:407322/141‑1 (MQ=255)
                                                                                              ttGTGAACCATCTGTACCACCGCCAGTAGGTGCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATc                                            <  2:703990/141‑1 (MQ=255)
                                                                                                       ttCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATCAACTTTAAc                                   <  2:536741/140‑1 (MQ=255)
                                                                                                                   gcgAATTGGTTCAGGTTTTGCTGGCATATCTAATTTGGCAACGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATCAACTTTAACTGGATAATTTTc                       <  2:414998/141‑1 (MQ=255)
                                                                                                                   gcgAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATCAACTTTAACTGGATAATTTTc                       <  1:886535/141‑1 (MQ=255)
                                                                                                                     gAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTACCATATTGAAATATTGATCCGATATATCAACTTTAACTGGATAATTTTCaa                     <  2:106031/141‑1 (MQ=255)
                                                                                                                        tAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATAACAACTTTAACTGGATAATTTTCAAAAt                  >  2:382157/1‑141 (MQ=255)
                                                                                                                         aGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAACCATTTATGGGGTATATTGAAATATTGATCCGATATATCAACTTTCACTGGATAATTTTCCAAATg                 >  1:448102/1‑141 (MQ=255)
                                                                                                                            ttCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATCAACTTTCACTGGATAATTTTCAAAATGGGc              >  2:815648/1‑141 (MQ=255)
                                                                                                                                ttgtTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTCTTCTGCCATATTGAAATATTGATCCGATATATCAACTTTAACTGGATAATTTTCAAAATGGGCATTg          <  2:721633/140‑1 (MQ=255)
                                                                                                                                 tgtTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATCAACTTTAACTGGATAATTTTCAAAATGGTCAttta         >  1:777168/1‑139 (MQ=255)
                                                                                                                                 tgtTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATCAACTTTAACTGGATAATTTTCAAAATGGGCATTGa         >  2:170548/1‑141 (MQ=255)
                                                                                                                                  gtTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATCAACTTTAACTGGATAATTTTCAAAATGGGCATTGat        <  1:205344/141‑1 (MQ=255)
                                                                                                                                   tttGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATCAACTTTAACTGGATAATTTTCAAAATGGGCATTGata       >  2:500449/1‑141 (MQ=255)
                                                                                                                                   tttGCGTGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATCAACTTTAACTGGATAATTTTCAAAATGGGCATTGata       <  2:926384/141‑1 (MQ=255)
                                                                                                                                     tGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATCAACTTTAACTGGATAATTTTCAAAATGGGCATTGatatg     >  2:785105/1‑140 (MQ=255)
                                                                                                                                        tGGAATATCTAATTTGGCAAAGACACGTTTAGGTATAGCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATCAACTTTAACTGGATAATTTTCAAAATGGGCATTGATATCGTc  <  2:945966/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                            
TTTCCATAACATCAATGGATGCATATTCATATGGGCCGTGGAAGTTACCGCAACCTGTAAAGATGTTTGGAGTTGGTAACCCCATAAATGACAATTGTGAACCATCTGTACCACCGCGAATAGGTTCAGTGTTTGCTGGAATATCTAATTTGGCAAAGACACGTTTAGGTATATCAATAATATGAGGCAATGGTAATATTTTTTCTGCCATATTGAAATATTGATCCGATATATCAACTTTAACTGGATAATTTTCAAAATGGGCATTGATATCGTC  >  CP000730/821165‑821441

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: ‑