Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,186,089 T→C I225T (ATC→ACC)  yqiG → putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,186,0890TC99.5% 685.6 / ‑5.5 218I225T (ATC→ACC) yqiGputative membrane protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (106/111);  total (107/111)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 224 total aligned reads displayed.

TACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCATTTTTCTGGCATTTCGTTGA  >  NC_000913/3186000‑3186178
                                                                                         |                                                                                         
tACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAAcccc                                                                                        >  1:1733490/1‑93 (MQ=255)
 aCTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCa                                                                                       <  1:2146367/93‑1 (MQ=255)
 aCTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCa                                                                                       <  2:2702969/93‑1 (MQ=255)
 aCTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCa                                                                                       <  1:2868519/93‑1 (MQ=255)
  cTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAg                                                                                      >  1:324495/1‑93 (MQ=255)
  cTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAg                                                                                      <  2:1552061/93‑1 (MQ=255)
  cTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAg                                                                                      <  2:1099795/93‑1 (MQ=255)
  cTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAg                                                                                      <  2:535982/93‑1 (MQ=255)
  cTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAg                                                                                      <  1:690057/93‑1 (MQ=255)
  cTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGACGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAg                                                                                      >  1:2350171/1‑93 (MQ=255)
    gATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCa                                                                                    >  1:4556221/1‑93 (MQ=255)
      ttACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCaat                                                                                  >  1:3515224/1‑93 (MQ=255)
      ttACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCaat                                                                                  <  1:2767020/93‑1 (MQ=255)
      ttACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCaat                                                                                  >  1:1905202/1‑93 (MQ=255)
         ccAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCaatgaa                                                                               <  2:1275522/93‑1 (MQ=255)
         ccAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCaatgaa                                                                               >  1:1002407/1‑93 (MQ=255)
         ccAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCaatgaa                                                                               <  2:2401638/93‑1 (MQ=255)
         ccAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCaatgaa                                                                               <  1:2023259/93‑1 (MQ=255)
          cAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCaatgaat                                                                              <  2:2173898/93‑1 (MQ=255)
            gTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGc                                                                            <  2:4245677/93‑1 (MQ=255)
                gATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCaaaa                                                                        <  2:1947383/93‑1 (MQ=255)
                gATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCaaaa                                                                        >  2:2010838/1‑93 (MQ=255)
                gATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCaaaa                                                                        <  2:4297140/93‑1 (MQ=255)
                 aTCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAGCCCCAGCAATGAATGCAAaac                                                                       <  2:2869187/93‑1 (MQ=255)
                 aTCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAaac                                                                       >  2:412886/1‑93 (MQ=255)
                 aTCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAaac                                                                       <  2:2244859/93‑1 (MQ=255)
                    aGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAaactaa                                                                    >  2:2134653/1‑93 (MQ=255)
                     gTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAaactaac                                                                   <  2:1710530/93‑1 (MQ=255)
                       cATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAactaacta                                                                 >  1:2977087/1‑93 (MQ=255)
                       cATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAactaacta                                                                 >  1:856063/1‑93 (MQ=255)
                       cATTTAACAACGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAactaacta                                                                 <  2:37832/93‑1 (MQ=255)
                         tttAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACtata                                                               >  1:1837550/1‑93 (MQ=255)
                         tttAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACtata                                                               <  1:4375956/93‑1 (MQ=255)
                         tttAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACtata                                                               <  2:437590/93‑1 (MQ=255)
                         tttAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACtata                                                               <  1:2594704/93‑1 (MQ=255)
                         tttAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACtata                                                               >  1:3605898/1‑93 (MQ=255)
                          ttAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAg                                                              >  2:2521989/1‑93 (MQ=255)
                          ttAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAg                                                              <  1:3759777/93‑1 (MQ=255)
                          ttAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAg                                                              >  1:3078004/1‑93 (MQ=255)
                          ttAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAg                                                              <  2:4504931/93‑1 (MQ=255)
                             aCAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTc                                                           <  1:784817/93‑1 (MQ=255)
                             aCAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTc                                                           >  2:3123080/1‑93 (MQ=255)
                               aaTGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCaa                                                         <  2:2800938/93‑1 (MQ=255)
                                 tGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATa                                                       <  2:449991/93‑1 (MQ=255)
                                  ggCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATAc                                                      <  1:4024260/93‑1 (MQ=255)
                                  ggCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATAc                                                      <  1:4407839/93‑1 (MQ=255)
                                  ggCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATAc                                                      <  2:748548/93‑1 (MQ=255)
                                  ggCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATAc                                                      <  1:3461990/93‑1 (MQ=255)
                                    cAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGa                                                    <  2:858913/93‑1 (MQ=255)
                                     aaaaGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGAt                                                   >  1:560990/1‑93 (MQ=255)
                                     aaaaGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGAt                                                   <  1:3259599/93‑1 (MQ=255)
                                     aaaaGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGAt                                                   <  2:3503869/93‑1 (MQ=255)
                                     aaaaGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGAt                                                   <  1:4218569/93‑1 (MQ=255)
                                      aaaGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATa                                                  >  2:1179148/1‑93 (MQ=255)
                                      aaaGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATa                                                  <  2:1380545/93‑1 (MQ=255)
                                        aGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTAAATACGATACt                                                <  1:4069167/93‑1 (MQ=255)
                                         gCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTg                                               >  2:358182/1‑93 (MQ=255)
                                           cAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGaa                                             >  1:394535/1‑93 (MQ=255)
                                           cAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGaa                                             <  2:2224841/93‑1 (MQ=255)
                                            aGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAAt                                            >  1:1317118/1‑93 (MQ=255)
                                            aGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAAt                                            >  1:3671742/1‑93 (MQ=255)
                                            aGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAAt                                            <  1:340877/93‑1 (MQ=255)
                                             ggCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATc                                           >  2:4542805/1‑93 (MQ=255)
                                             ggCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATc                                           <  1:2359930/93‑1 (MQ=255)
                                                 aCCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCctct                                       <  1:364540/93‑1 (MQ=255)
                                                  ccAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTa                                      >  2:4467044/1‑93 (MQ=255)
                                                   cAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAt                                     <  1:1053370/93‑1 (MQ=255)
                                                    aaCCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAtt                                    <  1:1545020/93‑1 (MQ=255)
                                                     aCCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAttt                                   <  2:3485478/93‑1 (MQ=255)
                                                      ccTGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAtttt                                  <  2:3088841/93‑1 (MQ=255)
                                                        tGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCg                                >  1:124301/1‑93 (MQ=255)
                                                        tGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCg                                >  2:3264213/1‑93 (MQ=255)
                                                        tGGATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATACGTCAATACGATACTGAATCCTCTATTTTCg                                >  2:4209394/1‑93 (MQ=255)
                                                          gATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGAc                              <  2:2891856/93‑1 (MQ=255)
                                                          gATTTTCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGAc                              <  2:3194970/93‑1 (MQ=255)
                                                            ttttCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGActc                            >  2:1879324/1‑93 (MQ=255)
                                                              ttCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCtt                          >  1:3720414/1‑93 (MQ=255)
                                                               tCCGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCttt                         <  1:2656962/93‑1 (MQ=255)
                                                                ccGCGTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTAAAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTc                        <  2:714954/93‑1 (MQ=255)
                                                                    gTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCAtt                    >  1:3167876/1‑93 (MQ=255)
                                                                    gTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCAtt                    <  1:1744543/93‑1 (MQ=255)
                                                                    gTATTTATTTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCAtt                    >  1:2538493/1‑93 (MQ=255)
                                                                     tatttattTGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCAttt                   >  2:2825824/1‑93 (MQ=255)
                                                                      ttttatttGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCAtttt                  <  2:3368755/92‑1 (MQ=255)
                                                                      atttatttGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCAtttt                  >  1:1117777/1‑93 (MQ=255)
                                                                       tttatttGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCAttttt                 <  1:1797587/93‑1 (MQ=255)
                                                                       tttatttGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCAttttt                 >  2:3863835/1‑93 (MQ=255)
                                                                       tttatttGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCAttttt                 <  2:3445174/93‑1 (MQ=255)
                                                                        ttatttGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCATTTTTc                >  1:1191648/1‑93 (MQ=255)
                                                                        ttatttGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCATTTTTc                >  1:2894985/1‑93 (MQ=255)
                                                                         tatttGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCATTTTTCt               <  1:1700358/93‑1 (MQ=255)
                                                                           tttGTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCATTTTTCTgg             >  1:1590025/1‑93 (MQ=255)
                                                                              gTTTCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCATTTTTCTGGCAt          >  1:3435032/1‑93 (MQ=255)
                                                                               tttCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCATTTTTCTGGCAtt         >  1:750643/1‑93 (MQ=255)
                                                                                 tCGCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCATTTTTCTGGCATTTc       <  1:1749048/93‑1 (MQ=255)
                                                                                   gCCCAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCATTTTTCTGGCATTTCGt     >  1:990468/1‑93 (MQ=255)
                                                                                    cccAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCATTTTTCTGGCATTTCGtt    <  2:4051070/93‑1 (MQ=255)
                                                                                     ccAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCATTTTTCTGGCATTTCGTTg   <  2:2305335/93‑1 (MQ=255)
                                                                                      cAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCATTTTTCTGGCATTTCGTTGa  <  2:1464613/93‑1 (MQ=255)
                                                                                      cAACCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCATTTTTCTGGCATTTCGTTGa  <  1:4264849/93‑1 (MQ=255)
                                                                                         |                                                                                         
TACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCATTTTTCTGGCATTTCGTTGA  >  NC_000913/3186000‑3186178

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: