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.G96.0% 254.7 / 2.5 74coding (120/1815 nt)yjgLSopA‑central‑domain‑like hexapeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base G (37/34);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

ACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑GGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGG  >  NC_000913/4475421‑4475683
                                                                                                                                         |                                                                                                                                
aCTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggg                                                                                                                              >  1:649305/1‑139 (MQ=255)
  tCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggg                                                                                                                            >  2:432963/1‑139 (MQ=255)
    aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                          <  1:432963/139‑2 (MQ=255)
    aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                          <  1:1230411/139‑2 (MQ=255)
    aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                          >  1:199646/1‑138 (MQ=255)
            cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATaa                                                                                                                   >  2:185410/1‑139 (MQ=255)
            cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATaa                                                                                                                   >  2:1563887/1‑139 (MQ=255)
               cATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggg                                                                                                                            <  1:1004235/126‑1 (MQ=255)
               cATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGCGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggg                                                                                                                            >  2:1004235/1‑126 (MQ=255)
                    gCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGa                                                                                                           <  2:649305/139‑1 (MQ=255)
                           atCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATAttct                                                                                                    <  1:1563887/139‑1 (MQ=255)
                                       ttattCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGc                                                                                        <  2:538535/139‑1 (MQ=255)
                                          ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCAtt                                                                                     >  2:388576/1‑139 (MQ=255)
                                          ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCAtt                                                                                     >  1:1367921/1‑139 (MQ=255)
                                          ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gggggggca                                                                                                                        >  1:54657/1‑100 (MQ=255)
                                          ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gggggggca                                                                                                                        <  2:54657/103‑4 (MQ=255)
                                               aaCACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTag                                                                                >  2:1423907/1‑139 (MQ=255)
                                                    attacattcGCCAACAATTATAAACAAAAAAGTGAAGTCCTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGTAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGttt                                                                           <  1:459469/139‑1 (MQ=255)
                                                    attacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGttt                                                                           >  1:523969/1‑139 (MQ=255)
                                                      tacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAg                                                                                                                <  2:1520572/100‑1 (MQ=255)
                                                      tacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAg                                                                                                                >  1:1520572/1‑100 (MQ=255)
                                                      tacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGttttt                                                                         >  1:1747366/1‑139 (MQ=255)
                                                      tacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGttttt                                                                         >  1:851879/1‑139 (MQ=255)
                                                      tacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggg                                                                                                                            >  1:1578846/1‑87 (MQ=255)
                                                      tacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggg                                                                                                                            <  2:1578846/87‑1 (MQ=255)
                                                         attaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATaa                                                                                                                   >  1:1464594/1‑94 (MQ=255)
                                                         attaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATaa                                                                                                                   <  2:1464594/94‑1 (MQ=255)
                                                          ttaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCa                                                                     <  2:1367921/139‑1 (MQ=255)
                                                           taGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaa                                                                    >  1:590592/1‑139 (MQ=255)
                                                             gCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAg                                                                                                                <  1:1651505/93‑1 (MQ=255)
                                                             gCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAg                                                                                                                >  2:1651505/1‑93 (MQ=255)
                                                             gCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAg                                                                  >  2:1045888/1‑139 (MQ=255)
                                                             gCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAg                                                                  >  2:1431230/1‑139 (MQ=255)
                                                             gCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg                                                                                                                               <  2:473619/77‑1 (MQ=255)
                                                             gCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑gg                                                                                                                               >  1:473619/1‑77 (MQ=255)
                                                              ccGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                   <  1:815555/137‑1 (MQ=255)
                                                              ccGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                   >  2:815555/1‑137 (MQ=255)
                                                                gACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                   <  1:821303/135‑1 (MQ=255)
                                                                gACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                   >  2:821303/1‑135 (MQ=255)
                                                                 aCAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGTAAAGatat                                                              <  2:1747366/139‑1 (MQ=255)
                                                                 aCAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                              <  2:590592/139‑1 (MQ=255)
                                                                  cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGACG‑‑GGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                              >  2:1141523/1‑138 (MQ=255)
                                                                  cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAgg                                                                                                            >  2:416495/1‑92 (MQ=255)
                                                                  cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAgg                                                                                                            <  1:416495/92‑1 (MQ=255)
                                                                  cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                              <  1:1141523/138‑1 (MQ=255)
                                                                  cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                          >  2:472863/1‑76 (MQ=255)
                                                                  cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                          <  1:472863/77‑2 (MQ=255)
                                                                   aaTTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                          <  2:209787/76‑2 (MQ=255)
                                                                   aaTTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑ggggggg                                                                                                                          >  1:209787/1‑75 (MQ=255)
                                                                               aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTAt                                                >  2:1005588/1‑139 (MQ=255)
                                                                                           ttAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAatatc                                                                                                      >  1:631373/1‑72 (MQ=255)
                                                                                           ttAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATAtt                                                                                                      <  2:631373/73‑1 (MQ=255)
                                                                                              aaTACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGa                                                                                                               <  1:587382/61‑1 (MQ=255)
                                                                                              aaTACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGa                                                                                                               >  2:587382/1‑61 (MQ=255)
                                                                                                tACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACg                               <  1:185410/139‑1 (MQ=255)
                                                                                                  ccTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCt                             >  2:491795/1‑139 (MQ=255)
                                                                                                   cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAg                                                                                                                <  1:271679/55‑1 (MQ=255)
                                                                                                   cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAg                                                                                                                >  2:271679/1‑55 (MQ=255)
                                                                                                   cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTc                            >  1:1520097/1‑139 (MQ=255)
                                                                                                     ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGa                                                                                                               >  2:1085032/1‑54 (MQ=255)
                                                                                                     ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGa                                                                                                               <  1:1085032/54‑1 (MQ=255)
                                                                                                     ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAg                                                                                                             <  1:567634/56‑1 (MQ=255)
                                                                                                     ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAg                                                                                                             >  2:567634/1‑56 (MQ=255)
                                                                                                          gtTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                              >  2:199577/1‑98 (MQ=255)
                                                                                                          gtTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                              <  1:199577/98‑1 (MQ=255)
                                                                                                          gtTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGaa                     <  1:388576/139‑1 (MQ=255)
                                                                                                             ggAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCg                  <  1:8560/139‑1 (MQ=255)
                                                                                                               aatgaatgATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCggg                <  2:851879/139‑1 (MQ=255)
                                                                                                                   aatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAg                                                                                         <  1:1016060/62‑1 (MQ=255)
                                                                                                                   aatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAg                                                                                         >  2:1016060/1‑62 (MQ=255)
                                                                                                                   aatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCaaa                                 >  1:950148/1‑118 (MQ=255)
                                                                                                                   aatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCaaa                                 <  2:950148/118‑1 (MQ=255)
                                                                                                                   aatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGa            <  2:199646/139‑1 (MQ=255)
                                                                                                                   aatgaTTTTGCCCGTATTGCC‑‑GGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACg                               >  1:1453024/1‑119 (MQ=255)
                                                                                                                   aatgaTTTTGCCCGTATTGCC‑‑GGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTCTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACg                               <  2:1453024/119‑1 (MQ=255)
                                                                                                                    atgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGaa           >  2:1652506/1‑139 (MQ=255)
                                                                                                                    atgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTCTCCTGGAg                                                                                         <  2:1107771/60‑1 (MQ=255)
                                                                                                                    atgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTCTCCTGGAg                                                                                         >  1:1107771/1‑60 (MQ=255)
                                                                                                                     tgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                              >  2:1198887/1‑87 (MQ=255)
                                                                                                                     tgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                              <  1:1198887/87‑1 (MQ=255)
                                                                                                                           tGCCCGTATTGCCGGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATatga     <  1:1399113/139‑1 (MQ=255)
                                                                                                                           tGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATatgat    >  2:438625/1‑139 (MQ=255)
                                                                                                                             cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCa                                                                                       <  2:998111/54‑1 (MQ=255)
                                                                                                                             cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCa                                                                                       >  1:998111/1‑54 (MQ=255)
                                                                                                                             cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCa                                                                     <  2:309398/72‑1 (MQ=255)
                                                                                                                             cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCa                                                                     >  1:309398/1‑72 (MQ=255)
                                                                                                                             cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGt               <  2:894277/126‑1 (MQ=255)
                                                                                                                             cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGt               >  1:894277/1‑126 (MQ=255)
                                                                                                                             cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATgg  >  1:633992/1‑139 (MQ=255)
                                                                                                                                  aTTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCt                             <  1:1525850/107‑1 (MQ=255)
                                                                                                                                  aTTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCt                             >  2:1525850/1‑107 (MQ=255)
                                                                                                                                         |                                                                                                                                
ACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCC‑‑GGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGG  >  NC_000913/4475421‑4475683

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