Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,632,444 A→G S146P (TCC→CCC)  rob ← right origin‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,632,4440AG98.9% 542.3 / ‑2.5 178S146P (TCC→CCC) robright origin‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (90/84);  total (90/88)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-02

Only 100 of 182 total aligned reads displayed.

CGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCG  >  NC_000913/4632355‑4632533
                                                                                         |                                                                                         
 ggAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAg                                                                                       >  1:3339202/1‑93 (MQ=255)
 ggAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAg                                                                                       <  2:444998/93‑1 (MQ=255)
   aaTGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGct                                                                                     >  1:2300398/1‑93 (MQ=255)
    aTGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGctc                                                                                    <  1:3344052/93‑1 (MQ=255)
    aTGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGctc                                                                                    >  2:3171248/1‑93 (MQ=255)
     tGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGctct                                                                                   >  1:4326690/1‑93 (MQ=255)
       gTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTgg                                                                                 >  1:773756/1‑93 (MQ=255)
       gTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGGAGCTCTgg                                                                                 >  2:3539977/1‑93 (MQ=255)
          ggCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTa                                                                              >  1:2503022/1‑93 (MQ=255)
            cgcgTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAac                                                                            <  2:1732287/93‑1 (MQ=255)
             gcgTTGCCGAGAAAATCGTGCCAGAGCTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAaca                                                                           >  1:339656/1‑93 (MQ=255)
              cgTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAacac                                                                          <  2:1465836/93‑1 (MQ=255)
              cgTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAacac                                                                          >  2:3246964/1‑93 (MQ=255)
              cgTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAacac                                                                          <  1:2750642/93‑1 (MQ=255)
                ttGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCa                                                                        >  1:3045300/1‑93 (MQ=255)
                ttGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCGCCAGCGAACAGGGGTAGCTCTGGGTAACACCa                                                                        <  2:1423813/93‑1 (MQ=255)
                  gCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAAt                                                                      <  1:2073484/93‑1 (MQ=255)
                  gCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAAt                                                                      >  2:4190519/1‑93 (MQ=255)
                   ccGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATc                                                                     <  2:789522/93‑1 (MQ=255)
                    cGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCa                                                                    <  1:648864/93‑1 (MQ=255)
                     gagaAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAg                                                                   <  2:596631/93‑1 (MQ=255)
                        aaaaTCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAgcgg                                                                >  1:2847269/1‑93 (MQ=255)
                         aaaTCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCACTCAgcggc                                                               >  1:1747755/1‑93 (MQ=255)
                          aaTCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAgcggcg                                                              <  1:2936675/93‑1 (MQ=255)
                          aaTCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAgcggcg                                                              >  2:2236174/1‑93 (MQ=255)
                            tCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGTGTAGCTCTGGGTAACACCAATCAGCGGCGTa                                                            <  2:2007231/93‑1 (MQ=255)
                            tCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTa                                                            <  1:3636909/93‑1 (MQ=255)
                            tCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTa                                                            >  2:3155086/1‑93 (MQ=255)
                            tCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTa                                                            <  2:1064325/93‑1 (MQ=255)
                            tCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTa                                                            <  1:2543552/93‑1 (MQ=255)
                            tCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTa                                                            >  2:533605/1‑93 (MQ=255)
                              gTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATc                                                          <  1:3599110/93‑1 (MQ=255)
                                 ccAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTc                                                       <  2:395286/93‑1 (MQ=255)
                                  cAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTcc                                                      >  1:169262/1‑93 (MQ=255)
                                   aGAACTGATAACGCATTTCATGGCTGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCa                                                     <  1:2279798/93‑1 (MQ=255)
                                   aGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGGGGCGTATCTTccc                                                     >  2:1539406/1‑92 (MQ=255)
                                    gAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAg                                                    >  2:271254/1‑93 (MQ=255)
                                    gAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAg                                                    <  1:3670046/93‑1 (MQ=255)
                                    gAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAg                                                    >  1:2239371/1‑93 (MQ=255)
                                    gAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAg                                                    <  1:3732129/93‑1 (MQ=255)
                                     aaCTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAgg                                                   <  2:2169145/93‑1 (MQ=255)
                                     aaCTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAgg                                                   >  1:388891/1‑93 (MQ=255)
                                      aCTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAggg                                                  >  2:3646460/1‑93 (MQ=255)
                                       cTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGt                                                 <  2:254151/93‑1 (MQ=255)
                                       cTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGt                                                 >  2:4641160/1‑93 (MQ=255)
                                        tGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTg                                                >  1:110834/1‑93 (MQ=255)
                                         gATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGa                                               >  1:1566325/1‑93 (MQ=255)
                                         gATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGa                                               >  2:4485144/1‑93 (MQ=255)
                                         gATAACGCATTTCACGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGa                                               <  1:4061212/93‑1 (MQ=255)
                                          aTAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGAc                                              >  2:1722115/1‑93 (MQ=255)
                                             aCGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACaaa                                           >  1:4150156/1‑93 (MQ=255)
                                             aCGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACaaa                                           >  1:3713383/1‑93 (MQ=255)
                                              cGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAAt                                          >  1:1565639/1‑93 (MQ=255)
                                                   ttCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTtgtg                                     >  1:1455792/1‑93 (MQ=255)
                                                    tCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGc                                    >  1:192751/1‑93 (MQ=255)
                                                      aTGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGctc                                  <  2:4299482/93‑1 (MQ=255)
                                                      aTGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGctc                                  <  2:3648090/93‑1 (MQ=255)
                                                       tGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGctct                                 >  2:2599424/1‑93 (MQ=255)
                                                        ggCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTg                                <  2:815586/93‑1 (MQ=255)
                                                        ggCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTg                                <  1:398203/93‑1 (MQ=255)
                                                         gCGGAAATCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTgg                               <  2:1793549/93‑1 (MQ=255)
                                                            gAAATCAGAGATTTGCTCCAGCGAACAGGGGGAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCAt                            >  2:3143370/1‑93 (MQ=255)
                                                             aaaTCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATa                           >  2:677348/1‑93 (MQ=255)
                                                             aaaTCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATa                           <  1:2143547/93‑1 (MQ=255)
                                                             aaaTCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATa                           <  1:2697477/93‑1 (MQ=255)
                                                             aaaTCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATa                           <  2:3422511/93‑1 (MQ=255)
                                                              aaTCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAg                          >  1:3815581/1‑93 (MQ=255)
                                                              aaTCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAg                          >  1:4282365/1‑93 (MQ=255)
                                                               aTCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGt                         <  1:102477/93‑1 (MQ=255)
                                                               aTCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGt                         <  1:3281056/93‑1 (MQ=255)
                                                                tCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAGCACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTg                        <  1:621767/93‑1 (MQ=255)
                                                                tCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTg                        >  2:4397193/1‑93 (MQ=255)
                                                                tCAGAGATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTg                        >  1:3385030/1‑93 (MQ=255)
                                                                  agagATTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACCAATTTTTGCTCTGGCATAGTGaa                      >  1:2618682/1‑93 (MQ=255)
                                                                   gagaTTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAAt                     >  2:204356/1‑93 (MQ=255)
                                                                   gagaTTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAAt                     >  1:1397374/1‑93 (MQ=255)
                                                                   gagaTTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAAt                     <  2:4274738/93‑1 (MQ=255)
                                                                   gagaTTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAAt                     >  2:3734450/1‑93 (MQ=255)
                                                                      aTTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCa                  <  1:1022947/93‑1 (MQ=255)
                                                                      aTTTGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCa                  >  1:4241125/1‑93 (MQ=255)
                                                                      aTTTGCTCCAGCGAACAGGGATAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCa                  <  2:2783802/93‑1 (MQ=255)
                                                                       tttGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCAc                 >  2:4565/1‑93 (MQ=255)
                                                                         tGCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCAccc               <  2:1162797/93‑1 (MQ=255)
                                                                          gCTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCa              <  1:161802/93‑1 (MQ=255)
                                                                           cTCCAGCGAACAGGGGTAGCTCTGGGTAACCCCAATCAGCTGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAg             <  1:96492/93‑1 (MQ=255)
                                                                           cTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAg             <  1:3149519/93‑1 (MQ=255)
                                                                           cTCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAg             >  1:2418637/1‑93 (MQ=255)
                                                                            tCCAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAgg            <  1:871079/93‑1 (MQ=255)
                                                                              cAGCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCg          >  1:2175852/1‑93 (MQ=255)
                                                                                gCGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTa        >  1:1807709/1‑93 (MQ=255)
                                                                                 cGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAg       <  1:4445147/93‑1 (MQ=255)
                                                                                 cGAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAg       <  2:1106979/93‑1 (MQ=255)
                                                                                  gAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAgc      <  2:2665908/93‑1 (MQ=255)
                                                                                  gAACAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAgc      >  1:2449507/1‑93 (MQ=255)
                                                                                    aCAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAgcgg    >  2:2601449/1‑93 (MQ=255)
                                                                                     cAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAgcggc   >  2:242478/1‑93 (MQ=255)
                                                                                     cAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAgcggc   <  1:9122/93‑1 (MQ=255)
                                                                                     cAGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAgcggc   <  1:4513778/93‑1 (MQ=255)
                                                                                      aGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAgcggcg  >  1:1906157/1‑93 (MQ=255)
                                                                                      aGGGGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAgcggcg  >  2:1556044/1‑93 (MQ=255)
                                                                                         |                                                                                         
CGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCG  >  NC_000913/4632355‑4632533

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: