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‑‑GGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTC  >  NC_000913/4475425‑4475686
                                                                                                                                    |                                                                                                                                     
aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                              <  2:995276/139‑2 (MQ=255)
aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                              <  2:511227/139‑2 (MQ=255)
aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCCGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gggggga                                                                                                                              >  1:897322/1‑138 (MQ=255)
   gAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gggggggcag                                                                                                                           >  1:1616710/1‑135 (MQ=255)
      ttCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGAt                                                                                                                        >  1:966721/1‑139 (MQ=255)
      accGTCATGAGCAACCAATCAGATTCAGATTATTCTCAACACATTACATTAGCCGACTTTTTTAAACAAAAAAGTGAAGTACTAACTGCCTTGCGTGTTGGAATGAATGATTTTGCCCGTATTACA‑‑GGTCGGGCAGGAt                                                                                                                        >  1:1588205/3‑139 (MQ=255)
        cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATaa                                                                                                                      >  2:858656/1‑139 (MQ=255)
        cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATaa                                                                                                                      >  2:792964/1‑139 (MQ=255)
        cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATaa                                                                                                                      >  2:1190853/1‑139 (MQ=255)
              gAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAg                                                                                                                >  1:535630/1‑139 (MQ=255)
                gCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAg                                                                                                                   <  1:1065217/134‑1 (MQ=255)
                gCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAg                                                                                                                   >  2:1065217/1‑134 (MQ=255)
                gCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGa                                                                                                              <  2:220831/139‑1 (MQ=255)
                gCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGa                                                                                                              <  1:792964/139‑1 (MQ=255)
                                aaaTTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gggggggca                                                                                                                            >  2:454295/1‑106 (MQ=255)
                                aaaTTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gggggggca                                                                                                                            <  1:454295/109‑4 (MQ=255)
                                  attattCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAg                                                                                            <  2:966721/139‑1 (MQ=255)
                                      ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGa                                                                                                              >  2:249982/1‑117 (MQ=255)
                                      ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGa                                                                                                              <  1:249982/117‑1 (MQ=255)
                                      ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCAtt                                                                                        >  1:746478/1‑139 (MQ=255)
                                          cAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTa                                                                                    <  2:1616710/139‑1 (MQ=255)
                                                cttacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGAATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGttt                                                                              >  1:1177511/2‑139 (MQ=255)
                                                      ttaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCa                                                                        <  2:535630/139‑1 (MQ=255)
                                                      ttaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCCGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑GGGGGGACAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCa                                                                        <  2:897322/139‑1 (MQ=255)
                                                       taGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGGATTTTTAGAGttttt                                                                            <  1:1255789/134‑1 (MQ=255)
                                                       taGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGttttt                                                                            >  2:1255789/1‑134 (MQ=255)
                                                         gCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAg                                                                     >  1:778749/1‑139 (MQ=255)
                                                            gACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGata                                                                  >  1:1559918/1‑139 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAg                                                                                                                   <  2:641138/88‑1 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAg                                                                                                                   >  1:641138/1‑88 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      >  1:496465/1‑133 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      <  2:496465/133‑1 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 >  1:71359/1‑138 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 <  2:71359/138‑1 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                              >  2:1835146/1‑76 (MQ=255)
                                                              cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                              <  1:1835146/77‑2 (MQ=255)
                                                               aaTTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGAATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATAttt                                                               <  2:1177511/139‑1 (MQ=255)
                                                                           aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                              <  1:1570029/64‑2 (MQ=255)
                                                                           aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                              >  2:1570029/1‑63 (MQ=255)
                                                                            aaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTAtg                                                  <  2:746478/139‑1 (MQ=255)
                                                                                  aaGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGAtt                                            <  1:1190853/139‑1 (MQ=255)
                                                                                   aGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGAttt                                           >  1:5700/1‑139 (MQ=255)
                                                                                              ccTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      >  1:575241/1‑101 (MQ=255)
                                                                                              ccTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      <  2:575241/101‑1 (MQ=255)
                                                                                               cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTc                               >  1:1359104/1‑139 (MQ=255)
                                                                                               cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                              >  2:1636405/1‑43 (MQ=255)
                                                                                               cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                              <  1:1636405/44‑2 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTa                                                                                    >  2:1876265/1‑85 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTa                                                                                    <  1:1876265/85‑1 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                              >  1:865956/1‑41 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                              >  2:364752/1‑41 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                              <  2:865956/42‑2 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                              <  1:364752/42‑2 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggg                                                                                                                                >  1:1031649/1‑40 (MQ=255)
                                                                                                 ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggg                                                                                                                                <  2:1031649/40‑1 (MQ=255)
                                                                                                       ttGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 >  1:416386/1‑97 (MQ=255)
                                                                                                       ttGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                 <  2:416386/97‑1 (MQ=255)
                                                                                                       ttGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTAtg                                                  >  1:303162/1‑112 (MQ=255)
                                                                                                       ttGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTAtg                                                  <  2:303162/112‑1 (MQ=255)
                                                                                                       ttGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAg                       >  2:131945/1‑139 (MQ=255)
                                                                                                       ttGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAg                       >  1:1794451/1‑139 (MQ=255)
                                                                                                       ttGGAATGAATGATTCTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAg                       <  1:742209/139‑1 (MQ=255)
                                                                                                           aatgaatgATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCggg                   <  1:131945/139‑1 (MQ=255)
                                                                                                            atgaatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAat                                                                                                            <  1:267493/49‑1 (MQ=39)
                                                                                                            atgaatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAat                                                                                                            >  2:267493/1‑49 (MQ=39)
                                                                                                               aatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      >  2:1795618/1‑84 (MQ=255)
                                                                                                               aatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      <  1:1795618/84‑1 (MQ=255)
                                                                                                               aatgaTTTTGCCCGTATTGCC‑‑GGGGGGCAGGATAACAGaa                                                                                                                 <  1:931064/40‑1 (MQ=255)
                                                                                                               aatgaTTTTGCCCGTATTGCC‑‑GGGGGGCAGGATAACAGaa                                                                                                                 >  2:931064/1‑40 (MQ=255)
                                                                                                                atgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGaa              >  2:1383228/1‑139 (MQ=255)
                                                                                                                    ttttGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTAcc                                                            <  2:571904/89‑1 (MQ=255)
                                                                                                                    ttttGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTAcc                                                            >  1:571904/1‑89 (MQ=255)
                                                                                                                         cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATgg     <  1:1383228/139‑1 (MQ=255)
                                                                                                                          ccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGc    <  2:778749/139‑1 (MQ=255)
                                                                                                                           cGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCt   <  2:1559918/139‑1 (MQ=255)
                                                                                                                            gTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTc  >  2:979964/1‑139 (MQ=255)
                                                                                                                             tATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATatg         >  1:1759331/1‑131 (MQ=255)
                                                                                                                             tATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATatg         <  2:1759331/131‑1 (MQ=255)
                                                                                                                              aTTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTa                                                                                    >  1:13221/1‑55 (MQ=255)
                                                                                                                              aTTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTa                                                                                    <  2:13221/55‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                     
AAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑GGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTC  >  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: