Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,885,627 A→G L130P (CTG→CCG)  rnd ← RNase D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,885,6270AG83.2% 532.5 / 110.2 250L130P (CTG→CCG) rndRNase D
Reads supporting (aligned to +/- strand):  ref base A (23/19);  new base G (105/103);  total (128/122)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-01

Only 100 of 251 total aligned reads displayed.

AACAGATACCAGACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAG  >  NC_000913/1885537‑1885719
                                                                                          |                                                                                            
aaCAGATACCAGACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGc                                                                                            <  1:460252/93‑1 (MQ=255)
 aCAGATACCAGACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCg                                                                                           <  2:5251251/93‑1 (MQ=255)
 aCAGATACCAGACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCg                                                                                           >  1:2684031/1‑93 (MQ=255)
     atcccAGACTCCCGCCCCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAAc                                                                                       <  1:2761831/90‑1 (MQ=255)
      tACCAGACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACg                                                                                      >  1:773488/1‑93 (MQ=255)
      tACCAGACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACg                                                                                      <  2:4286412/93‑1 (MQ=255)
       cccAGACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGc                                                                                     <  1:4050906/92‑1 (MQ=255)
       aCCAGACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGc                                                                                     <  2:288803/93‑1 (MQ=255)
        ccAGACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGcc                                                                                    >  1:30653/1‑93 (MQ=255)
        ccAGACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGcc                                                                                    >  1:1044132/1‑93 (MQ=255)
        ccAGACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGcc                                                                                    <  2:5355497/93‑1 (MQ=255)
         cAGACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCg                                                                                   >  1:4942063/1‑93 (MQ=255)
         cAGACACCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCg                                                                                   <  1:3212201/93‑1 (MQ=255)
          aGACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCgg                                                                                  >  1:5494105/1‑93 (MQ=255)
           gACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGa                                                                                 >  2:1081274/1‑93 (MQ=255)
           gACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGa                                                                                 <  2:5757560/93‑1 (MQ=255)
           gACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGa                                                                                 >  2:2640602/1‑93 (MQ=255)
           gACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGa                                                                                 >  1:4905097/1‑93 (MQ=255)
           gACATCCGCCCCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGa                                                                                 <  2:5306615/93‑1 (MQ=255)
            aCATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGaa                                                                                >  2:34522/1‑93 (MQ=255)
             cATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAAt                                                                               >  1:5643242/1‑93 (MQ=255)
               tCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATAc                                                                             <  2:4468630/93‑1 (MQ=255)
               tCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATAc                                                                             >  1:2294388/1‑93 (MQ=255)
               tCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATAc                                                                             >  1:3438078/1‑93 (MQ=255)
               tCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATAc                                                                             >  2:7120009/1‑93 (MQ=255)
                 cgccgcTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATActc                                                                           <  2:3843065/93‑1 (MQ=255)
                  gccgcTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATActct                                                                          >  1:3978344/1‑93 (MQ=255)
                   ccgcTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCtt                                                                         <  2:6432779/93‑1 (MQ=255)
                    cgcTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTc                                                                        >  1:5202936/1‑93 (MQ=255)
                    cgcTGCGTATTCACACTGACGTTCGGTCAGCGGCCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTc                                                                        >  1:4231846/1‑93 (MQ=255)
                     gcTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTcc                                                                       >  2:2761232/1‑93 (MQ=255)
                       tGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTccac                                                                     >  2:1373988/1‑93 (MQ=255)
                       tGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTccac                                                                     >  1:5624793/1‑93 (MQ=255)
                       tGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTccac                                                                     >  1:7272369/1‑93 (MQ=255)
                        gCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTccacc                                                                    >  2:1730280/1‑93 (MQ=255)
                         cGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTccacca                                                                   <  2:3286775/93‑1 (MQ=255)
                          gTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCAt                                                                  >  2:4607825/1‑93 (MQ=255)
                          gTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCAt                                                                  >  1:5908150/1‑93 (MQ=255)
                             ttCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGa                                                               >  1:144992/1‑93 (MQ=255)
                             tccccacGGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGa                                                               <  1:62554/89‑1 (MQ=255)
                              tCACACTGACGTTCGGTCATCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGaa                                                              <  2:3957409/93‑1 (MQ=255)
                              tCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGaa                                                              <  1:7004619/93‑1 (MQ=255)
                              tCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGGAACGCCGGAATACTCTTCCACCATGGaa                                                              >  2:4138112/1‑93 (MQ=255)
                               cacaCTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAg                                                             <  2:120386/93‑1 (MQ=255)
                                acacTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGc                                                            <  1:3666990/93‑1 (MQ=255)
                                 cacTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCg                                                           >  2:1970972/1‑93 (MQ=255)
                                    tGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGaaa                                                        <  2:2313757/93‑1 (MQ=255)
                                    tGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGaaa                                                        >  1:2372855/1‑93 (MQ=255)
                                       cGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAAccc                                                     >  1:1948772/1‑93 (MQ=255)
                                        gTTGGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAAcccc                                                    <  1:3809754/93‑1 (MQ=255)
                                         ttCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCa                                                   >  1:95598/1‑93 (MQ=255)
                                          tCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCAt                                                  <  2:5543905/93‑1 (MQ=255)
                                           cGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATg                                                 <  1:6531370/93‑1 (MQ=255)
                                             gTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGAc                                               <  2:3636680/93‑1 (MQ=255)
                                              tCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACa                                              >  1:5127436/1‑93 (MQ=255)
                                                aGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATc                                            >  1:697623/1‑93 (MQ=255)
                                                aGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATc                                            <  1:6336226/93‑1 (MQ=255)
                                                aGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATc                                            <  2:7363292/93‑1 (MQ=255)
                                                 gCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCg                                           >  2:6292396/1‑93 (MQ=255)
                                                 gCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCg                                           <  2:6752803/93‑1 (MQ=255)
                                                  cGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCgg                                          <  1:2934610/93‑1 (MQ=255)
                                                   ggTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCggg                                         <  2:7293670/93‑1 (MQ=255)
                                                    gTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGc                                        >  2:1391602/1‑93 (MQ=255)
                                                    gTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGc                                        >  2:5039468/1‑93 (MQ=255)
                                                       tGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTc                                     <  2:5956322/93‑1 (MQ=255)
                                                         gccagccaGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCg                                   <  1:7071470/93‑1 (MQ=255)
                                                            agccagTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAg                                <  1:1613153/93‑1 (MQ=255)
                                                             gccagTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGa                               <  2:1279822/93‑1 (MQ=255)
                                                              ccagTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGaa                              <  2:3232757/93‑1 (MQ=255)
                                                              ccagTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGaa                              <  1:6870137/93‑1 (MQ=255)
                                                               cagTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAg                             <  2:6997075/93‑1 (MQ=255)
                                                               cagTCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAg                             <  2:6843622/93‑1 (MQ=255)
                                                                  tCGGTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGgca                          <  2:2141422/93‑1 (MQ=255)
                                                                  tCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGgca                          <  2:2244695/93‑1 (MQ=255)
                                                                    ggTGCGCGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGgcagc                        <  2:5622990/93‑1 (MQ=255)
                                                                    ggTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGgcagc                        <  1:6160283/93‑1 (MQ=255)
                                                                       gcgcGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAg                     <  2:3913675/93‑1 (MQ=255)
                                                                       gcgcGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAg                     <  1:3740961/93‑1 (MQ=255)
                                                                       gcgcGATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAg                     >  1:453851/1‑93 (MQ=255)
                                                                          cgATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGAt                  <  1:3186799/93‑1 (MQ=255)
                                                                           gATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGAtt                 >  2:7527137/1‑93 (MQ=255)
                                                                           gATTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGAtt                 >  2:2439145/1‑93 (MQ=255)
                                                                           gATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGAtt                 >  2:4703848/1‑93 (MQ=255)
                                                                            aTTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGAttt                <  1:1688929/93‑1 (MQ=255)
                                                                            aTTCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGAttt                <  1:5088612/93‑1 (MQ=255)
                                                                            aTTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGAttt                >  1:1189286/1‑93 (MQ=255)
                                                                             ttCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTg               >  1:1251739/1‑93 (MQ=255)
                                                                             ttCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTg               <  1:2096245/93‑1 (MQ=255)
                                                                             ttCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTg               <  2:6506579/93‑1 (MQ=255)
                                                                              tCACTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGc              <  1:5564289/93‑1 (MQ=255)
                                                                              tCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGc              >  1:5257487/1‑93 (MQ=255)
                                                                                aCTCTTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCgt            <  1:3140252/93‑1 (MQ=255)
                                                                                 ctctTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCgtg           >  2:6319933/1‑93 (MQ=255)
                                                                                 ctctTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCgtg           <  2:19167/93‑1 (MQ=255)
                                                                                 ctctTGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCgtg           <  2:397142/93‑1 (MQ=255)
                                                                                   ctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtc         >  1:4149216/1‑93 (MQ=255)
                                                                                     tGTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaa       <  2:4354153/93‑1 (MQ=255)
                                                                                      gTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaat      <  1:5729015/93‑1 (MQ=255)
                                                                                      gTCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaat      <  1:2631817/93‑1 (MQ=255)
                                                                                       tCCGGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaatc     <  1:7397787/93‑1 (MQ=255)
                                                                                          |                                                                                            
AACAGATACCAGACATCCGCCGCTGCGTATTCACACTGACGTTCGGTCAGCGGTCTGGCCAGCCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAG  >  NC_000913/1885537‑1885719

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