Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,165 C→T V89V (GTG→GTA ppiC ← peptidyl‑prolyl cis‑trans isomerase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,957,1650CT99.7% 1302.4 / ‑5.8 372V89V (GTG→GTAppiCpeptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (198/173);  total (198/174)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

Only 100 of 385 total aligned reads displayed.

ATGCCCTCACCCCGA‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGG  >  NC_000913/3957077‑3957254
                                                                                            |                                                                                         
aTGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTCGCGAAGGCCTTGCTATTAGTTGCGGTACAGCa                                                                                          >  2:1372492/1‑93 (MQ=21)
     cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTa                                                                                     >  1:722109/1‑93 (MQ=25)
      tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTaa                                                                                    <  2:3416009/93‑1 (MQ=25)
      tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTaa                                                                                    <  2:2522403/93‑1 (MQ=25)
      tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTaa                                                                                    >  1:7648439/1‑93 (MQ=25)
       cACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGGGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAat                                                                                   >  1:1637316/1‑93 (MQ=25)
       cACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAat                                                                                   >  2:5248539/1‑93 (MQ=37)
             gaccctctCCCCCAGGGGGAGGGAGAAAACACTCAAAGCCTTCTCCTGGAGAAGGGCTTGCTATTAGTTGCGGGACAGTACCTTAATGAtttg                                                                             >  2:2251075/5‑93 (MQ=255)
             gaccctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGGCTTGCTATTAGTTGCGGTACAGTACCTTAATGAtgtg                                                                             >  2:7565311/5‑93 (MQ=38)
              accctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGa                                                                            <  2:4433766/90‑1 (MQ=38)
              accctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGa                                                                            >  1:7716766/4‑93 (MQ=38)
              accctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGa                                                                            <  1:3379505/90‑1 (MQ=38)
               ccctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGat                                                                           <  1:2279135/91‑1 (MQ=38)
                cctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGata                                                                          <  1:1109746/92‑1 (MQ=39)
                cctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGata                                                                          >  1:1708125/2‑93 (MQ=39)
                cctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGata                                                                          >  1:2094018/2‑93 (MQ=39)
                 ctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGatat                                                                         >  2:4808518/1‑93 (MQ=255)
                 ctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGatat                                                                         >  1:7339230/1‑93 (MQ=255)
                   ctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATcc                                                                       <  2:3925715/93‑1 (MQ=255)
                    tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCg                                                                      <  1:7290944/93‑1 (MQ=255)
                    tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCg                                                                      >  1:7474902/1‑93 (MQ=255)
                    tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCg                                                                      >  1:4190471/1‑93 (MQ=255)
                    tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGGTGCGGTACAGTACCTTAATGATGTGATATCCg                                                                      >  1:2900624/1‑93 (MQ=255)
                    ccccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCg                                                                      <  2:5868741/92‑1 (MQ=255)
                      ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGaa                                                                    <  1:675440/93‑1 (MQ=255)
                      ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGaa                                                                    <  2:5829797/93‑1 (MQ=255)
                       cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAAc                                                                   <  2:1866729/93‑1 (MQ=255)
                        acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACt                                                                  >  1:3592343/1‑93 (MQ=255)
                        acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACt                                                                  >  1:439506/1‑93 (MQ=255)
                          agggagagggagAAAACACTCAATGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTgg                                                                <  2:4274989/93‑1 (MQ=255)
                           gggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTggg                                                               <  1:3886679/93‑1 (MQ=255)
                           gggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTggg                                                               >  1:4215571/1‑93 (MQ=255)
                           gggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTggg                                                               >  1:5560521/1‑93 (MQ=255)
                           gggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTACTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTggg                                                               <  1:4654605/93‑1 (MQ=255)
                             gagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGgtg                                                             <  1:6502674/93‑1 (MQ=255)
                             gagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGgtg                                                             <  1:4976138/93‑1 (MQ=255)
                             gagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGgtg                                                             >  2:2598855/1‑93 (MQ=255)
                              agagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGgtgt                                                            >  1:3808719/1‑93 (MQ=255)
                              agagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGgtgt                                                            <  2:7219774/93‑1 (MQ=255)
                                 gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCa                                                         <  1:6335355/93‑1 (MQ=255)
                                     gaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg                                                     <  2:4500612/93‑1 (MQ=255)
                                       aagcacTCAAGGCGTCCTCCTGGAGACGGCATTGCTATGAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc                                                   <  1:5372885/90‑1 (MQ=255)
                                       aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc                                                   >  1:5500478/1‑93 (MQ=255)
                                       aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc                                                   >  2:6108260/1‑93 (MQ=255)
                                       aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc                                                   >  1:7512297/1‑93 (MQ=255)
                                       aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGCGGGGCAGCGGGc                                                   >  1:2471842/1‑93 (MQ=255)
                                         acacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCg                                                 <  2:5474092/93‑1 (MQ=255)
                                          cacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCgg                                                >  2:1683118/1‑93 (MQ=255)
                                           acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGt                                               <  2:2117424/93‑1 (MQ=255)
                                            cTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTc                                              >  2:958030/1‑93 (MQ=255)
                                                aaggCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGc                                           >  2:1176033/4‑93 (MQ=255)
                                                 aggCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTCGGTGTGCAGCGGGCCGGTCGGCt                                          >  2:4298993/3‑93 (MQ=255)
                                                  ggCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                                         >  2:4044649/2‑93 (MQ=255)
                                                  ggCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                                         >  1:5088601/2‑93 (MQ=255)
                                                   gCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc                                        <  1:5719506/93‑1 (MQ=255)
                                                   gCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc                                        >  2:7893679/1‑93 (MQ=255)
                                                    ccTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                       <  1:6157104/93‑1 (MQ=255)
                                                    ccTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                       <  1:4756836/93‑1 (MQ=255)
                                                     cTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAg                                      <  2:1939713/93‑1 (MQ=255)
                                                     cTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAg                                      >  1:3119567/1‑93 (MQ=255)
                                                     cTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAg                                      <  1:3251358/93‑1 (MQ=255)
                                                      ttCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGt                                     >  1:997001/1‑93 (MQ=255)
                                                       tctcCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTa                                    >  1:7035870/1‑93 (MQ=255)
                                                         tcCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAcc                                  >  1:7868876/1‑93 (MQ=255)
                                                             ggAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGaa                              >  1:2657489/1‑92 (MQ=255)
                                                               agaAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                            <  2:3148480/93‑1 (MQ=255)
                                                               agaAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                            >  1:3345208/1‑93 (MQ=255)
                                                                gaAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAag                           >  2:1517121/1‑93 (MQ=255)
                                                                gaAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAag                           >  1:1853281/1‑93 (MQ=255)
                                                                 aaGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAaga                          >  1:3137233/1‑93 (MQ=255)
                                                                 aaGGCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAaga                          >  1:3559149/1‑93 (MQ=255)
                                                                   ggCCTTGCTATTAGTTGTGGCACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                        >  1:6031266/1‑93 (MQ=255)
                                                                    gCCTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                       >  2:2870479/1‑93 (MQ=255)
                                                                     ccTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaa                      >  1:1002653/1‑93 (MQ=255)
                                                                      cTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaaa                     >  2:5120443/1‑93 (MQ=255)
                                                                      cTTGCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaaa                     >  1:5159689/1‑93 (MQ=255)
                                                                         gCTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCa                  >  1:928877/1‑93 (MQ=255)
                                                                          cTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc                 >  1:1477933/1‑93 (MQ=255)
                                                                          cTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc                 <  2:2745536/93‑1 (MQ=255)
                                                                          cTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc                 <  2:7757420/93‑1 (MQ=255)
                                                                          cTATTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc                 >  2:1461823/1‑93 (MQ=255)
                                                                            aTTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACtt               <  2:6986390/93‑1 (MQ=255)
                                                                            aTTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACtt               >  2:8869/1‑93 (MQ=255)
                                                                            aTTAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACtt               <  2:4617293/93‑1 (MQ=255)
                                                                             ttAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACttt              >  2:6386984/1‑93 (MQ=255)
                                                                             ttAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACttt              <  1:2035935/93‑1 (MQ=255)
                                                                              tAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTa             <  2:5667471/93‑1 (MQ=255)
                                                                              tAGTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTa             <  1:4057616/93‑1 (MQ=255)
                                                                                gTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc           >  1:693438/1‑93 (MQ=255)
                                                                                gTTGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc           <  2:2390838/93‑1 (MQ=255)
                                                                                 ttGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTCCCGGCCAAGAGAAAACCACTTTATca          >  1:282349/1‑92 (MQ=255)
                                                                                 ttGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCg          >  1:1338999/1‑93 (MQ=255)
                                                                                 ttGCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCg          >  2:4936082/1‑93 (MQ=255)
                                                                                   gCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGaa        >  2:6180874/1‑93 (MQ=255)
                                                                                   gCGGTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGaa        <  2:6561232/93‑1 (MQ=255)
                                                                                    cGGTACAGTACGTTACTGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAAc       <  2:3808877/93‑1 (MQ=255)
                                                                                     ggTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACg      >  2:7319981/1‑93 (MQ=255)
                                                                                     ggTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACg      <  2:3073772/93‑1 (MQ=255)
                                                                                      gTACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGc     >  1:5285600/1‑93 (MQ=255)
                                                                                       tACAGTACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGcc    >  1:3085028/1‑93 (MQ=255)
                                                                                            |                                                                                         
ATGCCCTCACCCCGA‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGG  >  NC_000913/3957077‑3957254

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