Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,597,072 T→C D73G (GAC→GGC)  b3460 ← LivJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,597,0720TC100.0% 800.5 / NA 245D73G (GAC→GGC) b3460LivJ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (107/138);  total (107/138)

Only 100 of 250 total aligned reads displayed.

GAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGTCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGG  >  NC_000913/3596982‑3597161
                                                                                          |                                                                                         
 aagaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCTGGTCACAGGCATCGCCat                                                                                        >  2:4461969/1‑93 (MQ=255)
 aagaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCat                                                                                        <  2:2338902/93‑1 (MQ=255)
  agaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCata                                                                                       >  1:1211421/1‑93 (MQ=255)
  agaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCata                                                                                       <  2:1518509/93‑1 (MQ=255)
  agaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCata                                                                                       <  1:5362611/93‑1 (MQ=255)
   gaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACGTCGTGTTTCGGGTCACAGGCATCGCCatat                                                                                      <  2:1317591/93‑1 (MQ=255)
    aaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAtt                                                                                     <  1:3421423/93‑1 (MQ=255)
     aCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAttt                                                                                    >  1:3963366/1‑93 (MQ=255)
       agagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTa                                                                                  <  1:738777/93‑1 (MQ=255)
       agagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTa                                                                                  >  2:658219/1‑93 (MQ=255)
       agagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTa                                                                                  >  2:1998821/1‑93 (MQ=255)
        gagTTGACCAATCACATATTTACTGCCGTCGTTAACGTCTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTAc                                                                                 <  1:2049280/93‑1 (MQ=255)
        gagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTAc                                                                                 >  1:1717461/1‑93 (MQ=255)
        aagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTAc                                                                                 <  2:4885650/92‑1 (MQ=255)
         agGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACg                                                                                >  1:5154980/1‑93 (MQ=255)
         agGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACg                                                                                >  1:6160693/1‑93 (MQ=255)
         agGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACg                                                                                >  1:1948251/1‑93 (MQ=255)
          ggTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGa                                                                               <  2:1174829/93‑1 (MQ=255)
          ggTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGa                                                                               >  1:1394779/1‑93 (MQ=255)
            tGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGAtt                                                                             <  1:5566265/93‑1 (MQ=255)
              aCCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTg                                                                           <  1:3328329/93‑1 (MQ=255)
              aCCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTg                                                                           >  1:4869009/1‑93 (MQ=255)
                 aaTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAg                                                                        >  2:6043893/1‑93 (MQ=255)
                 aaTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAg                                                                        >  2:2286340/1‑93 (MQ=255)
                 aaTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAg                                                                        >  2:4602333/1‑93 (MQ=255)
                 aaTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAg                                                                        >  2:5802762/1‑93 (MQ=255)
                 aaTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAg                                                                        >  1:4545118/1‑93 (MQ=255)
                  aTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGt                                                                       >  1:2982964/1‑93 (MQ=255)
                  aTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGt                                                                       <  1:3829903/93‑1 (MQ=255)
                  aTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGt                                                                       >  1:247682/1‑93 (MQ=255)
                    cacaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGttt                                                                     <  1:2336496/93‑1 (MQ=255)
                     acaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttg                                                                    <  1:4511990/93‑1 (MQ=255)
                      caTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACGGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgt                                                                   <  2:2785612/93‑1 (MQ=255)
                      caTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgt                                                                   >  1:4093078/1‑93 (MQ=255)
                         atTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGc                                                                >  1:1571189/1‑93 (MQ=255)
                           ttAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCt                                                              <  2:2762606/93‑1 (MQ=255)
                           ttAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCt                                                              <  2:541041/93‑1 (MQ=255)
                            tAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCtt                                                             <  2:324850/93‑1 (MQ=255)
                            tAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCtt                                                             >  2:2152926/1‑93 (MQ=255)
                             aaTGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCttt                                                            <  2:2524612/93‑1 (MQ=255)
                              aggCCGTCTTGAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTa                                                           <  2:4994328/91‑1 (MQ=255)
                               tGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCTCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTaa                                                          <  2:5632706/93‑1 (MQ=255)
                               tGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGTAGTTTGTTGCCTTTaa                                                          <  2:1464564/93‑1 (MQ=255)
                               tGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTaa                                                          <  2:2950008/93‑1 (MQ=255)
                               tGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTaa                                                          <  2:3376359/93‑1 (MQ=255)
                               tGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTaa                                                          <  2:5587824/93‑1 (MQ=255)
                                 ccGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATg                                                        <  1:2439448/93‑1 (MQ=255)
                                 ccGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATg                                                        <  1:2214168/93‑1 (MQ=255)
                                  cgtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgc                                                       >  2:192768/1‑93 (MQ=255)
                                   gtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgcc                                                      >  2:2242324/1‑93 (MQ=255)
                                    tcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgccg                                                     >  1:1530415/1‑93 (MQ=255)
                                       ttAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCt                                                  >  2:2061650/1‑93 (MQ=255)
                                       ttAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCt                                                  <  2:886643/93‑1 (MQ=255)
                                         aaCGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCttt                                                >  1:1012437/1‑93 (MQ=255)
                                         aaCGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCttt                                                >  2:801924/1‑93 (MQ=255)
                                         aaCGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCttt                                                >  2:5209270/1‑93 (MQ=255)
                                         aaCGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCttt                                                >  1:6109843/1‑93 (MQ=255)
                                         aaCGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCttt                                                <  1:5733594/93‑1 (MQ=255)
                                         aaCGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCttt                                                <  1:32286/93‑1 (MQ=255)
                                          aCGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTa                                               <  1:2386355/93‑1 (MQ=255)
                                          aCGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTa                                               <  2:3762196/93‑1 (MQ=255)
                                          aCGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTa                                               <  1:4286360/93‑1 (MQ=255)
                                           cGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAg                                              <  2:513649/93‑1 (MQ=255)
                                             aCTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCg                                            >  1:910787/1‑93 (MQ=255)
                                             aCTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCg                                            <  2:4163483/93‑1 (MQ=255)
                                              cTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGt                                           >  1:5970537/1‑93 (MQ=255)
                                               tttGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGtt                                          >  1:5097259/1‑93 (MQ=255)
                                                ttGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTg                                         <  1:973207/93‑1 (MQ=255)
                                                  gTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGat                                       <  2:4459436/93‑1 (MQ=255)
                                                   ttCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGata                                      >  1:399896/1‑93 (MQ=255)
                                                    tCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGatat                                     <  2:1738404/93‑1 (MQ=255)
                                                         aCCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCa                                >  2:68880/1‑93 (MQ=255)
                                                           cGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAAc                              <  2:5842926/93‑1 (MQ=255)
                                                           cGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAAc                              <  1:3756421/93‑1 (MQ=255)
                                                             cAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAAcag                            >  2:6197801/1‑91 (MQ=255)
                                                             cAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCg                            <  2:4753319/93‑1 (MQ=255)
                                                             cAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCg                            >  1:3012271/1‑93 (MQ=255)
                                                              aaCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGc                           <  1:624912/93‑1 (MQ=255)
                                                               aCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGcc                          >  2:4015471/1‑93 (MQ=255)
                                                               aCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGcc                          >  1:3122281/1‑93 (MQ=255)
                                                                ccGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCt                         >  2:3146722/1‑93 (MQ=255)
                                                                ccGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCt                         >  2:2866087/1‑93 (MQ=255)
                                                                    cTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctc                     <  1:416174/93‑1 (MQ=255)
                                                                    cTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctc                     <  1:4705290/93‑1 (MQ=255)
                                                                     tGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctct                    >  2:5295315/1‑93 (MQ=255)
                                                                      gTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTg                   <  1:547926/93‑1 (MQ=255)
                                                                      gTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTg                   <  1:2824714/93‑1 (MQ=255)
                                                                       tttCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTgc                  <  1:3565760/93‑1 (MQ=255)
                                                                       tttCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTgc                  <  1:3398759/93‑1 (MQ=255)
                                                                           gggTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCg              >  2:1602206/1‑93 (MQ=255)
                                                                              tCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTa           <  1:1807979/93‑1 (MQ=255)
                                                                               cacaGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTaa          >  2:2452020/1‑93 (MQ=255)
                                                                                 caGGCATCTCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATCCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGTTAAAc        <  2:1352313/93‑1 (MQ=255)
                                                                                 caGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAAc        <  1:525831/93‑1 (MQ=255)
                                                                                  aGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACt       <  1:3715973/93‑1 (MQ=255)
                                                                                  aGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACt       >  1:5207633/1‑93 (MQ=255)
                                                                                    gCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTcc     <  2:1103583/93‑1 (MQ=255)
                                                                                    gCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTcc     <  2:2358989/93‑1 (MQ=255)
                                                                                       tCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTgg  >  1:303796/1‑93 (MQ=255)
                                                                                       tCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTgg  >  1:5506419/1‑93 (MQ=255)
                                                                                          |                                                                                         
GAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGTCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGG  >  NC_000913/3596982‑3597161

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