Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,988,989 T→G A26A (GCA→GCC hemX ← putative uroporphyrinogen III methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,988,9890TG88.7% 457.5 / 28.1 170A26A (GCA→GCChemXputative uroporphyrinogen III methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (11/8);  new base G (71/78);  total (84/86)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 172 total aligned reads displayed.

TTACCCCAGCCATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTTGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCATTATGGCTTCCTGTTATGAG  >  NC_000913/3988891‑3989085
                                                                                                  |                                                                                                
   ccccAGCCATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGa                                                                                              >  2:2571028/1‑100 (MQ=255)
   ccccAGCCATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGa                                                                                              >  2:1698720/1‑100 (MQ=255)
   ccccAGCCATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGa                                                                                              >  1:1740716/1‑100 (MQ=255)
    cccAGCCATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATGGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGAc                                                                                             >  1:900167/1‑100 (MQ=255)
    cccAGCCATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGAc                                                                                             >  1:653539/1‑100 (MQ=255)
       aGCCATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAgg                                                                                          >  1:876610/1‑100 (MQ=255)
       aGCCATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAgg                                                                                          <  1:2255253/100‑1 (MQ=255)
        gCCATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGGATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGt                                                                                         >  1:3108402/1‑100 (MQ=255)
        gCCATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGt                                                                                         <  1:1155861/100‑1 (MQ=255)
        gCCATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGt                                                                                         >  1:497064/1‑100 (MQ=255)
        gCCATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGt                                                                                         <  1:3317521/100‑1 (MQ=255)
         ccATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGtt                                                                                        >  2:3875760/1‑100 (MQ=255)
         ccATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGtg                                                                                        >  2:2688436/1‑99 (MQ=255)
          cATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTtg                                                                                       <  2:1195184/100‑1 (MQ=255)
          cATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTtg                                                                                       >  1:3649538/1‑100 (MQ=255)
          cATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTtg                                                                                       <  2:3605312/100‑1 (MQ=255)
           aTACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTtgt                                                                                      <  1:951512/100‑1 (MQ=255)
           aTACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTtgt                                                                                      >  1:1620300/1‑100 (MQ=255)
           aTACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCACGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTtgt                                                                                      >  1:3954974/1‑100 (MQ=255)
             aCAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGa                                                                                    <  2:3359879/100‑1 (MQ=255)
              cAACCCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGTGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGAc                                                                                   <  2:2405891/100‑1 (MQ=255)
              cAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTTGCGACAGGTTGTGAc                                                                                   >  1:2145306/1‑100 (MQ=255)
              cAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGAc                                                                                   <  1:2074131/100‑1 (MQ=255)
              cAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGAc                                                                                   >  1:1186946/1‑100 (MQ=255)
                aaCCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACgt                                                                                 >  2:3833593/1‑100 (MQ=255)
                aaCCGATGCCCGACGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAACACAACCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACGGGTTGTGACgt                                                                                 >  2:385111/1‑100 (MQ=255)
                 aCCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACgtg                                                                                >  1:3476566/1‑100 (MQ=255)
                  ccGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTTGCGACAGGTTGTGACgtgg                                                                               <  2:1457264/100‑1 (MQ=255)
                  ccGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACgtgg                                                                               <  2:851248/100‑1 (MQ=255)
                  ccGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACgtgg                                                                               >  2:3523706/1‑100 (MQ=255)
                  ccGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACgtgg                                                                               >  1:3392691/1‑100 (MQ=255)
                     aTGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGgtgt                                                                            <  2:3879962/100‑1 (MQ=255)
                      tGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTc                                                                           >  2:2153192/1‑100 (MQ=255)
                      tGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTc                                                                           <  1:3437432/100‑1 (MQ=255)
                      tGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTc                                                                           <  2:3653169/100‑1 (MQ=255)
                        cccGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCa                                                                         >  1:2645696/1‑100 (MQ=255)
                         ccgccgCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTTTTTCTACAGGTTGTGACGTGGTGTCCAc                                                                        >  2:1456385/1‑100 (MQ=255)
                          cgccgcCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACg                                                                       <  1:738024/100‑1 (MQ=255)
                          cgccgcCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGGTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACg                                                                       >  2:2824529/1‑100 (MQ=255)
                             cgccAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGcc                                                                    <  2:2512457/100‑1 (MQ=255)
                               ccAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTc                                                                  <  1:1929511/100‑1 (MQ=255)
                                cAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTTTTTTTTTTCAGGCTGTGACGTGGTTTCCACGGCCTcc                                                                 >  2:3337776/1‑100 (MQ=255)
                                  gagCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTTTTTTTTTTCCGGATGTGACGTGGTGTACACGGCCTCCCt                                                               >  2:1497999/1‑100 (MQ=255)
                                    gCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTgg                                                             <  1:1872504/100‑1 (MQ=255)
                                    gCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTgg                                                             >  2:1643939/1‑100 (MQ=255)
                                    gCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTgg                                                             <  1:2318689/100‑1 (MQ=255)
                                    gCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTgg                                                             >  2:1145642/1‑100 (MQ=255)
                                    gCAATAGCGATAGCCACCGCGCCGAGAATCAATGCGGTATTGTCCTTACTCTTTTTTTCTGTCGCGACTGGTTGTAACGAGGACTCCACGGCCTCCATgg                                                             >  1:1257076/1‑100 (MQ=255)
                                      aaTAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGtc                                                           <  2:674234/100‑1 (MQ=255)
                                      aaTAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGtc                                                           <  2:3401077/100‑1 (MQ=255)
                                       aTAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGtct                                                          >  2:1976980/1‑100 (MQ=255)
                                         aGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTctct                                                        >  1:450965/1‑100 (MQ=255)
                                           cGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTTGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTc                                                      <  1:3481510/100‑1 (MQ=255)
                                            gttAGCCACCGCGGTGAGAACCAAAGCGGCATAGTTCTTACTCTTTTTTTCTGTAGCGATAGGTTGTGACGCGGTGTCCACGGCCTCCCTGGTCTCTTcg                                                     >  2:3037885/3‑99 (MQ=255)
                                            gATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCa                                                     <  1:1600792/100‑1 (MQ=255)
                                             aTAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGATGTGGTGTCCACGGCCTCCCTGGTCTCTTCaa                                                    <  1:1691549/100‑1 (MQ=255)
                                              tAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAAc                                                   <  1:2778716/100‑1 (MQ=255)
                                              tAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAAc                                                   <  1:3490035/100‑1 (MQ=255)
                                                gCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCa                                                 >  1:3133863/1‑100 (MQ=255)
                                                  cACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGACGGCAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCAcg                                               >  2:1266375/1‑100 (MQ=255)
                                                    ccGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCGCGGCCTCCCTGGTCTCTTCAACCAcggc                                             >  1:3380373/1‑100 (MQ=255)
                                                     cgcgCTGAGAATCAATGCGGTATTGTTCTTACGCTCTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCAcggcg                                            <  1:1247249/100‑1 (MQ=255)
                                                      gcgcTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCAcggcgg                                           >  1:2093402/1‑100 (MQ=255)
                                                         cTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAgg                                        >  1:1497735/1‑100 (MQ=255)
                                                           gagaATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGtt                                      >  1:2494105/1‑100 (MQ=255)
                                                             gaATCAATGCGGTATTGTTCTTACTCTTTTTTTTTTTTTTTTCAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGtttt                                    >  2:3739889/1‑100 (MQ=255)
                                                             gaATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGtttt                                    >  1:1582229/1‑100 (MQ=255)
                                                               aTCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGtttttt                                  <  1:3022697/100‑1 (MQ=255)
                                                                tCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTc                                 <  2:517808/100‑1 (MQ=255)
                                                                 cAATGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCCCGGCCTCCCTGGTCTCTTCAACCACGGCGGTGGTTTTTTCt                                >  2:2538699/1‑100 (MQ=255)
                                                                 cAAGGCGGTATTGTTCTTACTCTTTTTTGCTGGGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCt                                <  2:1168884/100‑1 (MQ=255)
                                                                  aaTGCGGTATTGTTCTTACTCTTTTTTTTTTTTTTTTCAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCtt                               >  2:3305528/1‑100 (MQ=255)
                                                                  aaTGCGGTATTGTTCTTACTCTTTTTTTCTGTTGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCtt                               <  1:2702500/100‑1 (MQ=255)
                                                                  aaTGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCtt                               >  1:3833752/1‑100 (MQ=255)
                                                                   aTGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTg                              <  2:2146165/100‑1 (MQ=255)
                                                                   aTGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTg                              <  1:2545562/100‑1 (MQ=255)
                                                                   aTGCGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTg                              <  2:1012299/100‑1 (MQ=255)
                                                                     gCGGTATTGTTCTTACTCTTTTTTTCTTTTGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGCTTTTTCTTGtt                            >  2:1913466/1‑100 (MQ=255)
                                                                      cGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGtta                           >  2:894509/1‑99 (MQ=255)
                                                                      cGGTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTc                           <  1:636986/100‑1 (MQ=255)
                                                                        gTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCg                         <  1:3585075/100‑1 (MQ=255)
                                                                        gTATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCg                         <  1:4029067/100‑1 (MQ=255)
                                                                         tATTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGt                        <  1:1673957/100‑1 (MQ=255)
                                                                          aTTGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTc                       <  2:2497256/100‑1 (MQ=255)
                                                                           ttGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCa                      <  1:1843043/100‑1 (MQ=255)
                                                                           ttGTTCTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCa                      <  1:3202150/100‑1 (MQ=255)
                                                                           ttGTTCTGACTCTTTTTTTCTGTTGCGACAGGTTGTGACGTGGTGTCCACGGCATCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCa                      <  2:3719593/100‑1 (MQ=255)
                                                                             gTTCTTACTCTTTTTTTTTTTTTTTTTAGGTTGTGACGTGGTGTCCACGGCCTCTCTGGTCTCTTAAACAAGGGCGGGGGTTTTTTCTTGTTCTGTCGtt                    >  2:524319/1‑100 (MQ=255)
                                                                                cTTACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCATTATg                 >  1:2909160/1‑100 (MQ=255)
                                                                                 ttACTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCATTATgg                <  1:1118411/100‑1 (MQ=255)
                                                                                   aCTCTTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCATTATGGCt              <  2:324286/100‑1 (MQ=255)
                                                                                      ctTTTTTTCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCATTATGGCTTcc           >  2:3405414/1‑100 (MQ=255)
                                                                                        ttttttCTGTGGCGACCGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCATTATGGCTTCCTg         <  1:3358927/100‑1 (MQ=255)
                                                                                        tcttgtCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCGGGTCTCTTCAACCACGGCGGAGGTTTTTTCTCGTTCCGTCATTATGGCTTCCTg         <  2:582039/95‑1 (MQ=255)
                                                                                           tttCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCATTATGGCTTCCTGTTa      <  2:407801/100‑1 (MQ=255)
                                                                                           tttCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCATTATGGCTTCCTGTTa      >  2:2495741/1‑100 (MQ=255)
                                                                                            ttCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCATTATGGCTTCCTGTTAt     <  1:1443306/100‑1 (MQ=255)
                                                                                             tCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCATTATGGCTTCCTGTTATg    <  2:543913/100‑1 (MQ=255)
                                                                                             tCTGTGGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCATTATGGCTTCCTGTTATg    <  2:639949/100‑1 (MQ=255)
                                                                                               tgtgGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTTGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCATTATGGCTTCCTGTTATgag  <  2:2656280/100‑1 (MQ=255)
                                                                                                  |                                                                                                
TTACCCCAGCCATACAAACCGATGCCCGCCGCCAGAGCAATAGCGATAGCCACCGCGCTGAGAATCAATGCGGTATTGTTCTTACTCTTTTTTTCTGTTGCGACAGGTTGTGACGTGGTGTCCACGGCCTCCCTGGTCTCTTCAACCACGGCGGAGGTTTTTTCTTGTTCCGTCATTATGGCTTCCTGTTATGAG  >  NC_000913/3988891‑3989085

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: