Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,141,015 T→C Q859R (CAG→CGG)  rne ← RNase E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,141,0150TC87.4% 673.2 / 77.0 270Q859R (CAG→CGG) rneRNase E
Reads supporting (aligned to +/- strand):  ref base T (14/20);  new base C (114/122);  total (128/142)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

Only 100 of 275 total aligned reads displayed.

ACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTAC  >  NC_000913/1140925‑1141103
                                                                                          |                                                                                        
acaacTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGcgccgc                                                                                        >  2:5911717/1‑93 (MQ=255)
  aacTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGctc                                                                                      <  1:1234048/93‑1 (MQ=255)
  aacTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGctc                                                                                      <  1:104426/93‑1 (MQ=255)
  aacTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGctc                                                                                      >  1:4838408/1‑93 (MQ=255)
   acTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGctct                                                                                     >  1:3892366/1‑93 (MQ=255)
     tGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTc                                                                                   <  1:5175576/93‑1 (MQ=255)
     tGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTc                                                                                   <  1:2402215/93‑1 (MQ=255)
         gcgcTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAAcc                                                                               >  1:6594579/1‑93 (MQ=255)
          cgcTAACAACCGGTTCGATAGCGGCGGTGACAGGGACCTCAGTAACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCt                                                                              <  1:371334/93‑1 (MQ=255)
          cgcTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCt                                                                              <  2:630986/93‑1 (MQ=255)
          cgcTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCt                                                                              <  1:3253720/93‑1 (MQ=255)
           gcTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTg                                                                             >  1:6309445/1‑93 (MQ=255)
            cTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGt                                                                            >  1:550371/1‑93 (MQ=255)
             tAACAGCCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTa                                                                           <  2:5546804/93‑1 (MQ=255)
              aacaacCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTAc                                                                          >  1:3786620/1‑93 (MQ=255)
               acaacCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACa                                                                         >  1:4634286/1‑93 (MQ=255)
               acaacCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACa                                                                         <  1:5985981/93‑1 (MQ=255)
                caacCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACAt                                                                        <  2:4082995/93‑1 (MQ=255)
                caacCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACAt                                                                        <  1:936345/93‑1 (MQ=255)
                  acCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCt                                                                      >  1:4701170/1‑93 (MQ=255)
                  acCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCt                                                                      >  1:2499680/1‑93 (MQ=255)
                   ccGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCtt                                                                     >  2:172611/1‑93 (MQ=255)
                   ccGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCtt                                                                     <  1:5698956/93‑1 (MQ=255)
                   ccGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCtt                                                                     >  1:2940983/1‑93 (MQ=255)
                   ccGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCtt                                                                     >  2:4901766/1‑93 (MQ=255)
                    cGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTg                                                                    >  2:3993864/1‑93 (MQ=255)
                     ggTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGc                                                                   >  2:5495943/1‑93 (MQ=255)
                     ggTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGc                                                                   >  2:6039138/1‑93 (MQ=255)
                     ggTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGc                                                                   >  2:668185/1‑93 (MQ=255)
                     ggTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGc                                                                   >  1:1839865/1‑93 (MQ=255)
                      gTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCg                                                                  <  2:6898857/93‑1 (MQ=255)
                      gTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCg                                                                  >  2:1000213/1‑93 (MQ=255)
                      gTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCg                                                                  <  1:6409949/93‑1 (MQ=255)
                      gTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCg                                                                  >  2:6190907/1‑93 (MQ=255)
                         cGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGAc                                                               <  2:3503335/93‑1 (MQ=255)
                         cGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGAc                                                               >  2:375883/1‑93 (MQ=255)
                           aTAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGt                                                             >  2:829840/1‑93 (MQ=255)
                            tAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTa                                                            <  1:398124/93‑1 (MQ=255)
                            tAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTa                                                            <  1:216937/93‑1 (MQ=255)
                            tAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTa                                                            >  2:5173304/1‑93 (MQ=255)
                            tAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCCGACGTa                                                            >  1:6125110/1‑93 (MQ=255)
                             aGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTAc                                                           <  1:4522657/93‑1 (MQ=255)
                                ggcggcgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAAt                                                        <  2:2426655/93‑1 (MQ=255)
                                 gcggcgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAAtt                                                       <  2:3208695/93‑1 (MQ=255)
                                   ggcgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTgg                                                     <  2:2666001/93‑1 (MQ=255)
                                    gcgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGa                                                    <  1:143105/93‑1 (MQ=255)
                                     cgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGAt                                                   <  1:3098040/93‑1 (MQ=255)
                                     cgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGAt                                                   <  1:6765052/93‑1 (MQ=255)
                                     cgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGAt                                                   <  2:6518241/93‑1 (MQ=255)
                                      gACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATa                                                  <  1:3397789/93‑1 (MQ=255)
                                       aCAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAg                                                 <  1:3029690/93‑1 (MQ=255)
                                       aCAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAg                                                 >  2:3642010/1‑93 (MQ=255)
                                       aCAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAg                                                 >  2:4359001/1‑93 (MQ=255)
                                       aCAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAg                                                 <  1:2094453/93‑1 (MQ=255)
                                       aCAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAg                                                 >  1:1316095/1‑93 (MQ=255)
                                          gggACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCgg                                              <  2:7208719/93‑1 (MQ=255)
                                          gggACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCgg                                              >  2:2091041/1‑93 (MQ=255)
                                           ggACCTCAGTCACCATCGGCTGCAGATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGa                                             <  2:4784907/93‑1 (MQ=255)
                                           ggACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGa                                             <  1:2586254/93‑1 (MQ=255)
                                            gACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGAt                                            <  1:1576946/93‑1 (MQ=255)
                                             aCCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATc                                           <  2:5091228/93‑1 (MQ=255)
                                              ccTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATcc                                          >  2:5672359/1‑93 (MQ=255)
                                               cTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCa                                         <  1:5325928/93‑1 (MQ=255)
                                                tCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAg                                        >  1:4069316/1‑93 (MQ=255)
                                                tCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAg                                        >  2:7279376/1‑93 (MQ=255)
                                                tCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAg                                        >  1:1921117/1‑93 (MQ=255)
                                                  aGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGAc                                      <  1:5264136/93‑1 (MQ=255)
                                                  aGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGAc                                      <  1:2078964/93‑1 (MQ=255)
                                                   gTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACt                                     <  1:501120/93‑1 (MQ=255)
                                                    tCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACtt                                    <  1:6198021/93‑1 (MQ=255)
                                                      aCCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTg                                  >  2:2350525/1‑93 (MQ=255)
                                                       ccATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGc                                 <  2:5787843/93‑1 (MQ=255)
                                                        cATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGcc                                >  2:6152232/1‑93 (MQ=255)
                                                          tCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCag                              <  1:5825193/93‑1 (MQ=255)
                                                          tCGGCTGCACATGGACTTTTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGGTCCAGACATTGCCag                              >  1:771424/1‑93 (MQ=255)
                                                           cGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCaga                             <  1:745224/93‑1 (MQ=255)
                                                            ggCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCagag                            <  1:4180142/93‑1 (MQ=255)
                                                            ggCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCagag                            <  2:1027356/93‑1 (MQ=255)
                                                            ggCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCagag                            >  1:2288519/1‑93 (MQ=255)
                                                            ggCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCagag                            >  2:4610943/1‑93 (MQ=255)
                                                              cTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGc                          <  1:7179795/93‑1 (MQ=255)
                                                                gCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCa                        <  2:1548449/93‑1 (MQ=255)
                                                                 cacaTGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAg                       >  2:5750083/1‑93 (MQ=255)
                                                                   caTGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGtt                     >  1:3581171/1‑93 (MQ=255)
                                                                   caTGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGtt                     >  1:6976948/1‑93 (MQ=255)
                                                                    aTGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTc                    <  1:4384882/93‑1 (MQ=255)
                                                                    aTGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTc                    >  1:5667776/1‑93 (MQ=255)
                                                                     tGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTcc                   <  1:495370/93‑1 (MQ=255)
                                                                     tGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTcc                   >  2:1397197/1‑93 (MQ=255)
                                                                     tGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTcc                   >  2:6376718/1‑93 (MQ=255)
                                                                      gTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCg                  >  2:1119287/1‑93 (MQ=255)
                                                                       tACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCgg                 >  1:2476604/1‑93 (MQ=255)
                                                                        aCTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGga                <  2:4063961/93‑1 (MQ=255)
                                                                        aCTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGga                <  1:2795645/93‑1 (MQ=255)
                                                                          ttcttcCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGgaga              <  2:1832227/93‑1 (MQ=255)
                                                                           tcttcCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGAc             >  2:7415141/1‑93 (MQ=255)
                                                                            cttcCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGAcg            <  2:1622007/93‑1 (MQ=255)
                                                                               ccTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACgcgc         <  2:7083343/93‑1 (MQ=255)
                                                                                cTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCa        >  2:1958420/1‑93 (MQ=255)
                                                                                    tCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCt    >  2:7367631/1‑93 (MQ=255)
                                                                                          |                                                                                        
ACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTAC  >  NC_000913/1140925‑1141103

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