Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,475,562 (G)6→7 coding (126/1815 nt) yjgL → SopA‑central‑domain‑like hexapeptide repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,475,5561.G97.0% 230.6 / ‑0.6 67coding (120/1815 nt)yjgLSopA‑central‑domain‑like hexapeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (30/35);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

CACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTC  >  NC_000913/4475420‑4475686
                                                                                                                                         |                                                                                                                                  
cACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCgg                                                                                                                                  <  1:1271944/139‑1 (MQ=255)
 aCTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggg                                                                                                                                 <  2:167409/139‑1 (MQ=255)
   tCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggg                                                                                                                               <  1:1351295/139‑1 (MQ=255)
   tCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggg                                                                                                                               >  2:1351295/1‑139 (MQ=255)
   tCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggg                                                                                                                               >  1:1269403/1‑139 (MQ=255)
     aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggggg                                                                                                                             <  2:747723/139‑2 (MQ=255)
          gTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGa                                                                                                                         <  1:1501652/139‑1 (MQ=255)
             cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATaa                                                                                                                      >  1:285195/1‑139 (MQ=255)
             cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATaa                                                                                                                      >  2:577211/1‑139 (MQ=255)
                 aTGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCATATTGCCGGGGGGGCAGGATAACAGa                                                                                                                  <  2:983018/139‑1 (MQ=255)
                  tGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGa                                                                                                                         >  2:369597/1‑131 (MQ=255)
                  tGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGa                                                                                                                         <  1:369597/131‑1 (MQ=255)
                    aGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggg                                                                                                                               >  2:1521198/1‑122 (MQ=255)
                    aGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggg                                                                                                                               <  1:1521198/122‑1 (MQ=255)
                     gCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggggg                                                                                                                             <  1:1280906/123‑2 (MQ=255)
                     gCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggggg                                                                                                                             >  2:1280906/1‑122 (MQ=255)
                              cAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATAttcttt                                                                                                     >  2:147117/1‑139 (MQ=255)
                                   ttAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCt                                                                                                <  1:1460774/139‑1 (MQ=255)
                                       attattCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAg                                                                                            <  1:1243919/139‑1 (MQ=255)
                                           ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCt                                                                                                <  2:71129/131‑1 (MQ=255)
                                           ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCt                                                                                                >  1:71129/1‑131 (MQ=255)
                                           ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCAtt                                                                                        >  1:653359/1‑139 (MQ=255)
                                               cAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTa                                                                                    <  1:147117/139‑1 (MQ=255)
                                                aaCACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCAt                                                                                         <  1:1469426/133‑1 (MQ=255)
                                                aaCACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCAt                                                                                         >  2:1469426/1‑133 (MQ=255)
                                                     attacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGttt                                                                              >  1:546209/1‑139 (MQ=255)
                                                              gCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGAt                                                                                                                        >  2:1145300/1‑88 (MQ=255)
                                                              gCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGAt                                                                                                                        <  1:1145300/88‑1 (MQ=255)
                                                               ccGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGa                                                                    >  2:1225180/1‑139 (MQ=255)
                                                               ccGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGa                                                                    >  2:1251368/1‑139 (MQ=255)
                                                               ccGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggggg                                                                                                                             <  1:1300983/81‑2 (MQ=255)
                                                               ccGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggggg                                                                                                                             >  2:1300983/1‑80 (MQ=255)
                                                               ccGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATATTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGa                                                                    >  1:34711/1‑139 (MQ=255)
                                                                 gACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGata                                                                  >  1:660611/1‑139 (MQ=255)
                                                                  aCAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 <  2:53648/139‑1 (MQ=255)
                                                                  aCAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 <  1:64512/139‑1 (MQ=255)
                                                                  aCAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 <  1:1251368/139‑1 (MQ=255)
                                                                   cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGa                                                                                                                  >  1:259865/1‑89 (MQ=255)
                                                                   cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGa                                                                                                                  <  2:259865/89‑1 (MQ=255)
                                                                   cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggggg                                                                                                                             <  1:246611/77‑2 (MQ=255)
                                                                   cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggggg                                                                                                                             >  2:246611/1‑76 (MQ=255)
                                                                              aaaaaaGTGAAGTTTTAAATACCTGGAGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGtt                                                     >  1:384891/1‑139 (MQ=255)
                                                                                aaaaGTGGAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCgggggggc                                                                                                                            >  1:1055610/1‑63 (MQ=255)
                                                                                aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGAt                                                                                                                        <  1:1273788/70‑1 (MQ=255)
                                                                                aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGAt                                                                                                                        >  2:1273788/1‑70 (MQ=255)
                                                                                aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCgggggggc                                                                                                                            <  2:1055610/65‑3 (MQ=255)
                                                                                aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggggg                                                                                                                             <  1:1050243/64‑2 (MQ=255)
                                                                                aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggggg                                                                                                                             >  2:1050243/1‑63 (MQ=255)
                                                                                   aGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTAtgtg                                                <  2:1269403/139‑1 (MQ=255)
                                                                                               aaTACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCaaa                                    <  1:1022012/139‑1 (MQ=255)
                                                                                               aaTACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCaaa                                    <  2:546209/139‑1 (MQ=255)
                                                                                               aaTACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCaaa                                    <  2:285195/139‑1 (MQ=255)
                                                                                                  aCCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGAt                                                                                                                        >  1:1567251/1‑52 (MQ=255)
                                                                                                  aCCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGAt                                                                                                                        <  2:1567251/52‑1 (MQ=255)
                                                                                                    cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAat                                                                                                            <  2:1081921/61‑1 (MQ=255)
                                                                                                    cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAat                                                                                                            >  1:1081921/1‑61 (MQ=255)
                                                                                                      ggCGTGTTGGAATGCATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 >  2:70399/1‑103 (MQ=255)
                                                                                                      ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCa                             >  1:961500/1‑139 (MQ=255)
                                                                                                      ggCGTGTTGGAATGAATGATGTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 <  1:70399/103‑1 (MQ=255)
                                                                                                           gtTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGaa                        <  1:680143/139‑1 (MQ=255)
                                                                                                           gtTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGaa                        <  2:384891/139‑1 (MQ=255)
                                                                                                            ttGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGa                                                                                                                         <  1:519050/41‑1 (MQ=255)
                                                                                                            ttGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGa                                                                                                                         >  2:519050/1‑41 (MQ=255)
                                                                                                            ttGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAata                                                                                                           >  2:908282/1‑55 (MQ=255)
                                                                                                            ttGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAata                                                                                                           <  1:908282/55‑1 (MQ=255)
                                                                                                                aatgaatgATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCggg                   <  2:653359/139‑1 (MQ=255)
                                                                                                                 atgaatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGt                  >  1:889118/1‑139 (MQ=255)
                                                                                                                 atgaatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGt                  <  2:1390053/139‑1 (MQ=255)
                                                                                                                 atgaatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGt                  >  1:1390053/1‑139 (MQ=255)
                                                                                                                 atgaatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGt                  >  1:296967/1‑139 (MQ=255)
                                                                                                                 atgaatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGt                  >  2:537460/1‑139 (MQ=255)
                                                                                                                         ttttGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCa                                                                        >  1:24885/1‑77 (MQ=255)
                                                                                                                         ttttGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCa                                                                        <  2:24885/77‑1 (MQ=255)
                                                                                                                          tttGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 >  1:1466037/1‑83 (MQ=255)
                                                                                                                          tttGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 <  2:1466037/83‑1 (MQ=255)
                                                                                                                          tttGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATatg         <  1:358481/139‑1 (MQ=255)
                                                                                                                           ttGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATatga        <  1:1225180/139‑1 (MQ=255)
                                                                                                                           ttGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATatga        <  2:889118/139‑1 (MQ=255)
                                                                                                                              cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCt                                                          <  1:173702/86‑1 (MQ=255)
                                                                                                                              cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCt                                                          >  2:173702/1‑86 (MQ=255)
                                                                                                                              cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCaaa                                    <  1:205681/108‑1 (MQ=255)
                                                                                                                              cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCaaa                                    >  2:205681/1‑108 (MQ=255)
                                                                                                                              cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATgg     >  2:573518/1‑139 (MQ=255)
                                                                                                                                 gTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTc  <  2:660611/139‑1 (MQ=255)
                                                                                                                                 gTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATAGTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTc  <  2:34711/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                  
CACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTC  >  NC_000913/4475420‑4475686

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