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.6% 847.4 / NA 232A199V (GCC→GTC) b1781putative aldehyde reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (111/120);  minor base C (1/0);  total (112/120)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

GCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGGCCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCGAGA  >  NC_000913/1862975‑1863154
                                                                                          |                                                                                         
gCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGAcc                                                                                         <  2:3461008/93‑1 (MQ=255)
gCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTCGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGAcc                                                                                         >  1:6622931/1‑93 (MQ=255)
cccAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGAcc                                                                                         <  2:1482729/92‑1 (MQ=255)
  caacaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTg                                                                                       >  1:6420068/1‑93 (MQ=255)
   aacaaTACTTGTGCCGCGCTGATATTGTGAGCATGTTCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTgg                                                                                      <  2:4058598/93‑1 (MQ=255)
   aacaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTgg                                                                                      >  1:4544638/1‑93 (MQ=255)
   aacaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTgg                                                                                      <  1:4984240/93‑1 (MQ=255)
   aacaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTgg                                                                                      >  1:1457437/1‑93 (MQ=255)
   aacaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTgg                                                                                      <  1:5060169/93‑1 (MQ=255)
   aacaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTgg                                                                                      >  2:2632344/1‑93 (MQ=255)
   aacaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTgg                                                                                      >  1:2148416/1‑93 (MQ=255)
   aacaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTgg                                                                                      >  2:1258397/1‑93 (MQ=255)
   aacaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTACATTGCGCAACCGCCCGACCTgg                                                                                      >  1:509089/1‑93 (MQ=255)
     caaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGc                                                                                    <  1:3844551/93‑1 (MQ=255)
     caaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGc                                                                                    >  1:537600/1‑93 (MQ=255)
     caaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGc                                                                                    <  1:5859169/93‑1 (MQ=255)
     caaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGc                                                                                    <  2:2982593/93‑1 (MQ=255)
       aTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTa                                                                                  >  2:4798680/1‑93 (MQ=255)
        tACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTaa                                                                                 >  2:2927156/1‑93 (MQ=255)
        tACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTaa                                                                                 <  1:1642566/93‑1 (MQ=255)
        tACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTaa                                                                                 >  1:4439180/1‑93 (MQ=255)
           ttGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACgg                                                                              >  1:5780789/1‑93 (MQ=255)
           ttGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACgg                                                                              >  2:6143031/1‑93 (MQ=255)
                  gcgcTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTaa                                                                       <  1:3377606/93‑1 (MQ=255)
                    gcTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGc                                                                     >  2:1969989/1‑93 (MQ=255)
                    gcTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGc                                                                     >  1:1057127/1‑93 (MQ=255)
                     cTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGcc                                                                    <  2:364221/93‑1 (MQ=255)
                      tGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCa                                                                   >  2:3152793/1‑93 (MQ=255)
                      tGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCa                                                                   <  1:5215945/93‑1 (MQ=255)
                       gATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCAt                                                                  <  2:374728/93‑1 (MQ=255)
                         tatTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCa                                                                >  2:6614645/1‑93 (MQ=255)
                          atTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAc                                                               <  2:5911838/93‑1 (MQ=255)
                           ttGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAcc                                                              >  1:4098780/1‑93 (MQ=255)
                           ttGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAcc                                                              >  1:2350204/1‑93 (MQ=255)
                             gtgAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCgg                                                            <  1:348109/93‑1 (MQ=255)
                              tgAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGc                                                           >  1:1849392/1‑93 (MQ=255)
                                aGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCAt                                                         >  2:4909874/1‑93 (MQ=255)
                                aGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCAt                                                         >  1:1105779/1‑93 (MQ=255)
                                 gCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATc                                                        >  1:2389303/1‑93 (MQ=255)
                                 gCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATc                                                        >  2:3397949/1‑93 (MQ=255)
                                  cATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATct                                                       <  1:208525/93‑1 (MQ=255)
                                   aTGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATctg                                                      <  1:1690622/93‑1 (MQ=255)
                                    tgtgCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATctgc                                                     >  1:3680488/1‑93 (MQ=255)
                                     gtgCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATctgct                                                    <  2:744665/93‑1 (MQ=255)
                                       gCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGt                                                  >  2:6524875/1‑93 (MQ=255)
                                        gaaTTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACAGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGtt                                                 <  1:5002152/92‑1 (MQ=255)
                                        cAATTTCGTTGACTACCGCGTTTTTTAACCGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGGAAGCCATCACCGGCATCTGCTGtt                                                 >  2:2758715/1‑93 (MQ=255)
                                          aTTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGc                                               >  1:2584457/1‑93 (MQ=255)
                                           tttCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCt                                              >  2:4484926/1‑93 (MQ=255)
                                           tttCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCt                                              <  1:3140755/93‑1 (MQ=255)
                                            ttCGTGGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTg                                             <  1:1612422/93‑1 (MQ=255)
                                                ttGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCAc                                         >  2:3071450/1‑93 (MQ=255)
                                                 tGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCAcc                                        <  1:5694687/93‑1 (MQ=255)
                                                  gACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCa                                       >  2:5420434/1‑93 (MQ=255)
                                                   aCTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAg                                      >  2:2480847/1‑93 (MQ=255)
                                                    cTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAgg                                     <  2:6604283/93‑1 (MQ=255)
                                                     tACCGCGTTTTTTAACAGTCCATTGCGCACCCGCCCGGCCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAggg                                    <  1:5537596/93‑1 (MQ=255)
                                                     tACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAggg                                    <  1:5307216/93‑1 (MQ=255)
                                                     tACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAggg                                    <  1:475572/93‑1 (MQ=255)
                                                      aCCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAgggg                                   >  1:2032726/1‑93 (MQ=255)
                                                      aCCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAgggg                                   >  2:561224/1‑93 (MQ=255)
                                                      aCCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAgggg                                   <  1:4716237/93‑1 (MQ=255)
                                                      aCCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGACATCACCGGCATCTTCTGTTGCTGGCACCAgggg                                   >  2:5130706/1‑93 (MQ=255)
                                                        cgcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAg                                 <  1:5396755/93‑1 (MQ=255)
                                                        cgcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCAGTTGCTGGCACCAGGGGAg                                 <  1:6554828/93‑1 (MQ=255)
                                                         gcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGt                                <  1:3863143/93‑1 (MQ=255)
                                                         gcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGt                                >  1:2765227/1‑93 (MQ=255)
                                                            ttttttAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGa                             >  2:1769243/1‑93 (MQ=255)
                                                             tttttAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGAt                            <  1:362470/93‑1 (MQ=255)
                                                                 tAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTa                        >  1:1887628/1‑93 (MQ=255)
                                                                 tAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTa                        >  1:6003244/1‑93 (MQ=255)
                                                                 tAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTa                        <  2:3492737/93‑1 (MQ=255)
                                                                  aaCAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTAc                       >  2:4392573/1‑93 (MQ=255)
                                                                   aCAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACt                      >  2:3573733/1‑93 (MQ=255)
                                                                    cAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTc                     >  2:1997029/1‑93 (MQ=255)
                                                                     aGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCa                    >  2:1414155/1‑93 (MQ=255)
                                                                     aGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCa                    >  2:1847898/1‑93 (MQ=255)
                                                                        ccATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAAtt                 <  2:50893/93‑1 (MQ=255)
                                                                          aTTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTcc               >  2:6244842/1‑93 (MQ=255)
                                                                          aTTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTcc               >  2:4124333/1‑93 (MQ=255)
                                                                          aTTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTcc               <  1:5925377/93‑1 (MQ=255)
                                                                          aTTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTcc               <  1:6120877/93‑1 (MQ=255)
                                                                           ttGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCt              <  2:5486172/93‑1 (MQ=255)
                                                                             gcgcAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCg            >  2:4288504/1‑93 (MQ=255)
                                                                              cgcAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGt           <  1:4903112/93‑1 (MQ=255)
                                                                              cgcAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGt           <  2:956157/93‑1 (MQ=255)
                                                                              cgcAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGt           <  2:3886968/93‑1 (MQ=255)
                                                                               gcAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTg          >  2:4706064/1‑93 (MQ=255)
                                                                                 aaCCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGaa        <  2:1080540/93‑1 (MQ=255)
                                                                                  aCCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAAc       >  2:58760/1‑93 (MQ=255)
                                                                                  aCCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAAc       >  1:1002263/1‑93 (MQ=255)
                                                                                  aCCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAAc       >  1:6158484/1‑93 (MQ=255)
                                                                                  aCCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAAc       <  1:859694/93‑1 (MQ=255)
                                                                                   ccGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAAcc      <  2:4460811/93‑1 (MQ=255)
                                                                                     gCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCga    >  1:309516/1‑93 (MQ=255)
                                                                                     gCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCga    >  2:2225329/1‑93 (MQ=255)
                                                                                      cccGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCgag   >  1:52668/1‑93 (MQ=255)
                                                                                      cccGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCgag   >  1:2605923/1‑93 (MQ=255)
                                                                                       ccGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCgaga  <  1:3346488/93‑1 (MQ=255)
                                                                                       ccGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCgaga  <  1:1340590/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGGCCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCGAGA  >  NC_000913/1862975‑1863154

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