Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,403,483 A→G A137A (GCA→GCG accB → acetyl‑CoA carboxylase BCCP subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,403,4830AG99.5% 585.4 / ‑6.5 189A137A (GCA→GCGaccBacetyl‑CoA carboxylase BCCP subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (106/82);  total (106/83)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 220 total aligned reads displayed.

GAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTA  >  NC_000913/3403393‑3403572
                                                                                          |                                                                                         
  aaGTCAACGTGGGCGATACCCTGTGCATCGTTGAATCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTc                                                                                       <  1:4366365/93‑1 (MQ=255)
  aaGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTc                                                                                       >  1:3978561/1‑93 (MQ=255)
  aaGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTc                                                                                       <  2:3685734/93‑1 (MQ=255)
  aaGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTc                                                                                       <  1:778621/93‑1 (MQ=255)
   aGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCt                                                                                      <  1:2908352/93‑1 (MQ=255)
    gTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCCTGAAAGCGATTCTg                                                                                     >  1:204430/1‑93 (MQ=255)
    atCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTg                                                                                     <  1:2939840/92‑1 (MQ=255)
     tCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTgg                                                                                    <  1:2235447/93‑1 (MQ=255)
     tCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTgg                                                                                    <  1:4187811/93‑1 (MQ=255)
     tCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTgg                                                                                    >  2:2935783/1‑93 (MQ=255)
      cAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGt                                                                                   <  2:2463397/93‑1 (MQ=255)
       aaCGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTc                                                                                  >  1:4224538/1‑93 (MQ=255)
       aaCGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGCGAAAGCGATTCTGGTc                                                                                  >  2:836701/1‑93 (MQ=255)
          gTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGaa                                                                               <  2:2034140/93‑1 (MQ=255)
             ggCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGt                                                                            <  2:1585494/93‑1 (MQ=255)
              gCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTg                                                                           <  1:61425/93‑1 (MQ=255)
                gATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGTCAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGa                                                                         >  2:4274566/1‑93 (MQ=255)
                gATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGa                                                                         <  2:861726/93‑1 (MQ=255)
                 aTACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGAc                                                                        >  1:3653438/1‑93 (MQ=255)
                  tACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACa                                                                       >  2:3969504/1‑93 (MQ=255)
                   aCCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACaa                                                                      >  1:4217897/1‑93 (MQ=255)
                   aCCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACaa                                                                      >  2:464337/1‑93 (MQ=255)
                     ccTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAAcc                                                                    <  2:2684073/93‑1 (MQ=255)
                     ccTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAAcc                                                                    >  1:3486659/1‑93 (MQ=255)
                      cTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCg                                                                   <  2:4218091/93‑1 (MQ=255)
                      cTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCg                                                                   >  1:3223428/1‑93 (MQ=255)
                       tgtgCATCGTTGAAGCCATGAAAATGATGAATCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGCTTCTGGTCGAAAGTGGACAACCgg                                                                  >  2:1872929/1‑93 (MQ=255)
                       tgtgCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCgg                                                                  >  1:4549867/1‑93 (MQ=255)
                         tgCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTa                                                                <  2:3620548/93‑1 (MQ=255)
                         tgCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTa                                                                <  2:2424248/93‑1 (MQ=255)
                         tgCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTa                                                                <  2:3540714/93‑1 (MQ=255)
                         tgCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTa                                                                >  2:4423414/1‑93 (MQ=255)
                          gCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAg                                                               >  2:1771703/1‑93 (MQ=255)
                          gCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGGGAAATCGCTTCTGGTCGAAAGTGGCAAACCggaca                                                               >  1:2643622/1‑90 (MQ=255)
                            aTCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGaa                                                             >  1:1545266/1‑93 (MQ=255)
                             tCGTTGAAGCGATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAAt                                                            >  2:1521002/1‑93 (MQ=255)
                             tCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAAt                                                            >  2:674281/1‑93 (MQ=255)
                             tCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAAt                                                            >  2:4102770/1‑93 (MQ=255)
                             tCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAAt                                                            >  1:1910459/1‑93 (MQ=255)
                               gTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAAttt                                                          >  2:172218/1‑93 (MQ=255)
                               gTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAAttt                                                          >  1:1224996/1‑93 (MQ=255)
                               gTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAAttt                                                          >  1:47412/1‑93 (MQ=255)
                                ttGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTg                                                         >  1:2132070/1‑93 (MQ=255)
                                 tGAAGCCATGAAAATGATGAACAAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGa                                                        <  1:4486825/93‑1 (MQ=255)
                                   aaGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACg                                                      >  1:2528955/1‑93 (MQ=255)
                                     gCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAg                                                    >  2:1877952/1‑93 (MQ=255)
                                          gAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGcg                                               >  2:3817776/1‑92 (MQ=255)
                                           aaaaTGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTg                                              >  2:181552/1‑93 (MQ=255)
                                            aaaTGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTgg                                             >  2:4233909/1‑93 (MQ=255)
                                             aaTGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGgt                                            <  1:1339213/93‑1 (MQ=255)
                                             aaTGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGgt                                            >  1:233332/1‑93 (MQ=255)
                                              atgatgAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGgtc                                           >  2:3522116/1‑93 (MQ=255)
                                               tgatgaACCAGATCGAAGCGGACAAATCCGGTCCCGGGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGgtcg                                          <  2:2972964/93‑1 (MQ=255)
                                               agatgaACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGgtcg                                          <  1:1658521/92‑1 (MQ=255)
                                                gatgaACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGgtcgt                                         >  2:3466553/1‑93 (MQ=255)
                                                gatgaACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGgtcgt                                         >  2:2886304/1‑93 (MQ=255)
                                                gatgaACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGgtcgt                                         >  2:3557021/1‑93 (MQ=255)
                                                  tgaACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGgtcgtct                                       >  1:270167/1‑92 (MQ=255)
                                                  tgaACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCa                                       <  1:1990616/93‑1 (MQ=255)
                                                  tgaACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCa                                       <  2:3650183/93‑1 (MQ=255)
                                                     aCCAGATCGAAGCGGACAAACCGGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCg                                    <  1:2487159/93‑1 (MQ=255)
                                                     aCCACATGGAAGCGGACAAATCCGGTCCCGTGAAAGCGATTCTGGTCGAAAGTGGACACCCGGTAGAATTTGACGAGCCGCTGGTCGTCATCg                                    <  2:3402797/93‑1 (MQ=255)
                                                        aGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGt                                 <  2:4529865/93‑1 (MQ=255)
                                                        aGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGt                                 >  1:3141675/1‑93 (MQ=255)
                                                        aGATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGt                                 >  1:3431372/1‑93 (MQ=255)
                                                        aGATCGAAGCGCAAAAATCCGGTACCGGGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGt                                 <  2:3093492/93‑1 (MQ=255)
                                                         gATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTa                                >  1:3062843/1‑93 (MQ=255)
                                                         gATCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTa                                >  1:71959/1‑93 (MQ=255)
                                                          aTCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTaa                               <  2:956708/93‑1 (MQ=255)
                                                          aTCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTaa                               >  2:1984729/1‑93 (MQ=255)
                                                           tCGAAGCGGACAACTCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAAc                              <  2:370927/93‑1 (MQ=255)
                                                           tCGAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAAc                              <  1:708296/93‑1 (MQ=255)
                                                           acggaGCGGACAAACCCGGTACCGTAAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAAc                              <  1:966218/89‑1 (MQ=255)
                                                             gagacGGACAAATAGCGTACTGTCAAAGATATTCTGGTGAAAAGGGAACACCCGGTAGAATTTGTTGAGCCTCTGGTCTTCATCGAGTAACGa                            <  2:937554/89‑1 (MQ=255)
                                                             gAAGCGGACAAGAGCGATACCGTGAAGGCGATTCTGGGCGAAAGTGGAAACCCGATAGAATTTGGCGAGCCGCTGGTCGTCATCGAGTAACGa                            <  1:831636/93‑1 (MQ=255)
                                                             gAAGCGGACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGa                            >  2:3751672/1‑93 (MQ=255)
                                                             gAAGCGCACAAGTCCGGAACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGa                            <  2:1995247/93‑1 (MQ=255)
                                                                gtggACAAATCGGGTACCGTGAAAGCGATTCGGGTCGAAAGTGGACAACCGGTAAAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGc                         <  1:1048158/91‑1 (MQ=255)
                                                                  ggACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGa                       >  2:3223623/1‑93 (MQ=255)
                                                                   gACAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGac                      >  1:1441755/1‑92 (MQ=255)
                                                                     cAAATCGGGTACCGTGAAAGCGATTCTGGCCGAAGGTCGACAACGGGAAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACa                    <  1:2731883/93‑1 (MQ=255)
                                                                     cAAATCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACa                    <  1:49830/93‑1 (MQ=255)
                                                                      aaaTCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACAt                   <  2:1325478/93‑1 (MQ=255)
                                                                        aTCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGc                 >  2:584720/1‑93 (MQ=255)
                                                                        aTCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGc                 >  1:1085729/1‑93 (MQ=255)
                                                                        aTCCGGTACCGTGAAAGCGATTCTGGTCGAAAGTGAACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGc                 <  1:4279723/93‑1 (MQ=255)
                                                                          ccGGTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTg               <  1:690713/93‑1 (MQ=255)
                                                                             gTACCGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGGCGTCATCGAGTAACGAGGCGAACATGCTgggt            >  2:3283280/1‑91 (MQ=255)
                                                                               tccGGGAAAGGGATTTTGGTGAAAAGTGGGAAACCGGTAAAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATaa          <  2:4482995/92‑1 (MQ=255)
                                                                                ccGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATaaa         >  1:1617973/1‑93 (MQ=255)
                                                                                 cGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATaaaa        <  2:2259979/93‑1 (MQ=255)
                                                                                 cGTGAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATaaaa        <  2:2184205/93‑1 (MQ=255)
                                                                                  gaaaaaGCGTTTTTGGTGAAAAGTGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAAtt      <  2:3079809/90‑1 (MQ=255)
                                                                                   tGAAAGCGATTCTGGTCGAAGGTGGACACCCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAAtt      <  2:1294399/93‑1 (MQ=255)
                                                                                    gAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTg     >  1:374930/1‑93 (MQ=255)
                                                                                    gAAAGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTg     >  2:1667504/1‑93 (MQ=255)
                                                                                      aaGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGtt   >  1:974721/1‑93 (MQ=255)
                                                                                      aaGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGtt   >  1:4098271/1‑93 (MQ=255)
                                                                                       aGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTa  <  1:185958/93‑1 (MQ=255)
                                                                                       aGCGATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTa  <  2:428444/93‑1 (MQ=255)
                                                                                          |                                                                                         
GAAAGTCAACGTGGGCGATACCCTGTGCATCGTTGAAGCCATGAAAATGATGAACCAGATCGAAGCGGACAAATCCGGTACCGTGAAAGCAATTCTGGTCGAAAGTGGACAACCGGTAGAATTTGACGAGCCGCTGGTCGTCATCGAGTAACGAGGCGAACATGCTGGATAAAATTGTTA  >  NC_000913/3403393‑3403572

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