Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,464,165 T→C E66E (GAA→GAG bfr ← bacterioferrin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,464,1650TC87.8% 771.1 / 90.0 305E66E (GAA→GAGbfrbacterioferrin
Reads supporting (aligned to +/- strand):  ref base T (16/21);  new base C (124/143);  total (140/165)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 306 total aligned reads displayed.

CCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGA  >  NC_000913/3464073‑3464255
                                                                                            |                                                                                          
ccATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCt                                                                                            <  1:7477312/93‑1 (MQ=255)
   tCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTcc                                                                                         >  2:804931/1‑93 (MQ=255)
   tCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTcc                                                                                         >  2:5286392/1‑93 (MQ=255)
    ccAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCa                                                                                        <  1:3937692/93‑1 (MQ=255)
       gCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGaa                                                                                     <  2:4199756/93‑1 (MQ=255)
       gCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGaa                                                                                     >  2:369450/1‑93 (MQ=255)
       gCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGaa                                                                                     <  1:522128/93‑1 (MQ=255)
        cTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGaaa                                                                                    >  1:1382097/1‑93 (MQ=255)
            aGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGa                                                                                <  1:5546493/93‑1 (MQ=255)
             gTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGaa                                                                               <  1:7297521/93‑1 (MQ=255)
              tGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAAt                                                                              >  1:2383166/1‑93 (MQ=255)
               gCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATg                                                                             <  1:1566051/93‑1 (MQ=255)
                ccAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATgc                                                                            <  1:1173456/93‑1 (MQ=255)
                ccAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATgc                                                                            <  1:4899934/93‑1 (MQ=255)
                  aGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATgcgc                                                                          <  2:5775518/93‑1 (MQ=255)
                   gATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCt                                                                         <  2:3511796/93‑1 (MQ=255)
                     tCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCa                                                                       >  2:7418576/1‑93 (MQ=255)
                      cAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCaa                                                                      <  2:5314633/93‑1 (MQ=255)
                       aGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAat                                                                     <  2:4292821/93‑1 (MQ=255)
                        gAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAata                                                                    <  2:2859003/93‑1 (MQ=255)
                         aaCGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAatat                                                                   >  2:1443022/1‑93 (MQ=255)
                          aCGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAAtata                                                                  >  2:4353934/1‑93 (MQ=255)
                           cGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATaa                                                                 <  2:2707282/93‑1 (MQ=255)
                           cGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATaa                                                                 >  1:7404304/1‑93 (MQ=255)
                             cagcaTTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACg                                                               >  2:1351698/1‑93 (MQ=255)
                              agcaTTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGa                                                              <  1:6744584/93‑1 (MQ=255)
                               gcaTTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGAt                                                             <  2:2941053/93‑1 (MQ=255)
                               gcaTTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGAt                                                             >  1:27447/1‑93 (MQ=255)
                                caTTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATc                                                            <  1:5141834/93‑1 (MQ=255)
                                caTTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATc                                                            >  2:7295881/1‑93 (MQ=255)
                                caTTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATc                                                            <  2:1033929/93‑1 (MQ=255)
                                  tttCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCgg                                                          >  2:1522222/1‑93 (MQ=255)
                                  tttCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCgg                                                          >  1:6512571/1‑93 (MQ=255)
                                   ttCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGc                                                         <  2:6525141/93‑1 (MQ=255)
                                    tCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCg                                                        >  1:3264916/1‑93 (MQ=255)
                                     ccTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCgt                                                       >  1:2237843/1‑93 (MQ=255)
                                     ccTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCgt                                                       >  2:6712117/1‑93 (MQ=255)
                                      cTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCgtg                                                      >  1:5433990/1‑93 (MQ=255)
                                       tCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCgtgt                                                     <  2:4036192/93‑1 (MQ=255)
                                       tCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCgtgt                                                     >  2:5142973/1‑93 (MQ=255)
                                        cAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGtt                                                    <  1:3024297/93‑1 (MQ=255)
                                          aCATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtc                                                  >  1:4319149/1‑93 (MQ=255)
                                           cATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtca                                                 >  1:4541269/1‑93 (MQ=255)
                                           cATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtca                                                 <  1:1874262/93‑1 (MQ=255)
                                           cATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtca                                                 <  1:342206/93‑1 (MQ=255)
                                           cATCTTCACCAATGTTCAGTTTGCCCAGGTCCTATAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtca                                                 <  2:6895975/93‑1 (MQ=255)
                                             tCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtcatc                                               >  2:2100209/1‑93 (MQ=255)
                                              cTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtcatct                                              <  1:6080897/93‑1 (MQ=255)
                                                tCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtcatctca                                            <  2:663653/93‑1 (MQ=255)
                                                tCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtcatctca                                            <  2:3163811/93‑1 (MQ=255)
                                                   ccAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCtcatca                                         >  2:5702071/1‑93 (MQ=255)
                                                    cAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCaa                                        >  1:176187/1‑93 (MQ=255)
                                                     aaTGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAAt                                       >  2:1392902/1‑93 (MQ=255)
                                                      aTGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATg                                      >  1:403521/1‑93 (MQ=255)
                                                      aTGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATg                                      <  1:5814547/93‑1 (MQ=255)
                                                      aTGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATg                                      >  1:4828186/1‑93 (MQ=255)
                                                      aTGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATg                                      <  1:952622/93‑1 (MQ=255)
                                                       tGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATgg                                     >  1:7205022/1‑93 (MQ=255)
                                                        gTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGa                                    >  1:81064/1‑93 (MQ=255)
                                                          tCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGAtt                                  >  2:5440879/1‑93 (MQ=255)
                                                          tCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGAtt                                  <  2:328647/93‑1 (MQ=255)
                                                            aGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCa                                <  2:7070700/93‑1 (MQ=255)
                                                            aGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCa                                <  2:2537020/93‑1 (MQ=255)
                                                            aGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCa                                >  2:6140370/1‑93 (MQ=255)
                                                              tttGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATg                              >  2:1804797/1‑93 (MQ=255)
                                                                tGCCCAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGAt                            <  2:1911695/93‑1 (MQ=255)
                                                                    cAGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTc                        >  1:335880/1‑93 (MQ=255)
                                                                     aGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTcc                       >  2:3096951/1‑93 (MQ=255)
                                                                     aGGTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTcc                       >  2:1163600/1‑93 (MQ=255)
                                                                      ggTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCca                      >  1:6585232/1‑93 (MQ=255)
                                                                      ggTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCca                      <  2:6348471/93‑1 (MQ=255)
                                                                       gTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCcac                     >  2:6045722/1‑93 (MQ=255)
                                                                       gTCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCcac                     <  1:1814676/93‑1 (MQ=255)
                                                                        tCCTGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCcaca                    <  1:6325393/93‑1 (MQ=255)
                                                                           tGTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCAcatca                 <  1:6846247/93‑1 (MQ=255)
                                                                            gTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCAcatcat                >  1:6673204/1‑93 (MQ=255)
                                                                            gTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCAcatcat                <  2:4668735/93‑1 (MQ=255)
                                                                            gTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCAcatcat                >  1:5100792/1‑93 (MQ=255)
                                                                            gTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCAcatcat                >  1:711642/1‑93 (MQ=255)
                                                                            gTAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCAcatcat                >  1:3239708/1‑93 (MQ=255)
                                                                             tAAGTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCAtt               <  1:3685572/93‑1 (MQ=255)
                                                                                gTTTGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTgag            >  1:7052574/1‑93 (MQ=255)
                                                                                 tttGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTgaga           <  1:1623978/93‑1 (MQ=255)
                                                                                  ttGGAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGAc          <  1:2517042/93‑1 (MQ=255)
                                                                                    ggAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGt        <  2:4539678/93‑1 (MQ=255)
                                                                                    ggAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGt        <  2:3430396/93‑1 (MQ=255)
                                                                                     gAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGtt       <  1:6781917/93‑1 (MQ=255)
                                                                                     gAAGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGtt       >  2:6213519/1‑93 (MQ=255)
                                                                                      aaGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGttt      <  1:4650980/93‑1 (MQ=255)
                                                                                      aaGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGttt      >  2:4827481/1‑93 (MQ=255)
                                                                                      aaGACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGttt      <  1:4112268/93‑1 (MQ=255)
                                                                                       aGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTg     <  2:3817072/93‑1 (MQ=255)
                                                                                        gACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTga    >  1:2068452/1‑93 (MQ=255)
                                                                                        gACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTga    >  2:2387181/1‑93 (MQ=255)
                                                                                        gACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTga    <  1:6363209/93‑1 (MQ=255)
                                                                                        gACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTga    <  2:5916503/93‑1 (MQ=255)
                                                                                        gACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTga    >  2:3955100/1‑93 (MQ=255)
                                                                                        gACCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTga    >  1:4220616/1‑93 (MQ=255)
                                                                                         aCCCTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTgag   >  2:5235194/1‑93 (MQ=255)
                                                                                          ccTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTgaga  <  1:5006042/93‑1 (MQ=255)
                                                                                            |                                                                                          
CCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGA  >  NC_000913/3464073‑3464255

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