Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,808,645 G→A S45S (AGC→AGT mutM ← formamidopyrimidine DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,808,6450GA85.1% 898.7 / 106.5 335S45S (AGC→AGTmutMformamidopyrimidine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (32/18);  new base A (162/123);  total (194/141)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

Only 100 of 336 total aligned reads displayed.

CCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTT  >  NC_000913/3808553‑3808736
                                                                                            |                                                                                           
    aCATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTa                                                                                         >  1:2316796/1‑93 (MQ=255)
    aCATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTa                                                                                         >  2:1018609/1‑93 (MQ=255)
      aTCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTaaa                                                                                       <  2:1234499/93‑1 (MQ=255)
      aTCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTaaa                                                                                       >  1:4589107/1‑93 (MQ=255)
      aTCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTaaa                                                                                       >  2:6170149/1‑93 (MQ=255)
       tCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGGCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAAc                                                                                      >  1:5781231/1‑93 (MQ=255)
        cccTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACg                                                                                     <  1:767661/93‑1 (MQ=255)
        cccTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACg                                                                                     <  2:2022031/93‑1 (MQ=255)
        cccTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACg                                                                                     <  1:5489966/93‑1 (MQ=255)
         ccTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACgg                                                                                    >  2:5785157/1‑93 (MQ=255)
         ccTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACgg                                                                                    >  2:4753465/1‑93 (MQ=255)
          cTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGt                                                                                   <  1:6347127/93‑1 (MQ=255)
          cTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGt                                                                                   <  2:4232196/93‑1 (MQ=255)
          cTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGt                                                                                   <  1:5581252/93‑1 (MQ=255)
           tAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTa                                                                                  >  2:2233290/1‑93 (MQ=255)
             aatgaatgATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGa                                                                                >  1:5275273/1‑93 (MQ=255)
             aatgaatgATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGa                                                                                <  1:6855528/93‑1 (MQ=255)
             aatgaatgATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGa                                                                                >  2:2724127/1‑93 (MQ=255)
             aatgaatgATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGa                                                                                >  1:7105336/1‑93 (MQ=255)
              atgaatgaTAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGAt                                                                               >  2:3390263/1‑93 (MQ=255)
               tgaatgaTAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGAtc                                                                              >  1:6012249/1‑93 (MQ=255)
               tgaatgaTAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGAtc                                                                              >  2:5907255/1‑93 (MQ=255)
               tgaatgaTAATCCAGCCCTCAGGCAGCTCAGGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGAtc                                                                              <  1:3906343/93‑1 (MQ=255)
                 aatgaTAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGAtctc                                                                            <  2:3797504/93‑1 (MQ=255)
                   tgaTAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCtctt                                                                          <  2:818962/93‑1 (MQ=255)
                    gaTAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCtcttc                                                                         >  2:1056904/1‑93 (MQ=255)
                    gaTAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCtcttc                                                                         <  2:4802383/93‑1 (MQ=255)
                    gaTAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCtcttc                                                                         <  1:3474487/93‑1 (MQ=255)
                     aTAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCtcttct                                                                        >  2:7413429/1‑93 (MQ=255)
                      tAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTg                                                                       >  2:6939493/1‑93 (MQ=255)
                      tAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTg                                                                       >  1:6417461/1‑93 (MQ=255)
                       aacccAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGa                                                                      <  1:5135443/90‑1 (MQ=255)
                        aTCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGaa                                                                     <  2:349607/93‑1 (MQ=255)
                        aTCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGaa                                                                     >  1:2234446/1‑93 (MQ=255)
                         tCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGaaa                                                                    <  2:60978/93‑1 (MQ=255)
                         tCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGaaa                                                                    >  1:1647278/1‑93 (MQ=255)
                          ccAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAAc                                                                   <  2:6699944/93‑1 (MQ=255)
                            aGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCg                                                                 >  2:1406645/1‑93 (MQ=255)
                            aGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCg                                                                 >  2:6484806/1‑93 (MQ=255)
                             gCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCgg                                                                <  1:2854839/93‑1 (MQ=255)
                              cccTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGc                                                               <  1:6450087/93‑1 (MQ=255)
                               ccTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGcc                                                              >  2:1681346/1‑93 (MQ=255)
                                cTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCa                                                             >  1:2600990/1‑93 (MQ=255)
                                   aGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAgcg                                                          >  2:5125369/1‑93 (MQ=255)
                                    ggCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAgcgc                                                         >  2:13075/1‑93 (MQ=255)
                                       aGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAg                                                      >  2:6598234/1‑93 (MQ=255)
                                        gCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGc                                                     >  2:2886/1‑93 (MQ=255)
                                        gCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGc                                                     >  2:4562569/1‑93 (MQ=255)
                                         cTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCg                                                    >  1:4813835/1‑93 (MQ=255)
                                           ccAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTc                                                  <  2:7329098/93‑1 (MQ=255)
                                           ccAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTc                                                  >  2:2541572/1‑93 (MQ=255)
                                              gcagcagATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGt                                               >  2:5345863/1‑93 (MQ=255)
                                              gcagcagATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGt                                               <  2:7262050/93‑1 (MQ=255)
                                                agcagATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTg                                             >  1:6726593/1‑93 (MQ=255)
                                                 gcagATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTgc                                            >  1:4653908/1‑93 (MQ=255)
                                                    gATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGca                                         >  2:3955180/1‑93 (MQ=255)
                                                    gATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGca                                         >  1:7252899/1‑93 (MQ=255)
                                                     atatTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGcac                                        >  1:1258893/1‑93 (MQ=255)
                                                      tatTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGcacc                                       >  1:6199139/1‑93 (MQ=255)
                                                      tatTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGcacc                                       <  2:5457228/93‑1 (MQ=255)
                                                        tttAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGcaccac                                     <  1:3057342/93‑1 (MQ=255)
                                                         ttAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACt                                    >  1:1005770/1‑93 (MQ=255)
                                                           aGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGc                                  <  1:4386871/93‑1 (MQ=255)
                                                            gCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCa                                 >  2:2419601/1‑93 (MQ=255)
                                                            gCCCGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCa                                 >  2:7103143/1‑93 (MQ=255)
                                                             cccGCCGCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCAt                                <  1:1401757/93‑1 (MQ=255)
                                                              ccgccgCTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATg                               <  1:6664813/93‑1 (MQ=255)
                                                                gccgcTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATgaa                             <  2:3662276/93‑1 (MQ=255)
                                                                 ccgcTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATgaag                            >  1:5937003/1‑93 (MQ=255)
                                                                  cgcTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATgaaga                           >  2:4115894/1‑93 (MQ=255)
                                                                   gcTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATgaagaa                          >  2:4321142/1‑93 (MQ=255)
                                                                   gcTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATgaagaa                          <  1:1169204/93‑1 (MQ=255)
                                                                   gcTGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATgaagaa                          >  1:2644211/1‑93 (MQ=255)
                                                                    cTGCACGCTAAGCACGGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAAt                         <  1:2376786/93‑1 (MQ=255)
                                                                     tGCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATg                        >  1:4406149/1‑93 (MQ=255)
                                                                      gCACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATgg                       >  2:147788/1‑93 (MQ=255)
                                                                       cACGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGt                      <  1:3968361/93‑1 (MQ=255)
                                                                        aCGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGtt                     >  2:2137481/1‑93 (MQ=255)
                                                                         cGCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTg                    <  1:4785133/93‑1 (MQ=255)
                                                                          gCTAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGc                   >  1:1947169/1‑93 (MQ=255)
                                                                            tAAGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCAc                 >  1:2529447/1‑93 (MQ=255)
                                                                             aaGCACTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCAcc                >  1:413641/1‑93 (MQ=255)
                                                                                 aCTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACg            >  1:4157103/1‑93 (MQ=255)
                                                                                  cTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACga           >  2:3128196/1‑93 (MQ=255)
                                                                                  cTGGTTGGTCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACga           >  2:4878098/1‑93 (MQ=255)
                                                                                   tggttggtCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACgag          <  2:6264070/93‑1 (MQ=255)
                                                                                   tggttggtCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACgag          >  2:6009846/1‑93 (MQ=255)
                                                                                   tggttggtCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACgag          >  1:3567235/1‑93 (MQ=255)
                                                                                    ggttggtCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACgaga         <  1:5019655/93‑1 (MQ=255)
                                                                                    ggttggtCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACgaga         <  2:5589825/93‑1 (MQ=255)
                                                                                    ggttggtCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACgaga         <  1:5446029/93‑1 (MQ=255)
                                                                                    ggttggtCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACgaga         <  1:607159/93‑1 (MQ=255)
                                                                                    ggttggtCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACgaga         <  2:7085552/93‑1 (MQ=255)
                                                                                    ggttggtCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACgaga         >  2:1708617/1‑93 (MQ=255)
                                                                                     gttggtCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGAt        >  1:3321431/1‑93 (MQ=255)
                                                                                       tggtCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGc      >  1:3788509/1‑93 (MQ=255)
                                                                                        ggtCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCg     <  2:3265938/93‑1 (MQ=255)
                                                                                         gtCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCgg    >  2:4988177/1‑93 (MQ=255)
                                                                                         gtCACTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCgg    <  2:2680619/93‑1 (MQ=255)
                                                                                           cGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGtt  <  2:878197/93‑1 (MQ=255)
                                                                                            |                                                                                           
CCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTT  >  NC_000913/3808553‑3808736

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