Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 860,412 C→T G473E (GGG→GAG)  ybiW ← putative formate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913860,4120CT99.5% 706.7 / ‑5.7 203G473E (GGG→GAG) ybiWputative formate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (99/103);  total (99/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 227 total aligned reads displayed.

CGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGACGGCGGTT  >  NC_000913/860323‑860501
                                                                                         |                                                                                         
 gTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGcc                                                                                       <  2:3017344/93‑1 (MQ=255)
     aTCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGGATCATGCTCGCCTTcc                                                                                   >  1:1235126/1‑93 (MQ=255)
     aTCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTcc                                                                                   >  1:448317/1‑93 (MQ=255)
      tCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAAAACTTTGCCGCTGGTGGCATCATGCTCGCCTTccc                                                                                  >  2:1762081/1‑92 (MQ=255)
       cACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAg                                                                                 <  2:1707100/93‑1 (MQ=255)
        ccTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAgc                                                                                <  2:1721914/92‑1 (MQ=255)
        aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTTCAgc                                                                                >  2:2425611/1‑93 (MQ=255)
        aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGATCGGGGCATCATGCTAGCCTTCCCgc                                                                                >  2:3613591/1‑93 (MQ=255)
        aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACCTTGCCGCTGGTGGCATCCTGCTCGCCTTCCAgc                                                                                >  2:2581474/1‑93 (MQ=255)
          ttCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTTCCGCTGGTGGCATCATGCTCGCCTTCCAgcgc                                                                              >  1:969898/1‑93 (MQ=255)
          ttCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAgcgc                                                                              <  2:6253609/93‑1 (MQ=255)
          ttCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGCGGCATCCTGCTCGCCTTCCAgcgc                                                                              >  1:1159771/1‑93 (MQ=255)
             aTCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGcgccgc                                                                           <  2:2734355/93‑1 (MQ=255)
              tCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCGTCCAGCgccgcc                                                                          >  1:612181/1‑93 (MQ=255)
                  aGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCa                                                                      <  1:849418/93‑1 (MQ=255)
                  aGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCa                                                                      <  1:5282503/93‑1 (MQ=255)
                  aGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCa                                                                      <  1:3538864/93‑1 (MQ=255)
                    ttgttgAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCCGCGGCGCCAGCATc                                                                    >  1:2893931/1‑93 (MQ=255)
                    ttgttgAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATc                                                                    <  2:6569218/93‑1 (MQ=255)
                    ttgttgAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATc                                                                    >  1:7202491/1‑93 (MQ=255)
                     tgttgAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCa                                                                   <  2:4457312/93‑1 (MQ=255)
                      gttgAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCAc                                                                  <  2:3481958/93‑1 (MQ=255)
                      gttgAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCAc                                                                  <  2:645182/93‑1 (MQ=255)
                      gttgAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCAc                                                                  >  2:1354674/1‑93 (MQ=255)
                       ttgAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCCTCAcg                                                                 >  2:4421879/1‑93 (MQ=255)
                        tgAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCAcgc                                                                <  1:2239026/93‑1 (MQ=255)
                        tgAAGTAACCAGCCGACAACGCTTTTTCTTGTGGCAGGAACACGTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCAcgc                                                                <  2:2857663/93‑1 (MQ=255)
                         gAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGGACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCCGCCTCCcgcg                                                               >  2:3301245/1‑93 (MQ=255)
                         gAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCAcgcg                                                               >  1:6896312/1‑93 (MQ=255)
                          aaGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCTGCTGGTGGCATCACGCTCGCCTTCCAGCGCCGCCAGCATCACGCgg                                                              <  1:6456858/93‑1 (MQ=255)
                          aaGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCgg                                                              <  2:3905343/93‑1 (MQ=255)
                            gTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGc                                                            >  2:4383452/1‑93 (MQ=255)
                               aCCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGaa                                                         <  1:6968074/93‑1 (MQ=255)
                                 cTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGt                                                       >  1:3816069/1‑93 (MQ=255)
                                 cTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGt                                                       <  2:1563009/93‑1 (MQ=255)
                                  tGCCGACAACGCTTTTTCTTGTGGCAGGTACACTTTACCGCTGGTGGCCTCATGCTCGCCTTCCAGCGCTGCCAGCATCACGCTGGCGACGtt                                                      >  2:5654633/1‑93 (MQ=255)
                                  tGCCGACAACGCTTTTTCTTGTGGCAGGGACACCTTTCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGCAggt                                                      >  1:4144424/1‑91 (MQ=255)
                                  cgCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCAGCACGCGGGCGAAGtt                                                      <  2:3553191/92‑1 (MQ=255)
                                    ccGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGGAGGTGa                                                    >  1:1071962/1‑93 (MQ=255)
                                    ccGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGa                                                    <  1:2990524/93‑1 (MQ=255)
                                     cGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGGAGTTGAt                                                   >  1:4054582/1‑93 (MQ=255)
                                        cAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATaaa                                                <  2:5906028/93‑1 (MQ=255)
                                        cAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATaaa                                                <  1:5661708/93‑1 (MQ=255)
                                        cAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATaaa                                                <  1:7884224/93‑1 (MQ=255)
                                          aCGCTTTTTCTTGTGGCAGGAACACTTTTCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGc                                              >  1:2736572/1‑93 (MQ=255)
                                          aCGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGc                                              >  1:1597492/1‑93 (MQ=255)
                                          aCGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGc                                              <  2:7245268/93‑1 (MQ=255)
                                          aCGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGc                                              >  1:7309008/1‑93 (MQ=255)
                                          aCGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGGGGCATCATGCTCGCCTTTCAGCGCCGCCAGCATCACGCGGGCGAAGGTGATAAAGc                                              >  2:1498406/1‑93 (MQ=255)
                                          aCGCTTTTTCTTGTGGCAGGAAAACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGc                                              >  1:5414371/1‑93 (MQ=255)
                                           cGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCt                                             <  2:3967511/93‑1 (MQ=255)
                                            gCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAACTc                                            >  1:4696046/1‑93 (MQ=255)
                                            gCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGGTGATAAAGCTc                                            >  2:7022794/1‑93 (MQ=255)
                                             cTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCa                                           <  1:771510/93‑1 (MQ=255)
                                             cTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCa                                           <  2:3521235/93‑1 (MQ=255)
                                             cTTTTCCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGGAGCTGATAAAGCTCa                                           >  2:7369157/1‑93 (MQ=255)
                                              tttttCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCAt                                          <  1:3381803/93‑1 (MQ=255)
                                              tttttCTTGTCGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCAt                                          <  1:1682697/93‑1 (MQ=255)
                                                 ttCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGcc                                       <  2:4242097/93‑1 (MQ=255)
                                                  tCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCg                                      <  1:4283637/93‑1 (MQ=255)
                                                    ttGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGt                                    <  2:131624/93‑1 (MQ=255)
                                                      gtgGCAGGAACACTTTGCCGCTTGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTAc                                  >  1:1349933/1‑93 (MQ=255)
                                                      gtgGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTAc                                  <  2:7760146/93‑1 (MQ=255)
                                                       tgGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACa                                 <  1:2473452/93‑1 (MQ=255)
                                                       tgGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACa                                 <  1:5924197/93‑1 (MQ=255)
                                                       tgGCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACa                                 <  1:2847500/93‑1 (MQ=255)
                                                        ggCAGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAg                                >  1:7768154/1‑93 (MQ=255)
                                                           aGGAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGa                             >  2:1050695/1‑93 (MQ=255)
                                                            ggAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCgct                            >  1:3138377/1‑91 (MQ=255)
                                                            ggAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGAt                            <  1:1456341/93‑1 (MQ=255)
                                                             gAACACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATa                           >  1:3563173/1‑93 (MQ=255)
                                                              aaCACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAg                          >  1:1861599/1‑93 (MQ=255)
                                                              aaCACTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAg                          <  2:5586490/93‑1 (MQ=255)
                                                               acacTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGc                         <  2:2390526/93‑1 (MQ=255)
                                                               acacTTTGCCGCTGGTGGCATCATGCTCGCCGTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGGGATAGc                         >  1:2088628/1‑93 (MQ=255)
                                                                 acTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGccc                       >  1:746580/1‑93 (MQ=255)
                                                                 acTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGccc                       >  2:5309174/1‑93 (MQ=255)
                                                                  cTTTGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGcccc                      >  1:7399450/1‑93 (MQ=255)
                                                                    ttGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCAt                    <  1:822676/93‑1 (MQ=255)
                                                                     tGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCAtt                   <  1:3140231/93‑1 (MQ=255)
                                                                     tGCCGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCAtt                   <  2:4510469/93‑1 (MQ=255)
                                                                       ccGCTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTg                 <  1:4808029/93‑1 (MQ=255)
                                                                          cTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGcaa              >  1:1462000/1‑91 (MQ=255)
                                                                          cTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCa              >  1:3044900/1‑93 (MQ=255)
                                                                          cTGGTGGCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCa              >  2:4868868/1‑93 (MQ=255)
                                                                           tggtggCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCAc             <  2:512885/93‑1 (MQ=255)
                                                                             gtggCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCg           <  2:2635834/93‑1 (MQ=255)
                                                                             gtggCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCg           >  1:4658822/1‑93 (MQ=255)
                                                                              gggCATCATGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCCCATGCCTGTACAGCGATAGCCCCATTTGCCACCGa          <  1:6597897/92‑1 (MQ=255)
                                                                                 catcatGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGAcgg       <  1:1135188/93‑1 (MQ=255)
                                                                                 catcatGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGAcgg       <  2:6337982/93‑1 (MQ=255)
                                                                                 catcatGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGAcgg       <  1:1431412/93‑1 (MQ=255)
                                                                                  atcatGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGAcggc      <  2:7162875/93‑1 (MQ=255)
                                                                                   tcatGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACATCGATAGCCCCATTTGACACCGAcggcg     >  1:5056214/1‑93 (MQ=255)
                                                                                   tcatGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGAcggcg     >  2:5807703/1‑93 (MQ=255)
                                                                                   tcatGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGAaggcg     >  1:4539717/1‑93 (MQ=255)
                                                                                    catGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGAcggcgg    <  2:7879841/93‑1 (MQ=255)
                                                                                      tGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGACGGCGGtt  >  1:2283820/1‑93 (MQ=255)
                                                                                      tGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGACGGCGGtt  <  2:7558147/93‑1 (MQ=255)
                                                                                      tGCTCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGCTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGACGGCGGtt  >  1:5858986/1‑93 (MQ=255)
                                                                                         |                                                                                         
CGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAGCGCCGCCAGCATCACGCGGGCGAAGTTGATAAAGCTCATGCCGGTACAGCGATAGCCCCATTTGCCACCGACGGCGGTT  >  NC_000913/860323‑860501

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