Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,863,065 G→A A199V (GCC→GTC)  b1781 ← putative aldehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,863,0650GA99.2% 415.8 / ‑5.5 121A199V (GCC→GTC) b1781putative aldehyde reductase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (55/64);  total (57/64)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 129 total aligned reads displayed.

GCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGGCCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCGAGA  >  NC_000913/1862975‑1863154
                                                                                          |                                                                                         
gCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGAcc                                                                                         <  1:3634959/93‑1 (MQ=255)
 ccAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCt                                                                                        >  1:3871157/1‑93 (MQ=255)
  caacaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTg                                                                                       >  1:602350/1‑93 (MQ=255)
    acaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGGCCTggg                                                                                     >  1:4171495/1‑93 (MQ=255)
    acaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTggg                                                                                     >  1:1093339/1‑93 (MQ=255)
    acaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTggg                                                                                     >  1:2732965/1‑93 (MQ=255)
     caaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGc                                                                                    <  2:946280/93‑1 (MQ=255)
     caaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGc                                                                                    <  1:3123741/93‑1 (MQ=255)
     caaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGc                                                                                    <  1:2857776/93‑1 (MQ=255)
     caaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGc                                                                                    <  2:1362615/93‑1 (MQ=255)
      aaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCt                                                                                   >  2:4399561/1‑93 (MQ=255)
       aTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTa                                                                                  >  2:3508934/1‑93 (MQ=255)
        tACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTaa                                                                                 <  2:3710776/93‑1 (MQ=255)
         aCTTGTGCCGCGCTTATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAAc                                                                                <  2:912934/93‑1 (MQ=255)
         aCTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAAc                                                                                <  1:2322225/93‑1 (MQ=255)
           ttGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACgg                                                                              >  2:742536/1‑93 (MQ=255)
            tgtgCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGa                                                                             >  1:3235632/1‑93 (MQ=255)
               gCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTg                                                                          >  2:4119276/1‑93 (MQ=255)
                ccGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGt                                                                         >  2:3643949/1‑93 (MQ=255)
                ccGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGt                                                                         <  1:3449333/93‑1 (MQ=255)
                ccGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGt                                                                         >  2:4594798/1‑93 (MQ=255)
                   cgcTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAg                                                                      <  2:4427260/93‑1 (MQ=255)
                     cTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGcc                                                                    <  2:4185893/93‑1 (MQ=255)
                      tGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCa                                                                   >  1:53272/1‑93 (MQ=255)
                       gATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCAt                                                                  >  2:378606/1‑93 (MQ=255)
                       gATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCAt                                                                  <  1:126868/93‑1 (MQ=255)
                        atatTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATc                                                                 >  2:1205168/1‑93 (MQ=255)
                         tatTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCa                                                                <  1:1427916/93‑1 (MQ=255)
                         tatTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCa                                                                >  1:4411342/1‑93 (MQ=255)
                          atTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAc                                                               >  1:71497/1‑93 (MQ=255)
                          atTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAc                                                               >  1:510317/1‑93 (MQ=255)
                          atTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAc                                                               <  1:1841395/93‑1 (MQ=255)
                          atTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAc                                                               <  2:150229/93‑1 (MQ=255)
                          atTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAc                                                               <  1:1016322/93‑1 (MQ=255)
                           ttGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAcc                                                              >  2:3330646/1‑93 (MQ=255)
                           ttGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAcc                                                              >  1:1054583/1‑93 (MQ=255)
                            tgtgAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCg                                                             <  2:1098302/93‑1 (MQ=255)
                             gtgAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCgg                                                            <  1:1472975/93‑1 (MQ=255)
                                aGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCAt                                                         <  1:2043194/93‑1 (MQ=255)
                                   aTGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATctg                                                      >  2:3703363/1‑93 (MQ=255)
                                    tgtgCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATctgc                                                     <  2:3932248/93‑1 (MQ=255)
                                      tgCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATctgctg                                                   >  2:3976940/1‑93 (MQ=255)
                                       gCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGt                                                  <  2:3356191/93‑1 (MQ=255)
                                           tttCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCt                                              <  1:1908594/93‑1 (MQ=255)
                                           tttCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCt                                              <  2:484252/93‑1 (MQ=255)
                                            ttCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTg                                             <  2:3219798/93‑1 (MQ=255)
                                            ttCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTg                                             >  2:1984661/1‑93 (MQ=255)
                                             tCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTgg                                            <  1:2944119/93‑1 (MQ=255)
                                              cGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGc                                           >  2:1592263/1‑93 (MQ=255)
                                               gTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCa                                          >  2:470960/1‑93 (MQ=255)
                                               gTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCa                                          >  2:2317385/1‑93 (MQ=255)
                                                  gACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCa                                       <  1:920079/93‑1 (MQ=255)
                                                   aCTACCGCGTTTTTTAACAGTCCATTGCGCACCCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAg                                      <  2:610729/93‑1 (MQ=255)
                                                   aCTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAg                                      >  2:1670630/1‑93 (MQ=255)
                                                   aCTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAg                                      <  1:3929818/93‑1 (MQ=255)
                                                     tACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAggg                                    <  1:3805571/93‑1 (MQ=255)
                                                     tACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAggg                                    <  1:3206689/93‑1 (MQ=255)
                                                       ccGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGa                                  <  2:2156036/93‑1 (MQ=255)
                                                        cgcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAg                                 >  1:1739060/1‑93 (MQ=255)
                                                        cgcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAg                                 <  2:3669081/93‑1 (MQ=255)
                                                        cgcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAg                                 >  2:2180478/1‑93 (MQ=255)
                                                        cgcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAg                                 <  1:1140662/93‑1 (MQ=255)
                                                         gcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGt                                >  1:2004028/1‑93 (MQ=255)
                                                         gcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGt                                <  1:3074879/93‑1 (MQ=255)
                                                           gTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGCGTAg                              >  2:204252/1‑93 (MQ=255)
                                                            gtttttAACAGTCCATTGCGCAACCGCCCGACGTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGa                             <  1:1506034/92‑1 (MQ=255)
                                                             tttttAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGAt                            <  1:3160584/93‑1 (MQ=255)
                                                             tttttAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGAt                            <  2:71679/93‑1 (MQ=255)
                                                             tttttAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGCGAGGAGAt                            >  1:1121468/1‑93 (MQ=255)
                                                               tttAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGGGTAGATCg                          <  2:2496426/93‑1 (MQ=255)
                                                                ttAACAGTCCATTGCGCAACCGCCCGTCCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCGGGCACCAGGGGAGTAGATCGt                         <  2:2051019/93‑1 (MQ=255)
                                                                ttAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATATGCTGTTGCTGGCACCAGGGGAGTAGATCGt                         >  1:202179/1‑93 (MQ=255)
                                                                ttAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATAACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGCTCGt                         >  2:2022047/1‑93 (MQ=255)
                                                                 tAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTa                        >  1:38271/1‑93 (MQ=255)
                                                                 tAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTa                        >  1:1134614/1‑93 (MQ=255)
                                                                 tAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTa                        <  1:2550426/93‑1 (MQ=255)
                                                                 tAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTa                        <  2:3158799/93‑1 (MQ=255)
                                                                   aCAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACt                      >  1:945130/1‑93 (MQ=255)
                                                                   aCAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACt                      >  2:1436309/1‑93 (MQ=255)
                                                                       tCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAAt                  <  1:678369/93‑1 (MQ=255)
                                                                         cATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTc                >  2:1715843/1‑93 (MQ=255)
                                                                         cATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTc                >  2:4306542/1‑93 (MQ=255)
                                                                          aTTGCCCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTcc               <  1:2461304/93‑1 (MQ=255)
                                                                           ttGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTCCTCAATTCCt              >  1:324707/1‑93 (MQ=255)
                                                                           ttGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCt              <  2:2161699/93‑1 (MQ=255)
                                                                            tGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTc             <  1:3593763/93‑1 (MQ=255)
                                                                             gcgcAACCGCCGGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCg            <  2:4051661/93‑1 (MQ=255)
                                                                             gcgcAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCg            <  1:2829444/93‑1 (MQ=255)
                                                                               gcAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTg          >  2:4033830/1‑93 (MQ=255)
                                                                                cAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGa         >  2:4613743/1‑93 (MQ=255)
                                                                                cAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGa         >  1:3956039/1‑93 (MQ=255)
                                                                                  aCCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAAc       >  1:1500752/1‑93 (MQ=255)
                                                                                  aCCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAAc       <  2:2313168/93‑1 (MQ=255)
                                                                                   ccGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAAcc      <  1:2794315/93‑1 (MQ=255)
                                                                                     gCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCga    <  2:302419/93‑1 (MQ=255)
                                                                                     gCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCga    <  1:4596262/93‑1 (MQ=255)
                                                                                      cccGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCACTTCCTCCTGAACCgag   >  2:714079/1‑93 (MQ=255)
                                                                                      cccGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCgag   <  2:1805083/93‑1 (MQ=255)
                                                                                       ccGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCgaga  >  1:3114398/1‑93 (MQ=255)
                                                                                       ccGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCgaga  <  2:3154295/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGGCCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCGAGA  >  NC_000913/1862975‑1863154

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