Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 356,492 A→G 77.7% T108A (ACG→GCG)  codA → cytosine/isoguanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913356,4920AG77.7% 334.4 / 79.8 164T108A (ACG→GCG) codAcytosine/isoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (20/17);  new base G (50/77);  total (70/94)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

GGAATCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGT  >  NC_000913/356397‑356589
                                                                                               |                                                                                                 
ggAATCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCt                                                                                               >  1:3192465/1‑100 (MQ=255)
ggAATCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCt                                                                                               >  2:3649837/1‑100 (MQ=255)
ggAATCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCt                                                                                               >  2:1711974/1‑100 (MQ=255)
 gAATCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTg                                                                                              <  1:2619089/100‑1 (MQ=255)
 gAATCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTg                                                                                              <  1:92541/100‑1 (MQ=255)
 gAATCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTg                                                                                              <  2:1106146/100‑1 (MQ=255)
  aaTCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGa                                                                                             >  1:3487319/1‑100 (MQ=255)
    tCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGaaa                                                                                           >  2:3170975/1‑100 (MQ=255)
       gtccggcACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATgg                                                                                        <  1:3538208‑M2/93‑1 (MQ=255)
       gtccggcACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATgg                                                                                        <  1:3907690‑M2/93‑1 (MQ=255)
        tccggcACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGc                                                                                       >  1:3825769‑M2/7‑100 (MQ=255)
         ccggcACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGc                                                                                       <  2:1304007‑M2/94‑1 (MQ=255)
           ggcACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGa                                                                                    >  1:1044782‑M2/4‑100 (MQ=255)
           ggcACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGAGCATGGCAAGCGCTGAAATGGCAGa                                                                                    >  2:3165021‑M2/4‑100 (MQ=255)
              aCGCTGTTTGAAGGCAGTGAACGCTGGGCAGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCGGAAATGGCAGAttg                                                                                 >  2:2038827‑M2/1‑97 (MQ=255)
               cGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGAttgc                                                                                <  1:81225‑M2/100‑5 (MQ=255)
                gCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGAttgcc                                                                               <  1:3356048‑M2/100‑6 (MQ=255)
                    tttGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACg                                                                           >  2:881517/1‑100 (MQ=255)
                     ttGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACgg                                                                          <  2:3565950/100‑1 (MQ=255)
                       gAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCa                                                                        <  1:1408567/100‑1 (MQ=255)
                       gAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCa                                                                        >  1:3004722/1‑100 (MQ=255)
                       gAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCa                                                                        <  1:2262265/100‑1 (MQ=255)
                        aaGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCAt                                                                       <  1:2154320/100‑1 (MQ=255)
                        aaGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCAt                                                                       >  1:1901729/1‑100 (MQ=255)
                          ggCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTc                                                                     <  1:3316870/100‑1 (MQ=255)
                          ggCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTc                                                                     <  1:1059084/100‑1 (MQ=255)
                            cATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGCCGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAg                                                                   >  2:2766365/1‑100 (MQ=255)
                            cATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAg                                                                   <  1:42455/100‑1 (MQ=255)
                             aTTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGc                                                                  >  2:2345321/1‑100 (MQ=255)
                              ttGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCa                                                                 <  2:63927/100‑1 (MQ=255)
                                gAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCAtg                                                               >  2:3077889/1‑100 (MQ=255)
                                  aCGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCAtgtg                                                             <  2:1794368/100‑1 (MQ=255)
                                  aCGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCAtgtg                                                             <  1:497831/100‑1 (MQ=255)
                                   cGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGc                                                            >  1:1396476/1‑100 (MQ=255)
                                    gCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCg                                                           <  1:2024769/100‑1 (MQ=255)
                                    gCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCg                                                           <  1:589530/100‑1 (MQ=255)
                                    gCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCg                                                           <  2:3824485/100‑1 (MQ=255)
                                    gCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCg                                                           <  2:3079000/100‑1 (MQ=255)
                                     cTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGt                                                          <  1:2577153/100‑1 (MQ=255)
                                      tGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTa                                                         <  2:2643770/100‑1 (MQ=255)
                                      tGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGATATGGCAGATTGCCAACGGCATTCAGCATGTGCg                                                           <  2:1242232/98‑1 (MQ=255)
                                        ggCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTAc                                                        <  1:1130708/99‑1 (MQ=255)
                                         gCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTAcc                                                       <  2:454833/99‑1 (MQ=255)
                                          cgatctgcAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCa                                                     <  1:3772192/94‑1 (MQ=255)
                                          ccGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCa                                                     <  1:3955498/100‑1 (MQ=255)
                                          ccGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCa                                                     >  1:3685925/1‑100 (MQ=255)
                                              gcgcAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTc                                                 <  1:3603481/100‑1 (MQ=255)
                                              gcgcAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTc                                                 <  1:850735/100‑1 (MQ=255)
                                                gcAACGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGCGCGTACCCCTGTCg                                                >  2:589384/1‑99 (MQ=255)
                                                gcAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCaagatcg                                               >  2:3772164/1‑94 (MQ=255)
                                                  aaaGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATg                                             <  1:3565967/100‑1 (MQ=255)
                                                  aaaGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATg                                             >  2:2470508/1‑100 (MQ=255)
                                                      cGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTc                                         <  2:3411007/100‑1 (MQ=255)
                                                      cGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTc                                         >  2:122594/1‑100 (MQ=255)
                                                       gTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCg                                        <  2:1380209/100‑1 (MQ=255)
                                                        ttattaACCCATGACGGTGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCgg                                       <  1:1225919/100‑1 (MQ=255)
                                                         tattaACCCATGACGATGTGAATCAACGCGCGTGGCAGGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGa                                      <  2:2810417/100‑1 (MQ=255)
                                                         tattaACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGa                                      <  1:992064/100‑1 (MQ=255)
                                                           ttaACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATg                                    <  1:618825/100‑1 (MQ=255)
                                                           ttaACCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATg                                    <  1:2377213/100‑1 (MQ=255)
                                                             aaCCCATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCa                                  >  2:1775220/1‑100 (MQ=255)
                                                             aaCCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCa                                  <  1:3757057/100‑1 (MQ=255)
                                                               cccATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAAc                                >  1:3708533/1‑100 (MQ=255)
                                                               cccATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGAttgccatgcgcgttgtttcacatcgtcatgggttaataacgctttgcgctcg                                <  1:3380042‑M2/100‑53 (MQ=255)
                                                               cccATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAAc                                >  2:887122/1‑100 (MQ=255)
                                                                ccATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATAGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACg                               >  2:1877022/1‑100 (MQ=255)
                                                                 cATGACGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGc                              >  1:1912493/1‑100 (MQ=255)
                                                                     aCGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAAc                          <  2:926213/100‑1 (MQ=255)
                                                                     aCGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAAc                          >  1:1468612/1‑100 (MQ=255)
                                                                      cGATGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTCTCGGATGCAACGGTAACt                         >  2:2949050/1‑100 (MQ=255)
                                                                        aTGTGAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTgc                       <  1:2543547/100‑1 (MQ=255)
                                                                         tgtgCAACAACGCGACTGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCACGAGCGTACAAGTGTCGATGTATCGGATGGGACGCTAACTg                        >  2:632238/1‑98 (MQ=255)
                                                                         tgtgAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTgcg                      <  1:3740009/100‑1 (MQ=255)
                                                                         tgtgAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTgcg                      >  1:2984905/1‑100 (MQ=255)
                                                                         tgtgAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATTTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTgcn                      <  2:2268982/100‑2 (MQ=255)
                                                                         tgtgAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTgcg                      >  2:3103433/1‑100 (MQ=255)
                                                                          gtgAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTgcgc                     >  2:568791/1‑100 (MQ=255)
                                                                          gtgAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTgcgc                     <  1:1979677/100‑1 (MQ=255)
                                                                           tgAAACAACGCGCTTGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCt                    >  1:2835337/1‑100 (MQ=255)
                                                                           tgAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTTCCCATGTCGATGTCTCGGATGCAACGCTAACTGCGCt                    >  2:4156475/1‑100 (MQ=255)
                                                                           tgAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCt                    >  2:3353932/1‑100 (MQ=255)
                                                                            gAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCt                    >  2:3999811/1‑99 (MQ=255)
                                                                            gAAACAACGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTg                   >  2:3283389/1‑100 (MQ=255)
                                                                            gAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATTTCGATGTTTCGGATGCACCGCTAACTgc                       >  2:1507640/1‑96 (MQ=255)
                                                                            gAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGAATGTGCGTACCCATGACGATGTTTCGGATGCAAAGCTAACTGCGCt                    >  1:4100783/1‑99 (MQ=255)
                                                                              aacaacGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGaa                 <  2:1205924/100‑1 (MQ=255)
                                                                              aacaacGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGAGGTTTCGGATGCAACGCTAACTGCGCTGaa                 >  1:3098779/1‑100 (MQ=255)
                                                                                 aacGCGCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGc              >  2:1722360/1‑100 (MQ=255)
                                                                                   cgcgCATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCaa            <  1:1931358/100‑1 (MQ=255)
                                                                                    gcgcATGGCAAGCGCTGAAATGGCAGAttgccatgcgcgttgtttcacatcgtcatgggttaataacgctttgcgctcggcccagcgttcaatgccttca           <  1:1163326‑M2/100‑74 (MQ=255)
                                                                                    gcgcATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCACGTGCGGACCCGGGTCGATGGTTGGGATGCAGCGCTAACTGCGCGGAAAGCaa            >  2:3477533/1‑99 (MQ=255)
                                                                                     cgcATGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATg          >  1:2220762/1‑100 (MQ=255)
                                                                                      gcATGGCAAGCGCTGAAATGGCAGAttgccatgcgcgttgtttcacatcgtcatgggttaataacgctttgcgctcggcccagcgttcaatgccttcaaa         <  1:462124‑M2/100‑76 (MQ=255)
                                                                                      gcATGGCAAGCGCTGAAATGGCAGAttgccatgcgcgttgtttcacatcgtcatgggttaataacgctttgcgctcggcccagcgttcaatgccttcaaa         <  1:2580799‑M2/100‑76 (MQ=255)
                                                                                         tGGCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTgg      <  2:4129037/100‑1 (MQ=255)
                                                                                          ggCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGa     >  1:3230773/1‑100 (MQ=255)
                                                                                          ggCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGa     <  2:806145/100‑1 (MQ=255)
                                                                                          ggCAAGCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGa     <  2:1158907/100‑1 (MQ=255)
                                                                                          ggCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGa     <  2:2012546/100‑1 (MQ=255)
                                                                                             aaaCGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGt  <  2:4105958/100‑1 (MQ=255)
                                                                                               |                                                                                                 
GGAATCAGTCCGGCACGCTGTTTGAAGGCATTGAACGCTGGGCCGAGCGCAAAGCGTTATTAACCCATGACGATGTGAAACAACGCGCATGGCAAACGCTGAAATGGCAGATTGCCAACGGCATTCAGCATGTGCGTACCCATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGT  >  NC_000913/356397‑356589

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: