Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,116,990 T→C 5.8% Q708R (CAG→CGG)  sslE ← putative secreted and surface‑associated lipoprotein mucinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,116,9900TC 5.8% 224.2 / 4.3 69Q708R (CAG→CGG) sslEputative secreted and surface‑associated lipoprotein mucinase
Reads supporting (aligned to +/- strand):  ref base T (25/40);  new base C (2/2);  total (27/42)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

AGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTA  >  NC_000913/3116864‑3117124
                                                                                                                              |                                                                                                                                      
aGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCtt                                                                                                                           <  2:61581/140‑1 (MQ=255)
aGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCtt                                                                                                                           >  2:1179495/1‑140 (MQ=255)
     tGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGtt                                                                                                                      >  1:374692/1‑140 (MQ=255)
              cAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTAt                                                                                                             >  1:177012/1‑140 (MQ=255)
                 cAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTg                                                                                                          >  2:1330169/1‑140 (MQ=255)
                   gACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCa                                                                                                        <  2:1074290/140‑1 (MQ=255)
                             aCCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGa                                                                                              <  1:82884/140‑1 (MQ=255)
                             aCCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGa                                                                                              <  1:1192268/140‑1 (MQ=255)
                             aCCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGa                                                                                              <  1:1237621/140‑1 (MQ=255)
                             aCCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGa                                                                                              <  2:356562/140‑1 (MQ=255)
                             aCCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGa                                                                                              <  1:1331673/140‑1 (MQ=255)
                             aCCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGa                                                                                              <  2:1577637/140‑1 (MQ=255)
                             aCCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGa                                                                                              <  2:1595264/140‑1 (MQ=255)
                                     cgcTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATaa                                                                                      <  1:3585/140‑1 (MQ=255)
                                      gcTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAg                                                                                     >  2:1621916/1‑140 (MQ=255)
                                      gcTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAg                                                                                     >  1:903468/1‑140 (MQ=255)
                                      gcTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAg                                                                                     >  2:1621770/1‑140 (MQ=255)
                                      gcTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCCGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGa                                                                                              >  2:484164/1‑131 (MQ=255)
                                      gcTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCCGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGa                                                                                              >  2:484764/1‑131 (MQ=255)
                                      gcTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCCGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGa                                                                                              <  1:484164/131‑1 (MQ=255)
                                      gcTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCCGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGa                                                                                              <  1:484764/131‑1 (MQ=255)
                                        tctctTCCTTAAGTTTGCTCCTCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGAAAGGCGATACGATCGGTCACATCCTGCTTAGTCAGCTTCGGTTTTCTGGAATCTGAACGAAGCGGGATTGCATTAgta                                                                                   >  2:751739/1‑138 (MQ=255)
                                         ctctTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGAt                                                                                  <  1:1558117/140‑1 (MQ=255)
                                         ctctTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGAt                                                                                  <  2:325761/140‑1 (MQ=255)
                                         ctctTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGAt                                                                                  <  2:327207/140‑1 (MQ=255)
                                         ctctTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGAt                                                                                  <  2:452662/140‑1 (MQ=255)
                                         ctctTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGAt                                                                                  <  2:64469/140‑1 (MQ=255)
                                         ctctTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGAt                                                                                  <  2:73450/140‑1 (MQ=255)
                                         ctctTCCTTAAGATTGCTCATCACGATTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGAt                                                                                  <  2:937662/140‑1 (MQ=255)
                                         ctcctcCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGAt                                                                                  <  2:327093/140‑1 (MQ=255)
                                             tCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGtt                                                                              <  1:740627/140‑1 (MQ=255)
                                             tCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGtt                                                                              <  1:441892/140‑1 (MQ=255)
                                                ttAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTAcc                                                                           <  1:1311998/140‑1 (MQ=255)
                                                   aGATTGCTCATCACGATTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCAc                                                                        <  1:613966/140‑1 (MQ=255)
                                                   aGATTGCTCATCACGATTTCCATGATCAGCACCGATCCACCTATGTTCAGATAGGCGATCAGATCGGACACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCAc                                                                        <  2:1135095/140‑1 (MQ=255)
                                                    gATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCAcc                                                                       >  2:1076414/1‑140 (MQ=255)
                                                       tGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGa                                                                    <  1:1011338/140‑1 (MQ=255)
                                                        gCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAg                                                                   >  1:1267351/1‑140 (MQ=255)
                                                        gCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAg                                                                   >  1:28541/1‑140 (MQ=255)
                                                         cTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGt                                                                  >  2:262806/1‑140 (MQ=255)
                                                                   tttCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCa                                                        >  1:366683/1‑140 (MQ=255)
                                                                   tttCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCa                                                        >  1:366684/1‑140 (MQ=255)
                                                                    ttCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCaa                                                       <  1:1436169/140‑1 (MQ=255)
                                                                     tCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCaaa                                                      <  1:1621916/140‑1 (MQ=255)
                                                                     tCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCaaa                                                      <  1:1474782/140‑1 (MQ=255)
                                                                     tCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCaaa                                                      <  1:1621770/140‑1 (MQ=255)
                                                                       cATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGc                                                    <  2:1294326/140‑1 (MQ=255)
                                                                       cATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGc                                                    <  2:1294234/140‑1 (MQ=255)
                                                                          gATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGt                                                 <  1:1179495/140‑1 (MQ=255)
                                                                             cAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTaa                                              <  2:465797/140‑1 (MQ=255)
                                                                             cAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTaa                                              <  2:1236306/140‑1 (MQ=255)
                                                                                     aTCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATcagc                                      >  1:585217/1‑140 (MQ=255)
                                                                                     aTCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATcagc                                      >  1:1227787/1‑140 (MQ=255)
                                                                                     aTCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATcagc                                      >  1:1227768/1‑140 (MQ=255)
                                                                                      tCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATcagca                                     <  1:1610363/140‑1 (MQ=255)
                                                                                          ccTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCgg                                 >  2:502186/1‑140 (MQ=255)
                                                                                          ccTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCgg                                 >  2:502434/1‑140 (MQ=255)
                                                                                          ccTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCgg                                 >  2:1264307/1‑140 (MQ=255)
                                                                                             ttGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAt                              >  1:1135650/1‑140 (MQ=255)
                                                                                                ttCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTc                           >  1:102002/1‑140 (MQ=255)
                                                                                                          cGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCg                 >  2:1028974/1‑140 (MQ=255)
                                                                                                            aTCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCgag               <  2:670574/140‑1 (MQ=255)
                                                                                                                  tCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGcc         <  1:1330169/140‑1 (MQ=255)
                                                                                                                   cGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCa        >  2:1150400/1‑140 (MQ=255)
                                                                                                                   cGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCa        >  2:1080029/1‑140 (MQ=255)
                                                                                                                       cacaTCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGc    <  1:1309558/140‑1 (MQ=255)
                                                                                                                       cacaTCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGc    <  2:1614900/140‑1 (MQ=255)
                                                                                                                       cacaTCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGc    <  2:1596764/140‑1 (MQ=255)
                                                                                                                         caTCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTa  <  2:241234/140‑1 (MQ=255)
                                                                                                                              |                                                                                                                                      
AGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATCAGCACCGATCCACCTTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTA  >  NC_000913/3116864‑3117124

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