Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,926,042 G→T 17.1% C232F (TGC→TTC)  queF → 7‑cyano‑7‑deazaguanine reductase (NADPH‑dependent)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,926,0420GT17.1% 260.5 / 40.0 103C232F (TGC→TTC) queF7‑cyano‑7‑deazaguanine reductase (NADPH‑dependent)
Reads supporting (aligned to +/- strand):  ref base G (42/43);  new base T (11/7);  total (53/50)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01

Only 100 of 103 total aligned reads displayed.

GCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAAC  >  NC_000913/2925947‑2926139
                                                                                               |                                                                                                 
gCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTg                                                                                               >  2:2145090/1‑100 (MQ=255)
gCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTg                                                                                               >  1:1678943/1‑100 (MQ=255)
gCTCCAAATTCAGTATCGTGGACGCCAAATTAACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTg                                                                                               <  2:1933132/100‑1 (MQ=255)
 cTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTgg                                                                                              >  2:2308027/1‑100 (MQ=255)
  tCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGa                                                                                             <  1:2164102/100‑1 (MQ=255)
   ccAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGaa                                                                                            >  2:621088/1‑100 (MQ=255)
    cAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGCACTGTTCGTGGAAc                                                                                           <  1:536400/100‑1 (MQ=255)
     aaaTTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCCTTCCGTCATCACAACGAGTTCCACGCACTGTTCGTGGAAct                                                                                          <  2:1303438/100‑2 (MQ=255)
     aaaTTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACg                                                                                          <  1:279464/100‑1 (MQ=255)
      aaTTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGCACTGTTCGTGGAActc                                                                                         <  1:21851/100‑3 (MQ=255)
      aaTTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGCACTGTTCGTGGAActc                                                                                         <  1:2150109/100‑3 (MQ=255)
      aaTTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGc                                                                                         >  2:1924378/1‑100 (MQ=255)
       aTTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGCACTGTTCGTGGAActcg                                                                                        <  1:782016/100‑4 (MQ=38)
       aTTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCa                                                                                        <  2:1141340/100‑1 (MQ=255)
        ttCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGCACTGTTCGTGGAActcgt                                                                                       <  1:1008788/100‑5 (MQ=38)
          cAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGCACTGTTCGTGGAActcgttg                                                                                     <  1:1746666/100‑7 (MQ=38)
          cAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAAAGCATc                                                                                      >  1:1843617/1‑99 (MQ=255)
           aGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCtt                                                                                    <  1:1094129/100‑1 (MQ=255)
           aGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCtt                                                                                    <  1:1077001/100‑1 (MQ=255)
            gTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCTTCATCACAACGAGATCCACGAACAGTGCGTGGAACGCATCttt                                                                                   >  2:429234/1‑100 (MQ=255)
            gTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCttt                                                                                   >  1:1752144/1‑100 (MQ=255)
            gTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCttt                                                                                   <  2:1410927/100‑1 (MQ=255)
                cGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATg                                                                               <  2:1370483/100‑1 (MQ=255)
                 gTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGa                                                                              >  1:1708625/1‑100 (MQ=255)
                      cGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGt                                                                         >  2:566360/1‑100 (MQ=255)
                       gCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGtt                                                                        >  1:1146527/1‑100 (MQ=255)
                        ccAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTa                                                                       >  1:2296034/1‑100 (MQ=255)
                        ccAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTa                                                                       >  1:891325/1‑100 (MQ=255)
                        ccAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTa                                                                       >  1:994279/1‑100 (MQ=255)
                         cAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCATCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTAc                                                                      >  2:1604561/1‑100 (MQ=255)
                          aaaTTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACg                                                                     <  2:1918851/100‑1 (MQ=255)
                          aaaTTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGGTCGTGGAACGCATCTTTAATGACCTGTTACg                                                                     >  1:1914647/1‑100 (MQ=255)
                           aaTTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGc                                                                    <  2:1249734/100‑1 (MQ=255)
                               gACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCt                                                                >  1:168844/1‑100 (MQ=255)
                                  agagAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgcc                                                             <  1:627870/100‑1 (MQ=255)
                                   gggaagaaCTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgcca                                                            <  2:714738/98‑1 (MQ=255)
                                   gagaAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgcca                                                            <  1:891621/100‑1 (MQ=255)
                                   gagaAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgcca                                                            <  2:772107/100‑1 (MQ=255)
                                   gagaAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgcca                                                            >  2:331891/1‑100 (MQ=255)
                                   gagaAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgcca                                                            >  2:2052856/1‑100 (MQ=255)
                                    agaAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgccag                                                           >  1:443925/1‑100 (MQ=255)
                                     gaAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgccagc                                                          >  1:1339647/1‑100 (MQ=255)
                                     gaAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgcaagc                                                          >  1:633456/1‑100 (MQ=255)
                                      aaaaaCTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgccagcc                                                         >  2:903939/1‑100 (MQ=255)
                                       aaaaCTGCTGTGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgccagcca                                                        >  1:17850/1‑100 (MQ=255)
                                           ctgctgCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGaaa                                                    >  2:490482/1‑100 (MQ=255)
                                                gcGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTTCGTGTGACGCATCTTTAATGACCGGTTACGCTTCTGCCAGCCAGAAAAAtt                                                >  2:1395074/1‑99 (MQ=255)
                                                gcGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTg                                               <  1:1527095/100‑1 (MQ=255)
                                                gcGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTg                                               <  1:490485/100‑1 (MQ=255)
                                                  gTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCTAGCCAGAAAAAGTGAg                                             <  2:1235105/100‑1 (MQ=255)
                                                  gTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAg                                             <  1:2212611/100‑1 (MQ=255)
                                                        tGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTa                                       <  1:574023/100‑1 (MQ=255)
                                                        tGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTa                                       >  1:1223548/1‑100 (MQ=255)
                                                           tttcATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTAc                                      <  1:641692/96‑1 (MQ=255)
                                                           tctcATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGc                                    >  1:2651/1‑100 (MQ=255)
                                                           tctcATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGc                                    >  1:1733688/1‑100 (MQ=255)
                                                            ctcATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCa                                   >  1:2195326/1‑100 (MQ=255)
                                                              cATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACg                                 >  2:1336534/1‑100 (MQ=255)
                                                                 tCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTta                              <  1:1057927/100‑1 (MQ=255)
                                                                 tCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTta                              <  2:462848/100‑1 (MQ=255)
                                                                 tCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTta                              <  1:2168342/100‑1 (MQ=255)
                                                                  ccGTCATAACAACGAGTTCCACGAAGAGTGCGTGGAACGAATCTTTAATGACCTGTTACGCTTATGCAAGACAGAAAAATTGAGCGTTTACGCACGtt                               <  2:2251881/98‑1 (MQ=255)
                                                                   cGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTtata                            >  1:634139/1‑100 (MQ=255)
                                                                     tcatcaCAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATAcc                          <  1:203052/100‑1 (MQ=255)
                                                                     tcatcaCAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATAcc                          <  1:452961/100‑1 (MQ=255)
                                                                     tcatcaCAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATAcc                          <  2:1519023/100‑1 (MQ=255)
                                                                       atcaCAACGAGTTCCACGAACAGTGCGTGCAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCcg                        <  2:1439571/100‑1 (MQ=255)
                                                                       atcaCAACGAGTTCCACGAACAGTGCATAGAACGCAACTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTTAGAGTTTACGCACGTTATAcc                          <  2:1504123/98‑1 (MQ=255)
                                                                         cacaACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCcgtc                      <  2:1992404/100‑1 (MQ=255)
                                                                         cacaACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTGACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCcgtc                      >  2:1067746/1‑100 (MQ=255)
                                                                            gatgcgTTCCACGCACTGTTCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTg                   >  1:1001026/6‑100 (MQ=38)
                                                                            gatgcgTTCCACGCACTGTTCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTg                   >  1:1505930/6‑100 (MQ=38)
                                                                             ctgcgTTCCACGCACTGTTCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTgg                  >  1:572333/5‑100 (MQ=38)
                                                                             atgcgTTCCACGCACTGTTTGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTgg                  >  2:1494153/5‑100 (MQ=255)
                                                                             atgcgTTCCACGCACTGTTCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTgg                  >  1:960260/5‑100 (MQ=38)
                                                                             atgcgTTCCACGCACTGTTCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTgg                  >  1:1086251/5‑100 (MQ=38)
                                                                             atgcgTTCCACGCACTGTTCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCAAGTTATACCCGTCATgg                  >  2:2018082/5‑100 (MQ=38)
                                                                             aCGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTgg                  >  1:794820/1‑100 (MQ=255)
                                                                              cGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGATTACGCACGTTATACCCGTCGTgg                  >  1:2219747/1‑99 (MQ=255)
                                                                              cGAGTGCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCTAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGc                 <  2:1663018/100‑1 (MQ=255)
                                                                                 gTTCCACGCACTGTTCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGt              >  2:1709909/1‑100 (MQ=255)
                                                                                  ttCCACGAACAGTGCGTGTAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTc             <  1:1395088/100‑1 (MQ=255)
                                                                                  ttCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTc             <  2:922058/100‑1 (MQ=255)
                                                                                  ttCCACGAACAGTGCGTGGAACGCATCTTTAATGAACTGTTACGCTTCTGCCAGCCAGAAAAATTTATCGTTTACTTACGTTATACTGGTCGTGGCg                <  2:526164/97‑1 (MQ=255)
                                                                                    ccACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTg           <  2:1375499/100‑1 (MQ=255)
                                                                                    ccACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTg           <  1:1228804/100‑1 (MQ=255)
                                                                                     cACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTgg          >  1:7249/1‑100 (MQ=255)
                                                                                     cACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTgg          >  1:656465/1‑100 (MQ=255)
                                                                                     cACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTgg          >  2:2290098/1‑100 (MQ=255)
                                                                                     cACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTgg          <  1:1336548/100‑1 (MQ=255)
                                                                                     cACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTgg          >  2:1536526/1‑100 (MQ=255)
                                                                                     cACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCCTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTAGCGGTCTg           >  1:866127/1‑99 (MQ=255)
                                                                                       cGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGAc        <  2:2090110/100‑1 (MQ=255)
                                                                                       cGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGAc        <  2:600430/100‑1 (MQ=255)
                                                                                         aaCAGTGCGTGGAACGCATCTTTAATGACCTTTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACAt      <  1:1782097/100‑1 (MQ=255)
                                                                                         aaCAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACa       <  1:750203/99‑1 (MQ=255)
                                                                                         aaCAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACAt      >  2:773236/1‑100 (MQ=255)
                                                                                           cAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTa    >  1:1125511/1‑100 (MQ=255)
                                                                                            aGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTaa   <  2:2403095/100‑1 (MQ=255)
                                                                                             gTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAAc  >  2:647290/1‑100 (MQ=255)
                                                                                               |                                                                                                 
GCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAAC  >  NC_000913/2925947‑2926139

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: