Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,341,955 G→A T15T (ACG→ACA melA → alpha‑galactosidase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,341,9550GA99.1% 373.0 / ‑4.7 110T15T (ACG→ACAmelAalpha‑galactosidase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (50/58);  total (52/58)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 110 total aligned reads displayed.

CTGCAGATTCGCCTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACGATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGG  >  NC_000913/4341859‑4342049
                                                                                                |                                                                                              
cTGCAGATTCGCCTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAAtt                                                                                             >  2:215267/1‑100 (MQ=255)
cTGCAGATTCGCCTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAAtt                                                                                             >  2:216124/1‑100 (MQ=255)
  gCAGATTCGCCTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAAtttt                                                                                           >  2:1276581/1‑100 (MQ=255)
  gCACATTCGCCTGCCATGATGAAGTTATTCAAACAAGCCAGGAGATCTGCATGATGTCTGCACCGAAAATTACCTTTATCGGCGCTTGTTCGACAAttt                                                                                            >  2:1128881/1‑99 (MQ=255)
    aGATTCGCCTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCg                                                                                         >  1:1534914/1‑100 (MQ=255)
    aGATTCGCCTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCg                                                                                         <  2:2911297/100‑1 (MQ=255)
     gATTCGCCTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGt                                                                                        <  1:143090/100‑1 (MQ=255)
      aTTCGCCTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGtt                                                                                       <  1:1696526/100‑1 (MQ=255)
       ttCGCCTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTa                                                                                      <  2:664513/100‑1 (MQ=255)
        tCGCCTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTaa                                                                                     <  2:2173685/100‑1 (MQ=255)
        tCGCCTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATGTTCGTTaa                                                                                     >  2:327473/1‑100 (MQ=255)
         cGCCTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTaaa                                                                                    <  2:2984846/100‑1 (MQ=255)
          gCCTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTaaaa                                                                                   >  2:1946168/1‑100 (MQ=255)
           ccTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTaaaaa                                                                                  <  1:1654193/100‑1 (MQ=255)
               ccATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATAtt                                                                              >  2:2235753/1‑100 (MQ=255)
                cATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTc                                                                             <  2:1547635/100‑1 (MQ=255)
                 atgatgAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCt                                                                            <  1:1097872/100‑1 (MQ=255)
                 atgatgAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCt                                                                            <  1:2403298/100‑1 (MQ=255)
                  tgatgaAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCtt                                                                           <  1:247867/100‑1 (MQ=255)
                  tgatgaAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCtt                                                                           <  1:1946186/100‑1 (MQ=255)
                   gatgaAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTg                                                                          <  1:1676640/100‑1 (MQ=255)
                    atgaAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTgg                                                                         <  1:185668/100‑1 (MQ=255)
                     tgaAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGt                                                                        <  2:388667/100‑1 (MQ=255)
                     tgaAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGt                                                                        <  2:1213299/100‑1 (MQ=255)
                     tgaAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGt                                                                        >  1:282229/1‑100 (MQ=255)
                       gaGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACGATTTTCGTTAAAAATATTCTTGGTGa                                                                      >  2:1667687/2‑100 (MQ=255)
                       aaGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGa                                                                      <  2:2423340/100‑1 (MQ=255)
                        aGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGAt                                                                     >  1:118086/1‑100 (MQ=255)
                        aGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGAt                                                                     >  2:1075876/1‑100 (MQ=255)
                        aGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGAt                                                                     >  1:1249476/1‑100 (MQ=255)
                        aGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGAt                                                                     >  1:2591195/1‑100 (MQ=255)
                        aGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGAt                                                                     >  1:1520320/1‑100 (MQ=255)
                           tATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGAtgtg                                                                  <  2:237696/100‑1 (MQ=255)
                           tATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGAtgtg                                                                  >  2:1580405/1‑100 (MQ=255)
                             ttCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGtt                                                                >  2:2164851/1‑100 (MQ=255)
                             ttCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGtt                                                                <  2:2558128/100‑1 (MQ=255)
                              tCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTc                                                               >  2:487746/1‑100 (MQ=255)
                               caagcaagCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTcc                                                              >  2:2213866/1‑100 (MQ=255)
                               caagcaagCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTcc                                                              <  1:1620175/100‑1 (MQ=255)
                               caagcaagCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTcc                                                              >  2:105541/1‑100 (MQ=255)
                               caagcaagCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTcc                                                              >  1:469887/1‑100 (MQ=255)
                                aagcaagcCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCa                                                             >  1:1136158/1‑100 (MQ=255)
                                 agcaagcCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCAt                                                            <  2:1617048/100‑1 (MQ=255)
                                  gcaagcCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATc                                                           <  1:1187350/100‑1 (MQ=255)
                                   caagcCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATcg                                                          >  1:1883307/1‑100 (MQ=255)
                                    aagcCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATcgc                                                         <  2:2838391/100‑1 (MQ=255)
                                     agcCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATcgcg                                                        >  2:572589/1‑100 (MQ=255)
                                     agcCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATcgcg                                                        >  2:527662/1‑100 (MQ=255)
                                       ccAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAg                                                      <  1:181856/100‑1 (MQ=255)
                                         aGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGgc                                                    >  2:1685328/1‑100 (MQ=255)
                                           gagaTCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGgcgc                                                  <  1:3004624/100‑1 (MQ=255)
                                            ngaTCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACCATTTTCGTTGAAAATATTCTTGGTAATGTGTTCCATCGCGAGGCGCt                                                 >  2:2047165/2‑100 (MQ=255)
                                             gaTCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTg                                                >  2:1376309/1‑100 (MQ=255)
                                              aTCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTTGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCt                                                 <  2:2523312/98‑1 (MQ=255)
                                              aTCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGGGTTCCATCGCGAGGCGCTGa                                               <  2:1851252/100‑1 (MQ=255)
                                               tCTGCATGATGTCTGCATCCAAAATTACATTTATCGGCGCTGGTTCGACAAGTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGaa                                              <  2:2744309/100‑1 (MQ=255)
                                               tCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATagcgcgaagctgc                                               >  1:2833426/1‑90 (MQ=255)
                                               tCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGaa                                              <  1:2561948/100‑1 (MQ=255)
                                                cTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTCTCGTTACAAATATTCTTGGTGATGTGTTCCATAGCGAGCCGCTGaaa                                             >  2:622828/1‑100 (MQ=255)
                                                 tGCATTATGTCGGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGgcg                                                   <  2:2591172/93‑1 (MQ=255)
                                                  gCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAAc                                           >  2:1221961/1‑100 (MQ=255)
                                                  gCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAAc                                           >  1:799522/1‑100 (MQ=255)
                                                   cATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGTGCTGAAAACg                                          >  2:1998840/1‑100 (MQ=255)
                                                   cATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACg                                          <  2:1427354/100‑1 (MQ=255)
                                                    atgatgTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACgg                                         <  2:282226/100‑1 (MQ=255)
                                                    atgatgTCTGCACCCAAAATTACATCTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACgg                                         <  2:1553506/100‑1 (MQ=255)
                                                     tgatgTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGgc                                        <  2:1803385/100‑1 (MQ=255)
                                                      gatgTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGgcg                                       >  1:1121927/1‑100 (MQ=255)
                                                      gatgTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGgcg                                       <  2:1897741/100‑1 (MQ=255)
                                                      gatgTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGgcg                                       >  1:436945/1‑100 (MQ=255)
                                                        tgTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCa                                     <  1:1961385/100‑1 (MQ=255)
                                                        tgTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCa                                     >  2:2817500/1‑100 (MQ=255)
                                                          tCTGCACCCAAACTTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCata                                   >  1:1987036/1‑100 (MQ=255)
                                                          tCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCata                                   <  1:1276593/100‑1 (MQ=255)
                                                          tCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCata                                   >  2:1580388/1‑100 (MQ=255)
                                                           cTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCatat                                  <  1:2819551/100‑1 (MQ=255)
                                                             gCACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTg                                >  2:498298/1‑100 (MQ=255)
                                                              cACCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGc                               <  2:2164648/100‑1 (MQ=255)
                                                               aCCGAAAATTACATTTATCGGCGCTGGTGCGACAATTTTCGTTAACAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGcc                              <  1:1402184/100‑1 (MQ=255)
                                                               aCCCAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGcc                              >  2:1643766/1‑100 (MQ=255)
                                                                cccAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGccc                             <  1:918916/100‑1 (MQ=255)
                                                                cccAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGccc                             >  2:133920/1‑100 (MQ=255)
                                                                 ccAAAATTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCt                            >  1:1696420/1‑100 (MQ=255)
                                                                     aaTTACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATg                        <  2:1534901/100‑1 (MQ=255)
                                                                       ttACATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGa                      <  1:1870516/100‑1 (MQ=255)
                                                                         aCATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACa                    <  2:1114656/100‑1 (MQ=255)
                                                                          cATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACAt                   >  1:1341177/1‑100 (MQ=255)
                                                                          cATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACAt                   <  1:946874/100‑1 (MQ=255)
                                                                          cATTTATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACAt                   >  2:2766220/1‑100 (MQ=255)
                                                                            tttATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTg                 <  1:2240565/100‑1 (MQ=255)
                                                                              tATCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTGAc               >  2:183230/1‑100 (MQ=255)
                                                                               aTCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTGAcc              <  1:2164870/100‑1 (MQ=255)
                                                                               aTCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCACTGAAAACGGCGCATATTGCCCTGATGGACATTGAcc              <  1:527668/100‑1 (MQ=255)
                                                                               aTCGGCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTAGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTGAcc              <  1:2258543/100‑1 (MQ=255)
                                                                                  ggCGCTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTGACCCCa           >  1:402518/1‑100 (MQ=255)
                                                                                   gcgcTGGTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTGACCCCAc          <  1:1108036/100‑1 (MQ=255)
                                                                                        ggTTCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTGACCCCACCCGcc     >  1:1603263/1‑100 (MQ=255)
                                                                                          ttCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTg   <  1:327476/100‑1 (MQ=255)
                                                                                          ttCGACAATTTTCGTTAAAAATATTCTTGGCGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTg   <  2:1341649/100‑1 (MQ=255)
                                                                                           tCGACAATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTgg  <  1:2619884/100‑1 (MQ=255)
                                                                                                |                                                                                              
CTGCAGATTCGCCTGCCATGATGAAGTTATTCAAGCAAGCCAGGAGATCTGCATGATGTCTGCACCCAAAATTACATTTATCGGCGCTGGTTCGACGATTTTCGTTAAAAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAACGGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGG  >  NC_000913/4341859‑4342049

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