Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,184,114 A→G 5.3% R281R (CGT→CGC USA300HOU_2058 ← S1 domain RNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,184,1140AG 5.3% 208.2 / 3.5 75R281R (CGT→CGCUSA300HOU_2058S1 domain RNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (38/33);  new base G (2/2);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

AACCTGTTCTAAATGCTGGATCTACGCCTAATATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGAAATCCTCTACTGATGTCGTATCGAATTCAAACTTTAC  >  CP000730/2183975‑2184254
                                                                                                                                           |                                                                                                                                            
aaCCTGTTCTAAATGCTGGCTCTACGCCTAATATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCAc                                                                                                                                             >  2:508968/1‑141 (MQ=255)
aaCCTGTTCTAAATGCTGGATCTACGCCTAATATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCAc                                                                                                                                             >  1:563127/1‑141 (MQ=255)
    tGTTCTAAATGCTGGATCTACGCCTAATATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGctc                                                                                                                                         <  1:138869/141‑1 (MQ=255)
    cgTTCTAAATGCTGGATCTACGCCTAATATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGctc                                                                                                                                         <  1:565281/140‑1 (MQ=255)
      ttCTAAATGCTGGATCTACGCCTAATATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTCAAATCAGAATGGATTTCACtctcgc                                                                                                                                       >  1:183375/1‑139 (MQ=255)
       tCTAAATGCTGGATCTACGCCTAATAGTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTACTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTAt                                                                                                                                      >  2:545318/1‑141 (MQ=255)
           aaTGCTGGATCTACGCCTAATATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGaa                                                                                                                                  >  1:917214/1‑141 (MQ=255)
             tGCTGGATCTACTCCTAGTATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAgg                                                                                                                                <  2:520106/141‑1 (MQ=255)
               cTGGATCTACGCCTAATATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACCCTCTATCGAAGGGa                                                                                                                              <  1:474873/141‑1 (MQ=255)
                  gATCTACGCCTAATATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACaa                                                                                                                           >  2:727489/1‑141 (MQ=255)
                  caacTACCCCTAATATTTGCTTACCTTTCATTGCAGGTCGCAGCAATAGATTTCTGAACTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAACCAGCATGGATTTCACGCTCTATCGAAGGGACaa                                                                                                                           <  2:240646/138‑1 (MQ=255)
                   aTCTACGCCTAATATTTGTTTACCTTTCATAGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAAt                                                                                                                          <  2:577655/141‑1 (MQ=255)
                        cGCCTAATATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCGCGCTCTATCGAAGGGACAATTAAGc                                                                                                                     >  2:740876/1‑141 (MQ=255)
                        cGCCTAATATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGACTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCGCGCTCTATCGAAGGGACAATTAAGc                                                                                                                     >  2:740860/1‑141 (MQ=255)
                          ccTAATATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGt                                                                                                                   >  2:868510/1‑141 (MQ=255)
                          ccTAATATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTAACTAAAAACATCTATTGCATGATTTTCAGCTTTTTAAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGt                                                                                                                   >  2:832691/1‑141 (MQ=255)
                              atatTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATCGATTTCACGCTCTATCGACGGTACAATTAAGCGTTTCa                                                                                                               <  1:508968/141‑1 (MQ=255)
                                 tttGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAACATGGATTTCACGCTCTATCGAAGGGACAATTAAGGGTTTCAAGc                                                                                                            >  2:166539/1‑141 (MQ=255)
                                    gTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGCTTTCAAGCTCTATCGAAGGGACACTTAAACGTTTCAAGCTAt                                                                                                         >  1:389348/1‑141 (MQ=255)
                                    gTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTAt                                                                                                         >  1:142087/1‑141 (MQ=255)
                                      ttACCTTTCATTGGCGGGTGCAACAATAGATTTCTTAAGGTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCACGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCt                                                                                                       <  2:5957/141‑1 (MQ=255)
                                       tACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCtt                                                                                                      <  2:559351/141‑1 (MQ=255)
                                       tACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCtt                                                                                                      >  2:567769/1‑141 (MQ=255)
                                       tACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCtt                                                                                                      >  2:567831/1‑141 (MQ=255)
                                       tACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCtt                                                                                                      <  2:609431/141‑1 (MQ=255)
                                       caCCTTTAAGTGGCGGGTGCAGTAATAGATTTCTTAAGTTTTCACGAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCtt                                                                                                      <  2:298529/140‑1 (MQ=255)
                                        aCCTTTCATTGGAGGTTGCAGTGATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCttt                                                                                                     <  2:92100/141‑1 (MQ=255)
                                        aCCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCttt                                                                                                     <  2:909650/141‑1 (MQ=255)
                                             tCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCg                                                                                                <  1:106053/141‑1 (MQ=255)
                                             tCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCg                                                                                                <  1:704013/141‑1 (MQ=255)
                                                 tGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCctt                                                                                            >  1:335260/1‑140 (MQ=255)
                                                    aGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTaa                                                                                         >  2:868406/1‑141 (MQ=255)
                                                        tGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGTATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAata                                                                                     <  2:146967/141‑1 (MQ=255)
                                                        tGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAGata                                                                                     >  1:611158/1‑141 (MQ=255)
                                                        ggCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGAGTTTCAGCTTTTTATGTTAAATCAGCATGGATTTCACGGTCTATCGAGGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAACata                                                                                     <  2:445175/140‑1 (MQ=255)
                                                             tAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACt                                                                                >  2:888343/1‑141 (MQ=255)
                                                                tAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCg                                                                             >  1:887420/1‑141 (MQ=255)
                                                                  gATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGAt                                                                           >  1:165778/1‑141 (MQ=255)
                                                                         tAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACatt                                                                    >  1:335438/1‑141 (MQ=255)
                                                                            gTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTGATCGGCTCTAAACTATAACTGCGATTTACattagt                                                                 >  1:42297/1‑139 (MQ=255)
                                                                            gTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATGGCCTCTAAAATATAACTGCGATTTACattgtt                                                                 >  1:824071/1‑141 (MQ=255)
                                                                              tttCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCGCGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATtattat                                                               <  1:740876/141‑1 (MQ=255)
                                                                              tttCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCGCGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATtattat                                                               <  1:740860/141‑1 (MQ=255)
                                                                               ttCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTAtg                                                              <  2:866209/141‑1 (MQ=255)
                                                                                tCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGCTTCAAGCTCTCTTTAATCGCCTCCAAAACACAACTCCGATTTACATTATTAtga                                                             >  2:144953/1‑141 (MQ=255)
                                                                                      aaaaCATCTATTGCATGATTTACAGTTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAATCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATgattgatt                                                       <  2:177583/141‑1 (MQ=255)
                                                                                       aaaCATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTGCACGCGCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCCCAAAAATATAACTGCGATTTACATTATTATGATTGAttt                                                      <  1:261022/141‑1 (MQ=255)
                                                                                          cATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCCCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCtt                                                   <  2:567622/141‑1 (MQ=255)
                                                                                            tgtATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAGGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGa                                                 <  2:197910/139‑1 (MQ=255)
                                                                                            tCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAGGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGa                                                 <  2:197744/141‑1 (MQ=255)
                                                                                                 tGCATGATTTTCAGCTCTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAAGCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGc                                            <  1:376002/141‑1 (MQ=255)
                                                                                                   cATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGGTTTATATTATTATGATTGATTTCTTGACGTGCaa                                          >  1:120199/1‑141 (MQ=255)
                                                                                                   cATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCaa                                          <  1:567769/141‑1 (MQ=255)
                                                                                                   cATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCaa                                          >  1:76685/1‑141 (MQ=255)
                                                                                                   cATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCaa                                          <  1:567831/141‑1 (MQ=255)
                                                                                                    aTGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCaag                                         <  2:790438/141‑2 (MQ=255)
                                                                                                       aTTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGaa                                      <  1:24722/141‑1 (MQ=255)
                                                                                                           tCAGCTTTTTCAGTTAAAACAGGAGGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCATTGCAATcc                                  >  1:218506/1‑141 (MQ=255)
                                                                                                            cAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGAAATCct                                 >  2:182892/1‑141 (MQ=255)
                                                                                                              gCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGAAATCctct                               >  1:855725/1‑141 (MQ=255)
                                                                                                               gcttttCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGGGATGGAATCctctc                              <  2:886143/139‑2 (MQ=255)
                                                                                                                      aGTTAAATCAGCATGGGTTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGAAATCCTCTACTGATGt                       >  1:357011/1‑141 (MQ=255)
                                                                                                                        ttAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGAAATCCTCTCCGGATGTct                     >  2:154989/1‑140 (MQ=255)
                                                                                                                        ctAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCCATGAAATCCTCTACTGTTGTct                     >  2:154998/2‑140 (MQ=255)
                                                                                                                         tAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGTGATTTACATTATTATGATTGATTTCTTGACGTGCAATGCAATCCTCTACTGATGTCGt                    <  1:545318/141‑1 (MQ=255)
                                                                                                                           aaTCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGAAATCCTCTACTGATGTCGTAt                  >  1:563264/1‑141 (MQ=255)
                                                                                                                            aTCAGCCTGGATTTCACGCGCTTTCTAAGGGACAATTAAGCGTTTCAAGCTATCTTTATTCGCCTCTAAATTATAACTGCGATTTACATTATTATGATTGATTTCTTGTCGTGCAATGAAATCCGCTACTGATGTCGTATc                 <  2:226986/141‑1 (MQ=255)
                                                                                                                             tCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGAAATCCTCTACTGCTGTCGTATCg                >  2:524161/1‑141 (MQ=255)
                                                                                                                             tCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGAAATCCTCTACTGATGTCGTATCg                >  1:819007/1‑141 (MQ=255)
                                                                                                                               aGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTCATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGAAATCCTCTACTGATGTCGTATCGat              >  1:809273/1‑140 (MQ=255)
                                                                                                                               aGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCACTGAAATCCTCTACTGATGTCGTATCGaa              >  1:502609/1‑141 (MQ=255)
                                                                                                                                 cATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGAAATCCTCTACTGATGTCGTATCGAAtt            >  2:57234/1‑141 (MQ=255)
                                                                                                                                  aTGGATTTCACGCTCTATCGAAGGGACAATTACGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTCGACGTGCAATCACATCATCTCCTGATGTCGTATCGACTTc           <  2:569745/141‑1 (MQ=255)
                                                                                                                                  aTGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGAAATCCTCTACTGATGTCGTATCGAATTc           <  1:868406/141‑1 (MQ=255)
                                                                                                                                         tCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGAAATCCTCTACTGGTGTCGTATCGAATTCAAACttt    >  1:554061/1‑141 (MQ=255)
                                                                                                                                           aCGCTCTATCGAAGGGACAATTAAGCGTATCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGAAATCCTCTACTGATGTCGTATCGAATTCAAACtttcc  >  1:576713/1‑139 (MQ=255)
                                                                                                                                           |                                                                                                                                            
AACCTGTTCTAAATGCTGGATCTACGCCTAATATTTGTTTACCTTTCATTGGAGGTTGCAGTAATAGATTTCTTAAGTTTTCACTAAAAACATCTATTGCATGATTTTCAGCTTTTTCAGTTAAATCAGCATGGATTTCACGCTCTATCGAAGGGACAATTAAGCGTTTCAAGCTATCTTTAATCGCCTCTAAAATATAACTGCGATTTACATTATTATGATTGATTTCTTGACGTGCAATGAAATCCTCTACTGATGTCGTATCGAATTCAAACTTTAC  >  CP000730/2183975‑2184254

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: