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.G95.1% 204.8 / 2.3 61coding (120/1815 nt)yjgLSopA‑central‑domain‑like hexapeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (30/28);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

AAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTC  >  NC_000913/4475425‑4475686
                                                                                                                                    |                                                                                                                                     
aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑ggggg                                                                                                                              <  2:995276/139‑2 (MQ=255)
aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑ggggg                                                                                                                              <  2:511227/139‑2 (MQ=255)
aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCCGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑gggga                                                                                                                              >  1:897322/1‑138 (MQ=255)
   gAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑gggggcag                                                                                                                           >  1:1616710/1‑135 (MQ=255)
      ttCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGAt                                                                                                                        >  1:966721/1‑139 (MQ=255)
      accGTCATGAGCAACCAATCAGATTCAGATTATTCTCAACACATTACATTAGCCGACTTTTTTAAACAAAAAAGTGAAGTACTAACTGCCTTGCGTGTTGGAATGAATGATTTTGCCCGTATTACA‑‑GGTCGGG‑CAGGAt                                                                                                                        >  1:1588205/3‑139 (MQ=255)
        cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATaa                                                                                                                      >  2:858656/1‑139 (MQ=255)
        cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATaa                                                                                                                      >  2:792964/1‑139 (MQ=255)
        cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATaa                                                                                                                      >  2:1190853/1‑139 (MQ=255)
              gAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAg                                                                                                                >  1:535630/1‑139 (MQ=255)
                gCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAg                                                                                                                   <  1:1065217/134‑1 (MQ=255)
                gCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAg                                                                                                                   >  2:1065217/1‑134 (MQ=255)
                gCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGa                                                                                                              <  2:220831/139‑1 (MQ=255)
                gCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGa                                                                                                              <  1:792964/139‑1 (MQ=255)
                                aaaTTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑gggggca                                                                                                                            >  2:454295/1‑106 (MQ=255)
                                aaaTTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑gggggca                                                                                                                            <  1:454295/109‑4 (MQ=255)
                                  attattCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAg                                                                                            <  2:966721/139‑1 (MQ=255)
                                      ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGa                                                                                                              >  2:249982/1‑117 (MQ=255)
                                      ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGa                                                                                                              <  1:249982/117‑1 (MQ=255)
                                      ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCAtt                                                                                        >  1:746478/1‑139 (MQ=255)
                                          cAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTa                                                                                    <  2:1616710/139‑1 (MQ=255)
                                                cttacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGAATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGttt                                                                              >  1:1177511/2‑139 (MQ=255)
                                                      ttaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCa                                                                        <  2:535630/139‑1 (MQ=255)
                                                      ttaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCCGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑GG‑GGGGACAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCa                                                                        <  2:897322/139‑1 (MQ=255)
                                                       taGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGGATTTTTAGAGttttt                                                                            <  1:1255789/134‑1 (MQ=255)
                                                       taGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGttttt                                                                            >  2:1255789/1‑134 (MQ=255)
                                                         gCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAg                                                                     >  1:778749/1‑139 (MQ=255)
                                                            gACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGata                                                                  >  1:1559918/1‑139 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAg                                                                                                                   <  2:641138/88‑1 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAg                                                                                                                   >  1:641138/1‑88 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      >  1:496465/1‑133 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      <  2:496465/133‑1 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 >  1:71359/1‑138 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 <  2:71359/138‑1 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑ggggg                                                                                                                              >  2:1835146/1‑76 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑ggggg                                                                                                                              <  1:1835146/77‑2 (MQ=255)
                                                               aaTTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGAATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATAttt                                                               <  2:1177511/139‑1 (MQ=255)
                                                                           aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑ggggg                                                                                                                              <  1:1570029/64‑2 (MQ=255)
                                                                           aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑ggggg                                                                                                                              >  2:1570029/1‑63 (MQ=255)
                                                                            aaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTAtg                                                  <  2:746478/139‑1 (MQ=255)
                                                                                  aaGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGAtt                                            <  1:1190853/139‑1 (MQ=255)
                                                                                   aGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGAttt                                           >  1:5700/1‑139 (MQ=255)
                                                                                              ccTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      >  1:575241/1‑101 (MQ=255)
                                                                                              ccTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      <  2:575241/101‑1 (MQ=255)
                                                                                               cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTc                               >  1:1359104/1‑139 (MQ=255)
                                                                                               cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑ggggg                                                                                                                              >  2:1636405/1‑43 (MQ=255)
                                                                                               cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑ggggg                                                                                                                              <  1:1636405/44‑2 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTa                                                                                    >  2:1876265/1‑85 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTa                                                                                    <  1:1876265/85‑1 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑ggggg                                                                                                                              >  1:865956/1‑41 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑ggggg                                                                                                                              >  2:364752/1‑41 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑ggggg                                                                                                                              <  2:865956/42‑2 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑ggggg                                                                                                                              <  1:364752/42‑2 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑ggg                                                                                                                                >  1:1031649/1‑40 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg‑ggg                                                                                                                                <  2:1031649/40‑1 (MQ=255)
                                                                                                       ttGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 >  1:416386/1‑97 (MQ=255)
                                                                                                       ttGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 <  2:416386/97‑1 (MQ=255)
                                                                                                       ttGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTAtg                                                  >  1:303162/1‑112 (MQ=255)
                                                                                                       ttGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTAtg                                                  <  2:303162/112‑1 (MQ=255)
                                                                                                       ttGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAg                       >  2:131945/1‑139 (MQ=255)
                                                                                                       ttGGAATGAATGATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAg                       >  1:1794451/1‑139 (MQ=255)
                                                                                                       ttGGAATGAATGATTCTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAg                       <  1:742209/139‑1 (MQ=255)
                                                                                                           aatgaatgATTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCggg                   <  1:131945/139‑1 (MQ=255)
                                                                                                            atgaatgaTTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAat                                                                                                            <  1:267493/49‑1 (MQ=39)
                                                                                                            atgaatgaTTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAat                                                                                                            >  2:267493/1‑49 (MQ=39)
                                                                                                               aatgaTTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      >  2:1795618/1‑84 (MQ=255)
                                                                                                               aatgaTTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      <  1:1795618/84‑1 (MQ=255)
                                                                                                               aatgaTTTTGCCCGTATTGCC‑‑GG‑GGGG‑CAGGATAACAGaa                                                                                                                 <  1:931064/40‑1 (MQ=255)
                                                                                                               aatgaTTTTGCCCGTATTGCC‑‑GG‑GGGG‑CAGGATAACAGaa                                                                                                                 >  2:931064/1‑40 (MQ=255)
                                                                                                                atgaTTTTGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGaa              >  2:1383228/1‑139 (MQ=255)
                                                                                                                    ttttGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTAcc                                                            <  2:571904/89‑1 (MQ=255)
                                                                                                                    ttttGCCCGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTAcc                                                            >  1:571904/1‑89 (MQ=255)
                                                                                                                         cccGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATgg     <  1:1383228/139‑1 (MQ=255)
                                                                                                                          ccGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGc    <  2:778749/139‑1 (MQ=255)
                                                                                                                           cGTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCt   <  2:1559918/139‑1 (MQ=255)
                                                                                                                            gTATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTc  >  2:979964/1‑139 (MQ=255)
                                                                                                                             tATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATatg         >  1:1759331/1‑131 (MQ=255)
                                                                                                                             tATTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATatg         <  2:1759331/131‑1 (MQ=255)
                                                                                                                              aTTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTa                                                                                    >  1:13221/1‑55 (MQ=255)
                                                                                                                              aTTGCCG‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTa                                                                                    <  2:13221/55‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                     
AAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑GG‑GGGG‑CAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTC  >  NC_000913/4475425‑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: ‑