Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,387,305 G→A P124P (CCG→CCA yfbM → DUF1877 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,387,3050GA98.5% 461.3 / ‑2.6 135P124P (CCG→CCAyfbMDUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (75/58);  total (76/59)
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.18e-01

Only 100 of 137 total aligned reads displayed.

AGTTAGCCGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCGGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTATGCTGCG  >  NC_000913/2387208‑2387400
                                                                                                 |                                                                                               
aGTTAGCCGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCagg                                                                                               <  1:825420/100‑1 (MQ=255)
 gTTAGCCGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGa                                                                                              >  2:615836/1‑100 (MQ=255)
 gTTAGCCGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGa                                                                                              >  2:483095/1‑100 (MQ=255)
 gTTAGCCGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGa                                                                                              >  1:2527783/1‑100 (MQ=255)
      ccGCGACAATCGACCGACTGGAATCACTGGATCGAGATGAACGGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGAc                                                                                         <  2:876246/100‑1 (MQ=255)
      ccGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGAc                                                                                         <  1:1208338/100‑1 (MQ=255)
      ccGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGAc                                                                                         <  2:1198258/100‑1 (MQ=255)
       cgcgACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGAcc                                                                                        <  1:1511494/100‑1 (MQ=255)
        gcgACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCt                                                                                       >  2:3128043/1‑100 (MQ=255)
          gACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCttt                                                                                     <  2:1348330/100‑1 (MQ=255)
            cAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTc                                                                                   <  2:2110860/100‑1 (MQ=255)
               tCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCaga                                                                                <  1:2958151/100‑1 (MQ=255)
                 gaccgacTGGAATCACTGGATGGAAATGCATTACGTAAACAGTTTTCCATTAAGCGGCTGAATGAGATGGATTTTTATCCAGGCGTGACCTTTTCagaag                                                                              <  2:877952/100‑1 (MQ=255)
                 gaccgacTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGGAAGGGAGATTTATCCAGGAGTGACCTTTTCagaag                                                                              <  2:2239143/100‑1 (MQ=255)
                 gaccgacTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCagaag                                                                              <  2:2239365/100‑1 (MQ=255)
                   ccgacTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAagag                                                                            <  2:1700820/100‑1 (MQ=255)
                   ccgacTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAagag                                                                            <  1:2814221/100‑1 (MQ=255)
                    caacTGGAATCACTGGATCGAAATGAACTGCGTAACCAGTTTTCCATTAAGCGACTCAATGAAATGTAGATTTATCCAGGAGAGACCTTTTCAGAAGAGt                                                                           >  2:2546015/3‑100 (MQ=255)
                       cTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTag                                                                        >  1:254077/1‑100 (MQ=255)
                       cTGGAATCACTGGAGCGAAATGAACTCCGTAAACAGTTTTCCATTAACCGACTCCATGAAATGGAGATTTATCCAGGAGTGACCTTATCAGAAGTGTTat                                                                        >  2:3733660/1‑99 (MQ=255)
                        tGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTaga                                                                       >  1:1176080/1‑100 (MQ=255)
                         ggAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTagag                                                                      <  1:3364569/100‑1 (MQ=255)
                         ggAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTagag                                                                      >  1:3633131/1‑100 (MQ=255)
                         ggAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTagag                                                                      <  2:1472174/100‑1 (MQ=255)
                          taaTCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCTACTTAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAgg                                                                     <  2:1828162/99‑1 (MQ=255)
                          gAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAgg                                                                     <  2:3226697/100‑1 (MQ=255)
                           aaTCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAggg                                                                    <  2:1104484/100‑1 (MQ=255)
                           aaTCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAggg                                                                    >  2:2113956/1‑100 (MQ=255)
                            aTCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAgggg                                                                   >  2:493589/1‑100 (MQ=255)
                            aTCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAgggg                                                                   >  2:3440475/1‑100 (MQ=255)
                            aTCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAgggg                                                                   <  2:2333704/100‑1 (MQ=255)
                            aTCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAgggg                                                                   >  2:2028022/1‑100 (MQ=255)
                              cACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCa                                                                 <  1:3387498/100‑1 (MQ=255)
                                cTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAAc                                                               <  1:965683/100‑1 (MQ=255)
                                cTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAAc                                                               <  2:2845029/100‑1 (MQ=255)
                                cTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAAc                                                               <  1:806162/100‑1 (MQ=255)
                                 tGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAAct                                                              >  2:2583211/1‑100 (MQ=255)
                                  ggATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAActc                                                             >  1:4029113/1‑100 (MQ=255)
                                   gATCGAAATTAACTGCGTAAACAGTTTTGCATTAAGCGAATCAAGGAAATGGAGTTTTATCTAAGAGTTATCTTTTCAGTAGAGTTAGAGGGGTAActgt                                                            >  2:824032/1‑98 (MQ=255)
                                    aTCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCGGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCtt                                                           >  1:1003416/1‑100 (MQ=255)
                                      cGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCGGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTg                                                         <  1:2550137/100‑1 (MQ=255)
                                      cGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTg                                                         <  2:1152997/100‑1 (MQ=255)
                                      cGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTg                                                         >  1:1765759/1‑100 (MQ=255)
                                       gAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGc                                                        <  1:3040703/100‑1 (MQ=255)
                                       gAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGc                                                        >  2:2673813/1‑100 (MQ=255)
                                       gAAAGGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGc                                                        <  2:1134845/100‑1 (MQ=255)
                                        aaaTGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGCGTGACCTTTTCAGAAGAGCTAGAGGGGCAACTCTTTGcc                                                       >  2:3446392/1‑100 (MQ=255)
                                        aaaTGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGcc                                                       <  2:1150195/100‑1 (MQ=255)
                                          aTGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTc                                                     >  2:611787/1‑100 (MQ=255)
                                           tGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCg                                                    <  2:2596120/100‑1 (MQ=255)
                                             aaCTGCGTAAACAGTTTTCCATTAAGCGACTCAACGAAATGGAGATTTATCCAAGAGTGACCATTTCAGAAGAGATAGAGGGGCAACTCTTTGCCTCGAt                                                  >  1:3793467/1‑100 (MQ=255)
                                               cTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCa                                                >  2:3576505/1‑100 (MQ=255)
                                                  cGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGc                                             <  2:2513400/100‑1 (MQ=255)
                                                  cGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGc                                             >  2:2961136/1‑100 (MQ=255)
                                                   gTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCt                                            >  2:1711556/1‑100 (MQ=255)
                                                     aaaCAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTg                                          >  2:2501870/1‑100 (MQ=255)
                                                     aaaCAGTTATCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTg                                          >  2:3833066/1‑100 (MQ=255)
                                                      aaCAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGa                                         <  2:1365800/100‑1 (MQ=255)
                                                        cAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGata                                       <  2:1643297/100‑1 (MQ=255)
                                                        cAGTTTTCCAGTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGata                                       <  2:570564/100‑1 (MQ=255)
                                                         aGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGatat                                      <  1:1670704/100‑1 (MQ=255)
                                                          gTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATg                                     >  2:699124/1‑100 (MQ=255)
                                                           ttttCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATgg                                    >  1:1993786/1‑100 (MQ=255)
                                                           ttttCCATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCGTGATATgg                                    >  2:653685/1‑100 (MQ=255)
                                                               ccATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGaaaa                                >  2:532276/1‑100 (MQ=255)
                                                                cATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAg                               >  2:2877066/1‑100 (MQ=255)
                                                                cATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAg                               >  2:3082597/1‑100 (MQ=255)
                                                                cATTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAg                               <  2:75305/100‑1 (MQ=255)
                                                                 aTTAAGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGc                              >  1:3754488/1‑100 (MQ=255)
                                                                    aaGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGa                           >  1:3148391/1‑100 (MQ=255)
                                                                     aGCGACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGAt                          <  1:3633373/100‑1 (MQ=255)
                                                                        gACTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAg                       >  1:3775994/1‑100 (MQ=255)
                                                                         aCTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAgc                      >  2:1086554/1‑100 (MQ=255)
                                                                          cTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAgcg                     <  1:2007877/100‑1 (MQ=255)
                                                                          cTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAgcg                     >  1:2619552/1‑100 (MQ=255)
                                                                          cTCAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAgcg                     <  1:2954667/100‑1 (MQ=255)
                                                                            cAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGcgcg                   >  1:3833660/1‑100 (MQ=255)
                                                                            cAATGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGcgcg                   >  1:3374530/1‑100 (MQ=255)
                                                                             aaTGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGt                  >  2:518226/1‑100 (MQ=255)
                                                                             aaTGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGt                  <  1:2309810/100‑1 (MQ=255)
                                                                             aaTGAAATGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGt                  >  1:1610834/1‑100 (MQ=255)
                                                                              aTGACATGGAGTTTTATCCAGGAGTGCCCTTTTCAGCAGTGTTAGAGGGGCAACTCTTTGCCTCGACCATGCTTGCTATGGAAACGCTGATTAGCGCGTa                 >  2:2975101/1‑100 (MQ=255)
                                                                                  aaTGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGa             >  2:3193733/1‑100 (MQ=255)
                                                                                  aaTGGAGATTTAGCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGa             <  2:2313077/100‑1 (MQ=255)
                                                                                   aTGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGAc            <  1:2540896/100‑1 (MQ=255)
                                                                                   aTGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGAc            <  2:2253990/100‑1 (MQ=255)
                                                                                    tGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACg           >  1:576039/1‑100 (MQ=255)
                                                                                    tGGAGATTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACg           <  1:2629828/100‑1 (MQ=255)
                                                                                     ggAGATTTATCCAGGAGTGACCTTGTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGt          >  2:187253/1‑100 (MQ=255)
                                                                                      gagaTTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTa         <  1:3002488/100‑1 (MQ=255)
                                                                                      gagaTTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTa         <  2:1737006/100‑1 (MQ=255)
                                                                                      gagaTTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTa         >  2:1639071/1‑100 (MQ=255)
                                                                                      gagaTTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTa         <  1:292139/100‑1 (MQ=255)
                                                                                      gagaTTTATCCAAGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCACGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTa         >  2:3851073/1‑100 (MQ=255)
                                                                                        gaTTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTAtg       >  2:3130549/1‑100 (MQ=255)
                                                                                         aTTTATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTAtgc      <  1:667620/100‑1 (MQ=255)
                                                                                           ttATCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTAtgctg    >  1:2891253/1‑100 (MQ=255)
                                                                                             aTCCAGGAGTGACCTTTTCAGAAGCGTTAGAGGGGCAACTCTTTGCCTCGATCACACTCGATATGGAAAAGCTGATTAGCGCGTATCGACGTATGCTGCg  >  2:3790302/1‑100 (MQ=255)
                                                                                             aTCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTATGCTGCg  <  1:3653914/100‑1 (MQ=255)
                                                                                             aTCCAGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGATGATTAGCGCGTATCGACGTATGCTGCg  <  1:2634327/100‑1 (MQ=255)
                                                                                                 |                                                                                               
AGTTAGCCGCGACAATCGACCGACTGGAATCACTGGATCGAAATGAACTGCGTAAACAGTTTTCCATTAAGCGACTCAATGAAATGGAGATTTATCCGGGAGTGACCTTTTCAGAAGAGTTAGAGGGGCAACTCTTTGCCTCGATCATGCTTGATATGGAAAAGCTGATTAGCGCGTATCGACGTATGCTGCG  >  NC_000913/2387208‑2387400

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