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.G98.5% 231.5 / ‑3.6 66coding (120/1815 nt)yjgLSopA‑central‑domain‑like hexapeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (31/34);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTC  >  NC_000913/4475423‑4475686
                                                                                                                                       |                                                                                                                                  
tCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggg                                                                                                                               >  2:995084/1‑139 (MQ=255)
tCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggg                                                                                                                               >  1:1356565/1‑139 (MQ=255)
tCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggg                                                                                                                               >  2:517314/1‑139 (MQ=255)
tCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggg                                                                                                                               <  2:1356565/139‑1 (MQ=255)
tCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggg                                                                                                                               >  2:301026/1‑139 (MQ=255)
tCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCgg                                                                                                                                  <  2:404734/136‑1 (MQ=255)
tCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCgg                                                                                                                                  >  1:404734/1‑136 (MQ=255)
  aaGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggggg                                                                                                                             <  2:1217778/139‑2 (MQ=255)
        ttCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGAt                                                                                                                        <  1:301026/139‑1 (MQ=255)
          cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATaa                                                                                                                      >  1:1424052/1‑139 (MQ=255)
                gAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGACTGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAg                                                                                                                >  2:401881/1‑139 (MQ=255)
                gAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggggg                                                                                                                             >  2:1605505/1‑124 (MQ=255)
                gAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCggggggg                                                                                                                             <  1:1605505/125‑2 (MQ=255)
                       atatCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGa                                                                                                                         <  1:1612411/123‑1 (MQ=255)
                       atatCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGa                                                                                                                         >  2:1612411/1‑123 (MQ=255)
                           cAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATAttcctt                                                                                                     >  2:593449/1‑139 (MQ=255)
                                    attattCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAg                                                                                            <  2:466214/139‑1 (MQ=255)
                                        ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCAtt                                                                                        >  2:887913/1‑139 (MQ=255)
                                        ttCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCAtt                                                                                        >  1:829254/1‑139 (MQ=255)
                                                  attacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGttt                                                                              >  1:996989/1‑139 (MQ=255)
                                                    tacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCgggggggcagg                                                                                                                         <  2:804316/93‑6 (MQ=255)
                                                    tacattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCgggggggcagg                                                                                                                         >  1:804316/1‑88 (MQ=255)
                                                     acattaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTg                                                                           >  2:1069881/1‑139 (MQ=255)
                                                        ttaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCa                                                                        <  2:829254/139‑1 (MQ=255)
                                                        ttaGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCa                                                                        <  2:929547/139‑1 (MQ=255)
                                                         taGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaa                                                                       >  1:202951/1‑139 (MQ=255)
                                                          aGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      <  2:996989/139‑1 (MQ=255)
                                                          aGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      <  2:699882/139‑1 (MQ=255)
                                                          aGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      <  2:202951/139‑1 (MQ=255)
                                                          aGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaaa                                                                      >  1:699882/1‑139 (MQ=255)
                                                          aGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCgg                                                                                                                                  >  2:911494/1‑78 (MQ=255)
                                                          aGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCgg                                                                                                                                  <  1:911494/78‑1 (MQ=255)
                                                           gCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAg                                                                     >  2:276044/1‑139 (MQ=255)
                                                            ccGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGa                                                                    >  1:677017/1‑139 (MQ=255)
                                                                cAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATAtt                                                                <  2:1600952/139‑1 (MQ=255)
                                                                   ttttAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTAc                                                             <  1:1069881/139‑1 (MQ=255)
                                                                        aaCAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTgg                                                        <  1:887913/139‑1 (MQ=255)
                                                                             aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTAcc                                                            <  1:263903/130‑1 (MQ=255)
                                                                             aaaaGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTAcc                                                            >  2:263903/1‑130 (MQ=255)
                                                                                            aaTACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGaa                                                                                                             <  2:52247/66‑1 (MQ=255)
                                                                                            aaTACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGaa                                                                                                             >  1:52247/1‑66 (MQ=255)
                                                                                                 cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGAt                                                                                                                        <  2:402444/50‑1 (MQ=255)
                                                                                                 cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGAt                                                                                                                        >  1:402444/1‑50 (MQ=255)
                                                                                                 cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATAttct                                                                                                       >  2:622656/1‑67 (MQ=255)
                                                                                                 cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATAttct                                                                                                       <  1:622656/67‑1 (MQ=255)
                                                                                                 cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCa                                                                        >  2:855493/1‑98 (MQ=255)
                                                                                                 cTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCa                                                                        <  1:855493/98‑1 (MQ=255)
                                                                                                   ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGAt                                                                                                                        >  1:688119/1‑48 (MQ=255)
                                                                                                   ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGAt                                                                                                                        <  2:688119/48‑1 (MQ=39)
                                                                                                   ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAat                                                                                                            >  1:812733/1‑60 (MQ=255)
                                                                                                   ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAat                                                                                                            <  2:812733/60‑1 (MQ=255)
                                                                                                   ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTAtg                                                  >  1:126571/1‑118 (MQ=255)
                                                                                                   ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTAtg                                                  <  2:126571/118‑1 (MQ=255)
                                                                                                   ggCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCa                             >  2:719680/1‑139 (MQ=255)
                                                                                                     cGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAAc                           >  2:47411/1‑139 (MQ=255)
                                                                                                         ttGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCCTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAg                       <  1:593449/139‑1 (MQ=255)
                                                                                                          tGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATAttcttt                                                                                                     <  1:245129/60‑1 (MQ=255)
                                                                                                          tGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATAttcttt                                                                                                     >  2:245129/1‑60 (MQ=255)
                                                                                                          tGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTAcc                                                            <  1:296116/101‑1 (MQ=255)
                                                                                                          tGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTAcc                                                            >  2:296116/1‑101 (MQ=255)
                                                                                                          tGGAATGAATGATTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGc                      >  2:517442/1‑139 (MQ=255)
                                                                                                              atgaatgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGt                  <  1:276044/139‑1 (MQ=255)
                                                                                                                   tgaTTTTGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGaaa             >  1:1442616/1‑139 (MQ=255)
                                                                                                                       tttGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATatg         <  2:1424052/139‑1 (MQ=255)
                                                                                                                        ttGCCCGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATatga        <  2:1442616/139‑1 (MQ=255)
                                                                                                                           cccGTATTGTCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCaaa                                    >  1:300721/1‑108 (MQ=255)
                                                                                                                           cccGTATTGTCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCaaa                                    <  2:300721/108‑1 (MQ=255)
                                                                                                                           cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTa                                                                                    >  2:969109/1‑60 (MQ=255)
                                                                                                                           cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTa                                                                                    <  1:969109/60‑1 (MQ=255)
                                                                                                                           cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCa                                      >  1:804448/1‑106 (MQ=255)
                                                                                                                           cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCa                                      <  2:804448/106‑1 (MQ=255)
                                                                                                                           cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGc                      >  1:334816/1‑122 (MQ=255)
                                                                                                                           cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGc                      <  2:334816/122‑1 (MQ=255)
                                                                                                                           cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATgg     <  2:677017/139‑1 (MQ=255)
                                                                                                                           cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATgg     >  1:23062/1‑139 (MQ=255)
                                                                                                                           cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATgg     <  1:47411/139‑1 (MQ=255)
                                                                                                                           cccGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATgg     <  1:719680/139‑1 (MQ=255)
                                                                                                                             cGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCaaa                                    >  1:1522830/1‑106 (MQ=255)
                                                                                                                             cGTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCaaa                                    <  2:1522830/106‑1 (MQ=255)
                                                                                                                              gTATTGCCGGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTc  <  1:995084/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                  
TCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTC  >  NC_000913/4475423‑4475686

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