Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,830,970 T→C R173R (CGT→CGC xthA → exonuclease III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,830,9700TC99.6% 847.6 / NA 236R173R (CGT→CGCxthAexonuclease III
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (115/120);  total (116/120)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 243 total aligned reads displayed.

TAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGTACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAA  >  NC_000913/1830880‑1831059
                                                                                          |                                                                                         
tAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCAc                                                                                         >  2:2432319/1‑93 (MQ=255)
tAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCAc                                                                                         <  2:4129918/93‑1 (MQ=255)
tAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCAc                                                                                         <  2:3215522/93‑1 (MQ=255)
tAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCAc                                                                                         <  2:3594965/93‑1 (MQ=255)
 aaaccGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCAcc                                                                                        <  1:5272920/90‑1 (MQ=255)
  aTCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCg                                                                                       >  2:1959386/1‑93 (MQ=255)
    ccGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGt                                                                                     <  2:1178423/93‑1 (MQ=255)
    ccGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGt                                                                                     <  2:5555795/93‑1 (MQ=255)
    ccGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGt                                                                                     <  1:2234803/93‑1 (MQ=255)
      ggTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTaa                                                                                   >  2:781466/1‑93 (MQ=255)
       gTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTaaa                                                                                  >  1:1958500/1‑93 (MQ=255)
       gTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTaaa                                                                                  <  2:1789760/93‑1 (MQ=255)
       gTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTaaa                                                                                  <  2:6232799/93‑1 (MQ=255)
         aCTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATg                                                                                >  2:4561775/1‑93 (MQ=255)
            gATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGctc                                                                             <  1:2463876/93‑1 (MQ=255)
             aTTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGctct                                                                            >  2:3779833/1‑93 (MQ=255)
             aTTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGctct                                                                            >  1:5661746/1‑93 (MQ=255)
              ttATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCtt                                                                           <  1:1521196/93‑1 (MQ=255)
               tATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCttt                                                                          >  2:3834608/1‑93 (MQ=255)
                 tGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTcc                                                                        <  2:485605/93‑1 (MQ=255)
                 tGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTcc                                                                        <  2:1515130/93‑1 (MQ=255)
                 tGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTcc                                                                        <  1:674514/93‑1 (MQ=255)
                  gggCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCt                                                                       >  1:3516288/1‑93 (MQ=255)
                   ggCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTg                                                                      <  1:4178728/93‑1 (MQ=255)
                    gCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGc                                                                     <  1:376172/93‑1 (MQ=255)
                    gCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGc                                                                     >  1:3093079/1‑93 (MQ=255)
                        tatGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGa                                                                 <  2:2155353/93‑1 (MQ=255)
                         atGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGaa                                                                <  1:5545276/93‑1 (MQ=255)
                         atGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGaa                                                                <  2:6037348/93‑1 (MQ=255)
                          tGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAag                                                               >  1:1591517/1‑93 (MQ=255)
                          tGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAag                                                               >  1:2704709/1‑93 (MQ=255)
                           gAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAaga                                                              <  2:4226406/93‑1 (MQ=255)
                            aaTATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAagag                                                             <  2:3114704/93‑1 (MQ=255)
                            aaTATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAagag                                                             <  1:5326076/93‑1 (MQ=255)
                             atatCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAgc                                                            >  1:1775303/1‑93 (MQ=255)
                              tatCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAgcg                                                           >  1:117073/1‑93 (MQ=255)
                              tatCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAgcg                                                           >  1:4524624/1‑93 (MQ=255)
                              tatCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAgcg                                                           >  1:5843610/1‑93 (MQ=255)
                              tatCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAgcg                                                           <  1:6339909/93‑1 (MQ=255)
                              tatCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAgcg                                                           >  1:6246211/1‑93 (MQ=255)
                               atCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAgcgc                                                          <  1:4465563/93‑1 (MQ=255)
                               atCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAgcgc                                                          <  1:5780204/93‑1 (MQ=255)
                                tCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGcgcg                                                         <  1:4208063/93‑1 (MQ=255)
                                  aGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGaa                                                       <  1:4906448/93‑1 (MQ=255)
                                    cccTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAatg                                                     <  2:535333/93‑1 (MQ=255)
                                     ccTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAatgg                                                    >  2:6099787/1‑93 (MQ=255)
                                        tcAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGGGCGAatggatg                                                 <  1:2947584/92‑1 (MQ=255)
                                        aCAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAatggatg                                                 <  2:751131/93‑1 (MQ=255)
                                          aGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAAtggatgga                                               <  1:2665921/93‑1 (MQ=255)
                                           gATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGAc                                              >  1:2901425/1‑93 (MQ=255)
                                           gATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGAc                                              <  2:4995384/93‑1 (MQ=255)
                                            aTCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACa                                             <  1:4350707/93‑1 (MQ=255)
                                             tCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAg                                            <  1:1401873/93‑1 (MQ=255)
                                             tCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAg                                            >  1:3633735/1‑93 (MQ=255)
                                              cTGGATATCGGCATTGGCGAAGAAACCCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAgg                                           <  2:2923323/93‑1 (MQ=255)
                                               tGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGc                                          >  2:5166405/1‑93 (MQ=255)
                                               tGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGc                                          >  1:2793632/1‑93 (MQ=255)
                                               tGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGc                                          <  1:165840/93‑1 (MQ=255)
                                                ggATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCt                                         <  1:6130408/93‑1 (MQ=255)
                                                 gATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCtg                                        <  2:4438529/93‑1 (MQ=255)
                                                 gATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCtg                                        >  2:6597020/1‑93 (MQ=255)
                                                   tatCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCtgat                                      >  1:5133902/1‑93 (MQ=255)
                                                    atCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCtgatg                                     <  2:3856587/93‑1 (MQ=255)
                                                      tggCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAg                                   <  2:5252541/92‑1 (MQ=255)
                                                      cGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAg                                   >  1:1452018/1‑93 (MQ=255)
                                                       ggCATTGTCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTACCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGc                                  <  2:5868867/93‑1 (MQ=255)
                                                       ggCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGc                                  >  2:5603376/1‑93 (MQ=255)
                                                        gCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCt                                 <  2:2037046/93‑1 (MQ=255)
                                                        gCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCt                                 >  2:737144/1‑93 (MQ=255)
                                                         cATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTg                                >  1:96727/1‑93 (MQ=255)
                                                          aTTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTgg                               >  2:1831049/1‑93 (MQ=255)
                                                           ttGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTggg                              <  1:1902459/93‑1 (MQ=255)
                                                           ttGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTggg                              >  1:46142/1‑93 (MQ=255)
                                                            tGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTgggg                             <  1:2736223/93‑1 (MQ=255)
                                                              gCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGt                           <  2:2726744/93‑1 (MQ=255)
                                                                gaagaaAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTg                         >  1:1753854/1‑93 (MQ=255)
                                                                gaagaaAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTg                         <  1:1326209/93‑1 (MQ=255)
                                                                  agaaAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGt                       >  1:5792670/1‑93 (MQ=255)
                                                                   gaaAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTc                      <  2:3370689/93‑1 (MQ=255)
                                                                   gaaAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTc                      >  1:4886430/1‑93 (MQ=255)
                                                                   gaaAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTc                      >  2:6523964/1‑93 (MQ=255)
                                                                   gaaAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTc                      <  2:1912170/93‑1 (MQ=255)
                                                                   gaaAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTc                      <  2:2354960/93‑1 (MQ=255)
                                                                     aaaCCGTAAGCGCTGGCTGTGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGa                    >  1:6496116/1‑93 (MQ=255)
                                                                     aaaCCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGa                    >  2:2493378/1‑93 (MQ=255)
                                                                      aaCCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGAt                   <  2:1302483/93‑1 (MQ=255)
                                                                      aaCCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGAt                   >  1:4943951/1‑93 (MQ=255)
                                                                        ccGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATAc                 >  2:6441915/1‑93 (MQ=255)
                                                                        ccGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATAc                 >  2:5465983/1‑93 (MQ=255)
                                                                           tAATCGCTTGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCtt              <  1:1212218/93‑1 (MQ=255)
                                                                           tAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCtt              <  2:5704016/93‑1 (MQ=255)
                                                                            aaGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTc             >  2:4405458/1‑93 (MQ=255)
                                                                            aaGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGGCGATACCTTc             >  2:1116965/1‑93 (MQ=255)
                                                                                gctggctgCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGcc         >  1:5975725/1‑93 (MQ=255)
                                                                                gctggctgCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACACGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGcc         >  2:4447234/1‑93 (MQ=255)
                                                                                 ctggctgCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCa        >  1:4216135/1‑93 (MQ=255)
                                                                                  tggctgCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCAt       <  1:5071897/93‑1 (MQ=255)
                                                                                  tggctgCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGCACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCAt       <  1:4584761/93‑1 (MQ=255)
                                                                                      tgCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGa   >  1:718019/1‑93 (MQ=255)
                                                                                       gcgcACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGaa  <  1:2096732/93‑1 (MQ=255)
                                                                                          |                                                                                         
TAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGTACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAA  >  NC_000913/1830880‑1831059

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: