Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,783,940 T→C 6.6% M381V (ATG→GTG)  thrS ← threonine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,783,9400TC 6.6% 223.7 / 4.7 76M381V (ATG→GTG) thrSthreonine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (35/36);  new base C (2/3);  total (37/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 = 8.31e-01

AGGGTCTCTATAACTTAATCTAAAGCTATAATCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTGGTTTATTCGCATAAATCATCATATGATGTGGACAGTTC  >  CP000730/1783800‑1784079
                                                                                                                                             |                                                                                                                                           
aGGGTCTCTATAACTTAATCTAAAGCTATAATCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCAt                                                                                                                                             <  2:324711/141‑1 (MQ=255)
     ctctATAACTTAATCTACAGCTATAATCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCAc                                                                                                                                        <  1:239915/141‑1 (MQ=255)
      tctATAACTTAATCTAAAGCTATAATCCTCGAAACCAAAGTCTTTATACACATCAGTAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCACACCACg                                                                                                                                       >  2:232038/1‑141 (MQ=255)
      tctATAACTTAATCTAAAGCTATAATCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTAGAACAAAGATATGTGAATCATTTAAAGTCAaacaaca                                                                                                                                       >  2:264404/1‑140 (MQ=255)
      tctATAACTTAATCTAAAGCTATAATCATAGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCATTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCACACCACg                                                                                                                                       >  2:232639/1‑141 (MQ=255)
        tataACTTAATCTAAAGCTATAATCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGaa                                                                                                                                     <  1:306148/141‑1 (MQ=255)
        tacaaCTTAATCTAAAGCTATAATCCTCGAAACCAAAGTCTTTATACACATAAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATCCCACGaa                                                                                                                                     <  1:507654/138‑1 (MQ=255)
            aCTTAATCTAAAGCTATAATCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGCGCTTTGGAAAGTCATACCACGAACACg                                                                                                                                 <  1:800543/141‑1 (MQ=255)
            aCTTAATCTAAAGCTATAATCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACg                                                                                                                                 <  1:132998/141‑1 (MQ=255)
            aCTTAATCTAAAGCTATAATCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACg                                                                                                                                 <  2:801728/141‑1 (MQ=255)
            aCTTAATCTAAAGCTACAATCCTCTAAACCAAAGTCTTTATACACATCAATACTCATGTTTACAACGCGTTTGAATTCTTCCTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACg                                                                                                                                 <  1:620308/141‑1 (MQ=255)
                aaTTTAAAGCTATAATCATCGAAACCAATGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGt                                                                                                                             <  2:609915/141‑1 (MQ=255)
                aaTCTAAAGCTATAATGCTCGAAACAAAAGTCTTTATACACATCAATATTCATGTTTACGATGCGATTGAATTCTTCCTTAATTTGCAAAAGTAGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTga                                                                                                                             >  2:397732/1‑140 (MQ=255)
                aaTCTAAAGCTATAATCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGTACAAAGATATTTGACTCTTTTAAATTCATACCTCGAACACGTTGt                                                                                                                             >  1:74227/1‑141 (MQ=255)
                aaTCTAAAGCTATAATCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGt                                                                                                                             >  2:715159/1‑141 (MQ=255)
                aaTCTAAAGCTATAATCCTCGAAACCAAAGTCTTTATACACATCAATAAACATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGt                                                                                                                             <  1:233440/141‑1 (MQ=255)
                aaTCTAAAGCTATAAACCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGt                                                                                                                             >  2:727347/1‑141 (MQ=255)
                 aTCTAAAGCTATAATCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTa                                                                                                                            >  1:765638/1‑141 (MQ=255)
                 aTCTAAAGCTATAATCCTCGAAACCAAAGTCTCTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGACTAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTa                                                                                                                            <  1:540853/141‑1 (MQ=255)
                           ataATCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATat                                                                                                                  >  1:354016/1‑140 (MQ=255)
                           ataATCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATAc                                                                                                                  >  2:460194/1‑141 (MQ=255)
                           ataATCCTCGAAACCAAAGTATTTATACACATCAATAATCATGTCTACAACGCGTTTGCATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTGAAGTCATGCCACGAACACATTGTAATCCTGATAc                                                                                                                  >  1:491399/1‑141 (MQ=255)
                              aTCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAGCGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGc                                                                                                               >  1:402501/1‑141 (MQ=255)
                              aTCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGc                                                                                                               >  1:56547/1‑141 (MQ=255)
                                ccccGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGcac                                                                                                             <  2:20583/138‑1 (MQ=255)
                                ccTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGcac                                                                                                             <  2:724552/141‑1 (MQ=255)
                                ccTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTAGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGcac                                                                                                             <  1:351633/141‑1 (MQ=255)
                                  ccGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGcacca                                                                                                           <  1:346217/140‑1 (MQ=255)
                                   cgataCCAAAGTCTTTACACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAATGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGcaccac                                                                                                          <  2:542518/137‑1 (MQ=255)
                                    gAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCCCt                                                                                                         >  2:444516/1‑141 (MQ=255)
                                      aaCCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTg                                                                                                       >  2:20066/1‑141 (MQ=255)
                                       aCCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGc                                                                                                      >  1:779101/1‑141 (MQ=255)
                                       aCCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGc                                                                                                      >  1:351244/1‑141 (MQ=255)
                                        ccAAAGTCTTTATACACATCCCTAACATCAGAGACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCt                                                                                                     >  1:799857/1‑141 (MQ=255)
                                          aaaGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTc                                                                                                   <  1:435542/141‑1 (MQ=255)
                                           aaGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCa                                                                                                  <  2:765638/141‑1 (MQ=255)
                                                  tataCACATCCATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTAt                                                                                           <  2:438286/141‑1 (MQ=255)
                                                  tataCACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTAt                                                                                           >  1:806556/1‑141 (MQ=255)
                                                  tataCACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTAt                                                                                           <  2:416321/141‑1 (MQ=255)
                                                  tagagccctccataatCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAGTCCTGATACAGCACCACTTGCTTCATATCTAt                                                                                           <  2:280868/130‑1 (MQ=255)
                                                   ataCACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATg                                                                                          <  2:783599/141‑1 (MQ=255)
                                                   ataCACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATg                                                                                          <  1:203978/141‑1 (MQ=255)
                                                        caTCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCg                                                                                     >  1:42706/1‑141 (MQ=255)
                                                                atCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGc                                                                             >  1:724910/1‑141 (MQ=255)
                                                                   aTGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATTTATGCATCGTTCCTAGCTCa                                                                          <  1:493178/141‑1 (MQ=255)
                                                                   aTGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCa                                                                          <  2:201109/141‑1 (MQ=255)
                                                                           cAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCAGATCTATCCATCGTTCCTAGCTAAGcactac                                                                   >  2:504363/1‑136 (MQ=255)
                                                                            aaCGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACg                                                                  >  2:484301/1‑141 (MQ=255)
                                                                                   ttGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGt                                                           <  1:757369/141‑1 (MQ=255)
                                                                                      aaTTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAAt                                                        <  1:577165/141‑1 (MQ=255)
                                                                                       aTTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGAAACGGTTAGGTaact                                                       >  2:153421/1‑139 (MQ=255)
                                                                                              tttAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGGCATCCCACGAACACGTTGTAAGCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGata                                                <  2:18727/141‑1 (MQ=255)
                                                                                                 aaTTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATTCGGATAGGTAATTCACGCTATGa                                             >  1:57596/1‑141 (MQ=255)
                                                                                                 aaTTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGa                                             >  1:441598/1‑141 (MQ=255)
                                                                                                  aTTTGACCAGGTCGAACAAAGATATGTGAATCAATTAAAGTCATACCACGAACACGTTGTATTCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTATCTCAGCGATACGGATAGGTAATTCACGATATGaa                                            <  2:551301/141‑1 (MQ=255)
                                                                                                   tttGAACAGGTCGAACAAAGATATGTGAATCATTTAAAGTCACACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAAt                                           <  1:232639/141‑1 (MQ=255)
                                                                                                   tctGATCATGTCGAACAAAGATATGTGAATCATTTAAAGTCACACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAAt                                           <  1:264404/139‑1 (MQ=255)
                                                                                                   tctGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCACACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAAt                                           <  1:232038/139‑1 (MQ=255)
                                                                                                    ttGATCAGGTCGAGCAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAAtg                                          <  1:8178/141‑1 (MQ=255)
                                                                                                    ttGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACCGCACCACTTGCTTCATATATATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAAtg                                          <  1:155155/141‑1 (MQ=255)
                                                                                                    ttGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAAtg                                          <  1:155201/141‑1 (MQ=255)
                                                                                                    ttGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAAtg                                          <  1:137159/141‑1 (MQ=255)
                                                                                                       aTCAGGTTGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTgg                                       >  2:681216/1‑141 (MQ=255)
                                                                                                          aGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTGGgtt                                    >  1:289300/1‑141 (MQ=255)
                                                                                                          aGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATAGGAATGTGGttt                                    >  1:324371/1‑141 (MQ=255)
                                                                                                           ggTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTGGTTTa                                   >  1:251298/1‑141 (MQ=255)
                                                                                                             tCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTTCGACACGATATGAATGTGGTTTAtt                                 <  2:779101/141‑1 (MQ=255)
                                                                                                             tCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTCTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTGGTTTAtt                                 <  2:351244/141‑1 (MQ=255)
                                                                                                                        tatGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCTTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTGGTTTATTCGCATAAatca                      >  1:546722/1‑141 (MQ=255)
                                                                                                                              aaTCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTGGTTTATTCGCATAAATCATCatat                <  1:123968/141‑1 (MQ=255)
                                                                                                                                  aTTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTGGTTTATTCGCATAAATCATCATatgatg            >  2:225593/1‑141 (MQ=255)
                                                                                                                                  aTTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTGGTTTATTCGCATAAATCATCATacgatg            >  2:783336/1‑141 (MQ=255)
                                                                                                                                  aTTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTGGTTTATTCGCATAAATCATCAGctgatt            >  2:34682/1‑140 (MQ=255)
                                                                                                                                     tAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTGGTTTATTCGCATAAATCATCATATGATGTgg         >  2:248036/1‑141 (MQ=255)
                                                                                                                                     cagattCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTGGTTTATTCGCATAAATCATCATATGATGTgg         <  2:806556/136‑1 (MQ=255)
                                                                                                                                       aaGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTGGTTTATTCGCATAAATCATCATATGATGAGGAc       >  2:640918/1‑141 (MQ=255)
                                                                                                                                            aTACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATCTGAATGTGGTTTATTCGCATAAATCATCATAAGATTTGGACAGATc  >  2:583517/1‑141 (MQ=255)
                                                                                                                                            aTACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTGGTTTATTAGCATAAAACATCATCTGATGTGGACAGGTc  >  2:602116/1‑141 (MQ=255)
                                                                                                                                             |                                                                                                                                           
AGGGTCTCTATAACTTAATCTAAAGCTATAATCCTCGAAACCAAAGTCTTTATACACATCAATAATCATGTTTACAACGCGTTTGAATTCTTCTTTAATTTGATCAGGTCGAACAAAGATATGTGAATCATTTAAAGTCATACCACGAACACGTTGTAATCCTGATACAGCACCACTTGCTTCATATCTATGCATCGTTCCTAGCTCAGCGATACGGATAGGTAATTCACGATATGAATGTGGTTTATTCGCATAAATCATCATATGATGTGGACAGTTC  >  CP000730/1783800‑1784079

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