Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 528,636 T→C intergenic (+282/+233) b0501 → / ← b0502 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913528,6360TC99.6% 851.1 / ‑5.7 264intergenic (+282/+233)b0501/b0502hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (127/136);  total (128/136)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 270 total aligned reads displayed.

TCAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAG  >  NC_000913/528546‑528725
                                                                                          |                                                                                         
tCAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACgcgt                                                                                         <  2:3307845/93‑1 (MQ=255)
 cAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGtt                                                                                        <  1:2642173/93‑1 (MQ=255)
  agagagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGttt                                                                                       <  1:3932887/93‑1 (MQ=255)
   gagagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTa                                                                                      >  2:4212770/1‑93 (MQ=255)
   gagagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTa                                                                                      <  1:2813923/93‑1 (MQ=255)
    agagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTaa                                                                                     <  2:1428430/93‑1 (MQ=255)
     gagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAt                                                                                    <  2:1828435/93‑1 (MQ=255)
     gagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAt                                                                                    <  1:5494164/93‑1 (MQ=255)
      agTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAtg                                                                                   <  1:897082/93‑1 (MQ=255)
      agTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAtg                                                                                   <  2:1697208/93‑1 (MQ=255)
       gTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAtgt                                                                                  <  1:4387379/93‑1 (MQ=255)
        tttACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAtgtg                                                                                 <  2:2102134/93‑1 (MQ=255)
           aCATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGt                                                                              >  2:4151828/1‑93 (MQ=255)
            cATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTc                                                                             <  2:4322613/93‑1 (MQ=255)
              tAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTccc                                                                           <  2:5313945/93‑1 (MQ=255)
              tAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTccc                                                                           >  1:1383958/1‑93 (MQ=255)
                gAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTa                                                                         >  1:5607327/1‑93 (MQ=255)
                gAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTa                                                                         <  1:1515714/93‑1 (MQ=255)
                 aaTTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTAt                                                                        <  1:783801/93‑1 (MQ=255)
                 aaTTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTAt                                                                        <  2:3499053/93‑1 (MQ=255)
                   ttGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCt                                                                      >  2:1561120/1‑93 (MQ=255)
                   ttGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCt                                                                      >  1:1016597/1‑93 (MQ=255)
                    tGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTg                                                                     >  1:23494/1‑93 (MQ=255)
                    tGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTg                                                                     <  1:4935371/93‑1 (MQ=255)
                    tGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTg                                                                     >  2:5517770/1‑93 (MQ=255)
                    tGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTg                                                                     <  2:2741333/93‑1 (MQ=255)
                     gCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGa                                                                    >  2:1061452/1‑93 (MQ=255)
                     gCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGa                                                                    <  2:3757398/93‑1 (MQ=255)
                       aCTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGaaa                                                                  >  1:5718481/1‑93 (MQ=255)
                       aCTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGaaa                                                                  >  2:2543552/1‑93 (MQ=255)
                         tGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAaaat                                                                >  1:4466988/1‑93 (MQ=255)
                         tGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAaaat                                                                >  2:3172066/1‑93 (MQ=255)
                         tGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAaaat                                                                >  1:1551643/1‑93 (MQ=255)
                         tGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAaaat                                                                >  2:3980057/1‑93 (MQ=255)
                          ggTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAaaata                                                               >  2:3907131/1‑93 (MQ=255)
                               tttACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAataaataa                                                          >  1:6229190/1‑93 (MQ=255)
                                ttACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAaataat                                                         <  2:3109716/93‑1 (MQ=255)
                                ttACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAaataat                                                         <  2:1458361/93‑1 (MQ=255)
                                 tACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATc                                                        >  1:4426814/1‑93 (MQ=255)
                                  aCGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATcc                                                       <  2:3886181/93‑1 (MQ=255)
                                   cGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCt                                                      >  2:1511287/1‑93 (MQ=255)
                                   cGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCt                                                      <  1:1731452/93‑1 (MQ=255)
                                   cGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCt                                                      <  1:5375373/93‑1 (MQ=255)
                                   cGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCt                                                      <  1:3405873/93‑1 (MQ=255)
                                   cGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTGAATGTGAGTCCCTATCTGAAAATAAATAATCCt                                                      <  2:5863207/93‑1 (MQ=255)
                                    gATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCtt                                                     <  2:1012870/93‑1 (MQ=255)
                                     atatCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTc                                                    <  1:3837/93‑1 (MQ=255)
                                      tatCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTcc                                                   <  2:1090945/93‑1 (MQ=255)
                                       atCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCg                                                  <  1:4165321/93‑1 (MQ=255)
                                        tCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCgg                                                 >  2:4857035/1‑93 (MQ=255)
                                        tCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCgg                                                 >  2:899551/1‑93 (MQ=255)
                                        tCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCgg                                                 >  2:3544826/1‑93 (MQ=255)
                                          tGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGAt                                               >  1:194638/1‑93 (MQ=255)
                                          tGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGAt                                               >  2:4292904/1‑93 (MQ=255)
                                            aCATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTa                                             >  1:50676/1‑93 (MQ=255)
                                            aCATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTa                                             >  2:2916849/1‑93 (MQ=255)
                                             cATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTaa                                            >  2:6189384/1‑93 (MQ=255)
                                             cATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTaa                                            >  2:3070138/1‑93 (MQ=255)
                                               ttGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTaaaa                                          >  2:6082418/1‑93 (MQ=255)
                                                tgtgTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAaaat                                         <  2:1303852/93‑1 (MQ=255)
                                                 gtgtAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAaaata                                        >  1:2098313/1‑93 (MQ=255)
                                                 gtgtAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAaaata                                        >  2:5658006/1‑93 (MQ=255)
                                                  tgtAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAaaataa                                       >  1:1241867/1‑93 (MQ=255)
                                                   gtAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAaaataaa                                      <  1:4884655/93‑1 (MQ=255)
                                                        aCATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCtt                                 <  1:3883090/93‑1 (MQ=255)
                                                            atTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCg                             >  1:329400/1‑93 (MQ=255)
                                                             ttCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCgg                            >  2:508085/1‑93 (MQ=255)
                                                             ttCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCgg                            <  1:3134929/93‑1 (MQ=255)
                                                                 gcaaTGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaa                        <  1:2269888/90‑1 (MQ=255)
                                                                 gCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCgggcag                        >  2:4991937/1‑90 (MQ=255)
                                                                 gCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaa                        <  1:662971/93‑1 (MQ=255)
                                                                 gCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaa                        <  1:2028709/93‑1 (MQ=255)
                                                                 gCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaa                        <  2:1572409/93‑1 (MQ=255)
                                                                   cATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaaga                      >  2:364403/1‑93 (MQ=255)
                                                                    aTGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaagaa                     >  2:2130390/1‑93 (MQ=255)
                                                                    aTGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaagaa                     <  1:5812174/93‑1 (MQ=255)
                                                                    aTGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaagaa                     >  1:2566431/1‑93 (MQ=255)
                                                                    aTGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaagaa                     <  1:6037563/93‑1 (MQ=255)
                                                                       cATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAaagaaaga                  <  1:3001410/93‑1 (MQ=255)
                                                                       cATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAaagaaaga                  <  1:6192500/93‑1 (MQ=255)
                                                                         ttaattaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTTTTCCCGGGAAAGAAAGAgg                >  2:5406464/1‑93 (MQ=255)
                                                                         ttaattaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAgg                >  1:5243504/1‑93 (MQ=255)
                                                                         ttaattaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAgg                <  1:4114018/93‑1 (MQ=255)
                                                                          taattaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGa               <  2:2418611/93‑1 (MQ=255)
                                                                             ttaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAAt            <  2:6245433/93‑1 (MQ=255)
                                                                             ttaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAAt            >  2:795381/1‑93 (MQ=255)
                                                                             ttaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAAt            >  1:4091242/1‑93 (MQ=255)
                                                                              taaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATa           >  1:5970604/1‑93 (MQ=255)
                                                                               aaaaTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATaa          >  1:3335408/1‑93 (MQ=255)
                                                                               aaaaTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATaa          <  1:3730741/93‑1 (MQ=255)
                                                                                aaaTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATaaa         <  1:4287539/93‑1 (MQ=255)
                                                                                 aaTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAAc        <  1:5899291/93‑1 (MQ=255)
                                                                                 aaTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAAc        <  2:2793381/93‑1 (MQ=255)
                                                                                  aTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAAcc       <  2:5999848/93‑1 (MQ=255)
                                                                                    gTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAt     <  1:3503146/93‑1 (MQ=255)
                                                                                    gTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAt     >  1:4766201/1‑93 (MQ=255)
                                                                                     ttACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAtt    >  1:4015499/1‑93 (MQ=255)
                                                                                     ttACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAtt    <  1:84094/93‑1 (MQ=255)
                                                                                     ttACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAtt    >  2:2477110/1‑93 (MQ=255)
                                                                                     ttACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAtt    >  1:1982913/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAG  >  NC_000913/528546‑528725

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