Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,582,538 A→G G336G (GGT→GGC yffG ← putative oxidoreductase Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,582,5380AG98.4% 380.7 / ‑3.3 123G336G (GGT→GGCyffGputative oxidoreductase Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (52/69);  total (53/70)
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.04e-01

Only 100 of 129 total aligned reads displayed.

GAGCAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGACCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTGAAA  >  NC_000913/2582448‑2582626
                                                                                           |                                                                                        
gAGCAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGcc                                                                                        <  1:1487301/93‑1 (MQ=255)
gAGCAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGcc                                                                                        >  1:1583572/1‑93 (MQ=255)
gAGCAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACCCCATTGCGGGCCAGAACATCCGCACAGGCCAGCCCGGCCGGGcc                                                                                        >  1:1531583/1‑93 (MQ=255)
 aGCAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCt                                                                                       <  1:2845851/93‑1 (MQ=255)
  gCAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTg                                                                                      >  1:1089080/1‑93 (MQ=255)
  gCAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTg                                                                                      >  1:3907632/1‑93 (MQ=255)
     aGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCAc                                                                                   >  2:15294/1‑93 (MQ=255)
     aGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCAc                                                                                   <  2:59577/93‑1 (MQ=255)
         aCCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGAt                                                                               >  1:1577208/1‑93 (MQ=255)
           cGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATaa                                                                             <  1:4273802/93‑1 (MQ=255)
            gATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAAt                                                                            >  1:734739/1‑93 (MQ=255)
            gATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAAt                                                                            >  2:2063861/1‑93 (MQ=255)
            gATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAAt                                                                            <  1:3155904/93‑1 (MQ=255)
            gATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAAt                                                                            <  2:3659665/93‑1 (MQ=255)
            gATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAAt                                                                            <  2:4247452/93‑1 (MQ=255)
            caTTTCTGGATGACGCTCTTACACCGTCCCCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAAt                                                                            <  2:4122065/92‑1 (MQ=25)
                  tGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCAc                                                                      >  2:163432/1‑93 (MQ=255)
                  tGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCAc                                                                      <  2:1285911/93‑1 (MQ=255)
                  tGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCAc                                                                      <  1:166805/93‑1 (MQ=255)
                   ggAGGGCGATCGTACACCGTCCCCCCCACGCCATTGCGGGTCAGACCACCCGCACAGGCAAGCCCTGCCGGACCTGCCCCGATAATCGCCAcc                                                                     <  2:3929696/93‑1 (MQ=255)
                    gATGGCGATCGTCCACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCAccc                                                                    <  1:29105/93‑1 (MQ=255)
                    gATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCAccc                                                                    <  1:1814324/93‑1 (MQ=255)
                     aTGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCg                                                                   >  1:2416095/1‑93 (MQ=255)
                      tGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGc                                                                  <  1:2426830/93‑1 (MQ=255)
                       ggCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCt                                                                 <  2:2105869/93‑1 (MQ=255)
                       ggCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCt                                                                 >  2:2688/1‑93 (MQ=255)
                           aTCGTACACCGTCCCCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTc                                                             <  1:2610256/93‑1 (MQ=255)
                           aTCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTc                                                             >  1:3592792/1‑93 (MQ=255)
                           aTCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTc                                                             >  2:4042925/1‑93 (MQ=255)
                             cGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCa                                                           >  2:3672107/1‑93 (MQ=255)
                             cGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCa                                                           <  2:1839465/93‑1 (MQ=255)
                              gTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCAc                                                          <  2:3047804/93‑1 (MQ=255)
                                acacCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACtt                                                        >  2:566034/1‑93 (MQ=255)
                                 cacCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACttt                                                       <  2:3763121/93‑1 (MQ=255)
                                  acCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTg                                                      >  2:343722/1‑93 (MQ=255)
                                   ccGTCGCCCCCGCGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACGTTgg                                                     <  2:3029562/93‑1 (MQ=255)
                                   ccGTCACCCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTg                                                      <  1:1200422/93‑1 (MQ=255)
                                   ccGTCACCCCCTCGCCAGTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCTTGCACCGATAATTGCCACCCGCTTGTCCACTTTgg                                                     <  1:836475/93‑1 (MQ=255)
                                   ccGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTgg                                                     >  2:876740/1‑93 (MQ=255)
                                   ccGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTgg                                                     >  2:2165492/1‑93 (MQ=255)
                                   ccGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTgg                                                     <  1:1618662/93‑1 (MQ=255)
                                    cGTCCCCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGt                                                    <  2:1039141/93‑1 (MQ=255)
                                      tCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGa                                                  >  2:2760466/1‑93 (MQ=255)
                                       cACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGAc                                                 <  1:1315375/93‑1 (MQ=255)
                                        aCCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGACCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACa                                                >  1:4029463/1‑93 (MQ=255)
                                          cccccACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACAt                                               <  1:3723639/93‑1 (MQ=255)
                                           ccccACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGACCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACGTg                                              >  1:4339294/1‑93 (MQ=255)
                                            cccACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATgg                                             >  2:107402/1‑93 (MQ=255)
                                             ccACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGc                                            <  1:2117233/93‑1 (MQ=255)
                                              cACGCCATTGCGGGTCAGAACATCCGCACAGGACAGCCCTGCCGGGCCTGCACCGATAATCGCAACCCGCTTGTCCACTTTGGATACATGGCt                                           >  2:1654915/1‑93 (MQ=255)
                                               aCGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCtt                                          >  2:1898061/1‑93 (MQ=255)
                                               aCGCCATGGCGGGTCAGAACAGCCGCACATGACAGCCATGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCtt                                          <  1:730760/93‑1 (MQ=255)
                                                cGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTa                                         <  1:719916/93‑1 (MQ=255)
                                                  ccATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAg                                       <  1:1855921/93‑1 (MQ=255)
                                                  ccATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAg                                       <  2:1575215/93‑1 (MQ=255)
                                                  ccATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAg                                       <  1:2670626/93‑1 (MQ=255)
                                                     ttGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCa                                    >  2:634709/1‑93 (MQ=255)
                                                     ttGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCa                                    >  1:2303631/1‑93 (MQ=255)
                                                      tGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAg                                   <  2:2795779/93‑1 (MQ=255)
                                                        cGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCaagc                                 >  1:2823291/1‑90 (MQ=255)
                                                        cGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGa                                 <  1:300728/93‑1 (MQ=255)
                                                        cGGGCCAGAACATCCCCCCAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGa                                 <  2:4083226/93‑1 (MQ=255)
                                                         gggTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGAc                                <  2:1315715/93‑1 (MQ=255)
                                                          ggTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACg                               >  2:1049187/1‑93 (MQ=255)
                                                          ggTCAGAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACg                               >  1:2590070/1‑93 (MQ=255)
                                                               gAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAAc                          >  1:899047/1‑93 (MQ=255)
                                                               gAACATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAAc                          >  2:2405038/1‑93 (MQ=255)
                                                                aaCATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAAcc                         >  1:3731590/1‑93 (MQ=255)
                                                                aaCATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAAcc                         >  1:1782894/1‑93 (MQ=255)
                                                                aaCATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAAcc                         <  2:3113055/93‑1 (MQ=255)
                                                                aaCATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAAcc                         >  2:3104467/1‑93 (MQ=255)
                                                                 aCATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCt                        >  2:2476135/1‑93 (MQ=255)
                                                                  cATCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCtt                       <  1:419282/93‑1 (MQ=255)
                                                                   acccGCACAGGCCGGCCCTGCCGGGCGTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCttt                      <  2:2686723/91‑1 (MQ=38)
                                                                   aTCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCttt                      <  1:118041/93‑1 (MQ=255)
                                                                   aTCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCttt                      <  1:1019917/93‑1 (MQ=255)
                                                                   aTCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCttt                      <  1:2651300/93‑1 (MQ=255)
                                                                   aTCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCttt                      <  1:2649350/93‑1 (MQ=255)
                                                                    tCCGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTc                     <  2:4308327/93‑1 (MQ=255)
                                                                      cGCACAGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGc                   >  1:3226821/1‑93 (MQ=255)
                                                                         ccaGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCaa                <  1:1740049/92‑1 (MQ=255)
                                                                          caGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAAc               >  1:2254546/1‑93 (MQ=255)
                                                                          caGGCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACGTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAAc               <  2:408499/93‑1 (MQ=255)
                                                                           aGGCGAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACg              <  2:815477/93‑1 (MQ=255)
                                                                            ggCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGc             >  2:960314/1‑93 (MQ=255)
                                                                            ggCCAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGc             >  2:3772555/1‑93 (MQ=255)
                                                                              ccAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCt           <  1:3303590/93‑1 (MQ=255)
                                                                               cAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTg          <  1:1975559/93‑1 (MQ=255)
                                                                               cAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTg          <  2:2183396/93‑1 (MQ=255)
                                                                               cAGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTg          >  1:2043379/1‑93 (MQ=255)
                                                                                aGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGa         <  1:36463/93‑1 (MQ=255)
                                                                                aGCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGa         >  1:357416/1‑93 (MQ=255)
                                                                                 gCCCTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGAt        >  2:2695285/1‑93 (MQ=255)
                                                                                  ccacgCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCGCTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATc       <  2:345014/89‑1 (MQ=255)
                                                                                   ccTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCt      >  2:4389869/1‑93 (MQ=255)
                                                                                   ccTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCt      <  1:3730869/93‑1 (MQ=255)
                                                                                   ccTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCt      >  1:3606575/1‑93 (MQ=255)
                                                                                    cTGCCGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTg     >  1:4270989/1‑93 (MQ=255)
                                                                                       ccGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTGaaa  <  1:2197097/93‑1 (MQ=255)
                                                                                       ccGGGCCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTGaaa  <  1:4059534/93‑1 (MQ=255)
                                                                                           |                                                                                        
GAGCAAGCCACCGATTTCTGGATGGCGATCGTACACCGTCACCCCCACGCCATTGCGGGTCAGAACATCCGCACAGGCCAGCCCTGCCGGACCTGCACCGATAATCGCCACCCGCTTGTCCACTTTGGTGACATGGCTTAAGTCAGGACGCCAACCTTTCGCCAACGCCTGATCTGAAA  >  NC_000913/2582448‑2582626

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