Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,509,477 C→T 10.3% pseudogene (469/471 nt) insO → KpLE2 phage‑like element; IS911B regulator fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,509,4770CT10.3% 138.6 / 4.8 66pseudogene (469/471 nt)insOKpLE2 phage‑like element; IS911B regulator fragment
Reads supporting (aligned to +/- strand):  ref base C (23/32);  new base T (3/3);  total (27/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 = 3.70e-01

CTCGCCGTTGTTCTCGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAA  >  NC_000913/4509342‑4509615
                                                                                                                                       |                                                                                                                                          
cTCGCCGTTGTTCTCGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCACGCATGAGCTTGTGATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGAt                                                                                                                                     >  1:1269278/1‑143 (MQ=255)
   gCCGTTGTTCTCGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACGCGTGGTAAGCCCGTCGGGGGGATGTTCCAAGCGAtgga                                                                                                                                  >  1:99841/1‑140 (MQ=255)
   gCCGTTGTTCTCGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCaa                                                                                                                                  <  1:708959/143‑1 (MQ=255)
    ccGTTCTTCTCGACCTGCTCGCAAGAAACCCACTGGGCTGGCCCATGTCATTCTCGCCGCACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAg                                                                                                                                 <  1:858258/143‑1 (MQ=255)
       ttgttCTCGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCa                                                                                                                              <  2:431324/143‑1 (MQ=255)
          ttCTCGACGGGGTCGCAAGAAACCCAGTGGGCGGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTc                                                                                                                           <  2:757525/143‑1 (MQ=255)
             tcGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGCAAGGGCATGGGAACCCCGTCGTAAGCCCGGCGGGGAGAGGTTCGAAGCGATCAAGGCAgccgca                                                                                                                        >  2:475323/1‑138 (MQ=255)
             tcGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGCCGGGGGGATGTTCCAAGCGCTCAAGGCAGTCAtt                                                                                                                        >  2:589995/1‑143 (MQ=255)
                     ttCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATGAGACGAgc                                                                                                                >  1:22472/1‑143 (MQ=255)
                      gcGCCAGAAAACCAGTGGGCGGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAgca                                                                                                               <  1:458884/142‑1 (MQ=255)
                       cGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGATTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCAATATACGAgcag                                                                                                              >  2:171067/1‑143 (MQ=255)
                        gCAAGAAAACCAGTGGGCAGGGCCATGGCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAGATGGCAGGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAgcagg                                                                                                             <  1:355977/143‑1 (MQ=255)
                           aGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAgcaggcag                                                                                                          >  1:908630/1‑143 (MQ=255)
                                  cAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGc                                                                                                   >  1:7519/1‑143 (MQ=255)
                                    gtaggctgggcCATGTCGTTCTCGGCGGAGCGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAg                                                                                                 <  2:1158828/140‑1 (MQ=255)
                                    gccggctgcgcCATGTCGGTCCCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAg                                                                                                 <  1:85990/140‑1 (MQ=255)
                                         ctcggcCGTGACGGTGTCGCCGGACAGCAGGCTTACCGTGAAAGCACTGGAAATGGCATGGGAAACCCGCGGTAAGCCCGGCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACt                                                                                            <  2:353956/140‑1 (MQ=255)
                                             gcCGTGGCGTTCTCGCCGGGACGCATGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTgg                                                                                        <  2:609044/143‑1 (MQ=255)
                                             gcCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGGAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTgg                                                                                        <  1:1265691/143‑1 (MQ=255)
                                                     gTTCTCGCCGGACAGCAGGCTTAACATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGACGGGGCGAAGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTCACTGTGGCGATGCCg                                                                                >  2:413244/1‑143 (MQ=255)
                                                      ttCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCgg                                                                               >  1:996881/1‑143 (MQ=255)
                                                      ttCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTAGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCgg                                                                               >  1:636820/1‑143 (MQ=255)
                                                        cgcgcgGGACGGCAGGCTTACCAGGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGGCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGAt                                                                             <  2:1189232/141‑1 (MQ=255)
                                                           ggcGGACCGCAGGCATAGCATGGAAGCGCTGGACATGGCATGGGAAACCCGTGGTGAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAg                                                                          <  2:480507/141‑1 (MQ=255)
                                                           gCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAg                                                                          >  2:1156009/1‑143 (MQ=255)
                                                            ccGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTGCTGTGGCGGCACCGGAGCAgg                                                                         >  2:359733/1‑143 (MQ=255)
                                                             cGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCGCGTCGAGGTGATGTTCCAAGCGATCAGGGCAGTCATTATACGAGCAGGCAGTGCCGGCAGTTGCTGTGGCGAGAGAGGCTCAGGc                                                                        >  2:884319/1‑143 (MQ=255)
                                                                aCAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATCTGCCAAGCTATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCaga                                                                     <  1:1212930/143‑1 (MQ=255)
                                                                  agcagGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACTGGATCAGGCAGAGt                                                                   >  1:91452/1‑143 (MQ=255)
                                                                   gcagGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCGGGCCGTTACTGTGGCGATACCGGATCAGGCAGAGTa                                                                  >  1:232887/1‑143 (MQ=255)
                                                                    cagGCTTACCATGAAAGCACTGGAAATGGCATGGGAATCCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTAt                                                                 <  2:671186/143‑1 (MQ=255)
                                                                      ggCTTACCATGAAATCACTGGAAATGGCATGGGAAACCTGTTGTAAGCCCGTCGGGGTGATTTTCTAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGa                                                               <  1:738311/143‑1 (MQ=255)
                                                                           aCCATGATAGCACTGGATATGGCATGTGAAACCCTTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCgg                                                          <  1:1012391/143‑1 (MQ=255)
                                                                            ccATGAAAGCACTGGAAATGGCATAGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGc                                                         <  1:838394/143‑1 (MQ=255)
                                                                              aTGCAAGCACTGGAAAGGGCATGGGAGACCCGTGGTAAGCCCGTCGGGGTGAGGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGt                                                       <  2:519860/143‑1 (MQ=255)
                                                                                gAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTgg                                                     >  1:380499/1‑143 (MQ=255)
                                                                                  gaGCACTGGGAATGGCATGGGAGACCCGTTGTAAGCGCGTCGGGGTGCTGCGTCACGAGATCAAGGCAGTCATCATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGaa                                                   <  2:1241704/142‑1 (MQ=255)
                                                                                   aGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGACTATGAGTCGGCGTGGaaa                                                  >  1:237362/1‑143 (MQ=255)
                                                                                   aGCACTGGAAAGGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGaaa                                                  <  2:16137/143‑1 (MQ=255)
                                                                                        tGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCGGATACAGTGCCGATACCGGATCAGGCAGAGTATGAGTCGGCGGGGCAActgct                                             >  2:1177339/1‑143 (MQ=255)
                                                                                         ggAAACAGCATGGGAAACCCGTGGTAAGCGCGTCGGGGTGAGGTTCCAAGCGATCGAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAActgctg                                            <  2:136591/143‑1 (MQ=255)
                                                                                         cgatagggCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGAGGTTCCAGGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAActgctg                                            <  2:692115/137‑1 (MQ=255)
                                                                                            aaTGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGa                                         <  2:1220908/143‑1 (MQ=255)
                                                                                             aTGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGAt                                        <  2:996881/143‑1 (MQ=255)
                                                                                                gCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGAGCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGCCGGAGTGGAAACGCCTGGGATAAc                                     >  2:1116626/1‑143 (MQ=255)
                                                                                                 cATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACa                                    <  1:287840/143‑1 (MQ=255)
                                                                                                   tGGGACACCCGTAGTAAGCCCCTCGGGGTGATGTTCCAAGCGATCAAGGCATTCATTCACCCCGCCCCTTTCTCTGCCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGc                                  <  1:998464/143‑1 (MQ=255)
                                                                                                   tGGGAAACCCGTGGTAAGCCCGACGGGGAGAAATAATAATCGGTGAATGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGCTGACAGc                                  >  1:543543/1‑143 (MQ=255)
                                                                                                         aCCCGTGGTAAGCCCGTCGGGGTGATGTCCGAAGCGATCAAGGCTGTCATTATACTGGCACGCCGTCCACCGAGTTCCTGTGGCGAGTCCGGACTAGGCACAGTATGGGCGGGCGTCGAATCTGCTGGGATAACAGCCCAATg                            <  2:312706/143‑1 (MQ=255)
                                                                                                           ccGTGGTAAGCCCGTCGGGGTGATGATCCAAACGATCAAGGCAGTAATTATACGAGCAGGCAGTTCCGGCAGTTGCTGTGGCGATACCGGATCAGGCCTAGTATGCGGCGGCGTAGACAGTGCGGGGATAACAGCCGAGTGGa                          >  2:732104/1‑143 (MQ=255)
                                                                                                            cGTGGTAAGCCCGTCGGGGTGATGTTCTAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAActgc                                              <  2:405489/122‑1 (MQ=255)
                                                                                                            cGTGGTAAGCCCGTCGGGGTGATGTTCTAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAActgc                                              >  1:405489/1‑122 (MQ=255)
                                                                                                            cGTGGTAAGCCCGTCGGGGTGATGTTCTAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAActgc                                              <  2:402464/122‑1 (MQ=255)
                                                                                                            cGTGGTAAGCCCGTCGGGGTGATGTTCTAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAActgc                                              >  1:402464/1‑122 (MQ=255)
                                                                                                             gTGGTAAGCGCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCGGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCAGCGTGGAAACTGCTGGGATAACAGCCCAATGGAgc                        <  2:628968/143‑1 (MQ=255)
                                                                                                             gTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGGGGCAACTGCTGGGATACCAGCCCAATGGAgc                        >  1:13735/1‑143 (MQ=255)
                                                                                                             gTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGGATGAGTGGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAgc                        >  2:880447/1‑143 (MQ=255)
                                                                                                                 gtaggcggTCGGGGGGACGTGCGAGGCGAGCAGGGCGGTCGTTGGACGAGCGGGGAGTTCGGGCAGTTGCTGTGGGGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGctt                    <  2:949374/136‑1 (MQ=255)
                                                                                                                     cccGTCGGGGTGATGTTCCAAGCGATGAAGGCAGTCTTTATCTGTGCAGTCGGTTGCGGCAGTTCCTGTGGCGATACCGGAGCAGGCAGAGTGTGAGTCGGCGTGGAAACTGCTGGGATAACCGCCCAATGGAGCGCttcttc                <  2:761959/143‑1 (MQ=255)
                                                                                                                       cGTCGGGGTGATGTGCGACGCGGTAGAGGCAGTCATTATACGAGCAGGCAGTGCGGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAg              <  2:437061/143‑1 (MQ=255)
                                                                                                                         gcGGGGTGATGTTCCAAGCGATCCAGGCAGTCACTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGa            <  1:403039/142‑1 (MQ=255)
                                                                                                                            gcgTGATGTCCGACGCGATCAAGGCAGTCATTATACGGGCAGGCCGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAGCTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTc         <  2:945201/141‑1 (MQ=255)
                                                                                                                              gTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTg       >  1:1013160/1‑143 (MQ=255)
                                                                                                                               tGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTTGCTTAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTga      <  1:1005449/143‑1 (MQ=255)
                                                                                                                                gATGCTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCGGCTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTgaa     <  1:1260583/143‑1 (MQ=255)
                                                                                                                                  tGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTC‑GCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTgaagaa  <  1:1156009/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                          
CTCGCCGTTGTTCTCGACCTGTTCGCAAGAAAACCAGTGGGCTGGGCCATGTCGTTCTCGCCGGACAGCAGGCTTACCATGAAAGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAA  >  NC_000913/4509342‑4509615

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