Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,580,434 (G)5→4 coding (357/1158 nt) ydeM ← putative anaerobic sulfatase maturation enzyme YdeM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,580,4300G.98.7% 341.1 / ‑2.2 77coding (361/1158 nt)ydeMputative anaerobic sulfatase maturation enzyme YdeM
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (35/41);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGAT  >  NC_000913/1580291‑1580561
                                                                                                                                           |                                                                                                                                   
tAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTg                                                                                                                                     >  1:1035017/1‑140 (MQ=255)
tAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTg                                                                                                                                     >  2:856030/1‑140 (MQ=255)
tAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTg                                                                                                                                     >  1:1461863/1‑140 (MQ=255)
      aCATTAACGTTATTAATGACGGGTAACGTATTAAACTCTACTTGATATGATTGCAGGCGCTCGAGGGCTGCTATCTCTTTGGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTggggcca                                                                                                                               <  2:1223988/140‑4 (MQ=255)
        aTTAACGTTATTAATGACGGTTAACGTATTAAACTCTGCTTGATATGATTTCAGACGCTCGATGGCTGCTTTCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCatc                                                                                                                            <  1:1171403/140‑1 (MQ=255)
        aTTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCatc                                                                                                                            <  2:1416070/140‑1 (MQ=255)
        aTTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCatc                                                                                                                            <  2:556709/140‑1 (MQ=255)
        aTTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCatc                                                                                                                            >  1:260016/1‑140 (MQ=255)
        aTTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCatc                                                                                                                            <  2:75766/140‑1 (MQ=255)
        aTTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCatc                                                                                                                            >  1:1397194/1‑140 (MQ=255)
          tAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCatcga                                                                                                                          >  2:1182520/1‑140 (MQ=255)
          tAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCatcga                                                                                                                          >  2:747804/1‑140 (MQ=255)
          tAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCatcga                                                                                                                          <  2:902459/140‑1 (MQ=255)
           aaCGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCatcgat                                                                                                                         <  1:391080/140‑1 (MQ=255)
                 attaATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGAt                                                                                                                   >  1:563170/1‑140 (MQ=255)
                 attaATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGAt                                                                                                                   >  1:1070034/1‑140 (MQ=255)
                  ttaATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATa                                                                                                                  >  2:805749/1‑140 (MQ=255)
                    aaTGACGGTTAACGTATTAATCTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATacc                                                                                                                <  1:1228438/140‑1 (MQ=255)
                      tGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATaccaa                                                                                                              >  2:1522419/1‑140 (MQ=255)
                        aCGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATaccaacc                                                                                                            >  1:110144/1‑140 (MQ=255)
                        aCGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGCTCATGTAACTCCTGGGGGCCATCGATCGAGATaccaac                                                                                                             >  2:1536714/1‑140 (MQ=255)
                            ttAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGaa                                                                                                        >  1:1035269/1‑140 (MQ=255)
                            ttAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGaa                                                                                                        >  2:1023954/1‑140 (MQ=255)
                            ttAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGaa                                                                                                        >  2:790121/1‑140 (MQ=255)
                               aCGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAAtt                                                                                                     >  1:287243/1‑140 (MQ=255)
                               aCGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAAtt                                                                                                     >  1:52202/1‑140 (MQ=255)
                                   aTTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATg                                                                                                 <  2:1461863/140‑1 (MQ=255)
                                      aaaCTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTc                                                                                              <  1:1323128/140‑1 (MQ=255)
                                       aaCTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCt                                                                                             <  2:1282433/140‑1 (MQ=255)
                                           ctACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTga                                                                                         >  2:1526955/1‑140 (MQ=255)
                                            tACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTgag                                                                                        >  1:881368/1‑140 (MQ=255)
                                            tACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTgag                                                                                        <  1:1512595/140‑1 (MQ=255)
                                            tACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTgag                                                                                        >  2:1568535/1‑140 (MQ=255)
                                            tACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTgag                                                                                        >  2:192831/1‑140 (MQ=255)
                                            tACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTgag                                                                                        >  2:852007/1‑140 (MQ=255)
                                              cTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGaa                                                                                      >  2:1221087/1‑140 (MQ=255)
                                                   tatGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGca                                                                                 <  2:1029300/140‑1 (MQ=255)
                                                     tGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGcaca                                                                               <  1:747804/140‑1 (MQ=255)
                                                     tGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGcaca                                                                               <  1:856030/140‑1 (MQ=255)
                                                     tGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGcaca                                                                               <  1:576515/140‑1 (MQ=255)
                                                       aTTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACAcc                                                                             <  2:375483/140‑1 (MQ=255)
                                                       aTTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACAcc                                                                             <  2:1052592/140‑1 (MQ=255)
                                                          tCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACcatt                                                                          <  1:729719/140‑1 (MQ=255)
                                                           cAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACcattc                                                                         >  2:666079/1‑140 (MQ=255)
                                                           cAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACcattc                                                                         >  2:144592/1‑140 (MQ=255)
                                                            aGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACcattca                                                                        <  2:1035017/140‑1 (MQ=255)
                                                            aGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACcattca                                                                        <  1:1276396/140‑1 (MQ=255)
                                                            aGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACcattca                                                                        <  2:795981/140‑1 (MQ=255)
                                                            aGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACcattca                                                                        <  1:805749/140‑1 (MQ=255)
                                                             gACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACcattcat                                                                       >  2:384093/1‑140 (MQ=255)
                                                              aCGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACcattcatt                                                                      >  2:978692/1‑140 (MQ=255)
                                                                              aTCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCa                                                      <  2:623422/140‑1 (MQ=255)
                                                                              aTCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCa                                                      <  1:1023954/140‑1 (MQ=255)
                                                                              aTCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCa                                                      <  2:17447/140‑1 (MQ=255)
                                                                              aTCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCa                                                      <  2:287243/140‑1 (MQ=255)
                                                                                    tttGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGtt                                                >  1:542407/1‑140 (MQ=255)
                                                                                     ttGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGttt                                               >  1:552757/1‑140 (MQ=255)
                                                                                      tGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTg                                              >  1:714048/1‑140 (MQ=255)
                                                                                        cAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTa                                            <  2:1070034/140‑1 (MQ=255)
                                                                                        cAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTa                                            <  1:934235/140‑1 (MQ=255)
                                                                                         aaaaGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTaa                                           <  2:563170/140‑1 (MQ=255)
                                                                                           aaGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATg                                         <  1:1526955/140‑1 (MQ=255)
                                                                                                      ccTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATa                              <  1:192831/140‑1 (MQ=255)
                                                                                                        tGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACg                            <  1:790121/140‑1 (MQ=255)
                                                                                                           aTTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGttt                         <  2:52202/140‑1 (MQ=255)
                                                                                                           aTTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGttt                         <  1:53579/140‑1 (MQ=255)
                                                                                                           aTTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGttt                         <  2:850065/140‑1 (MQ=255)
                                                                                                                tGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACTAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAGAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGc                    <  1:161165/140‑1 (MQ=255)
                                                                                                                tGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGc                    <  2:277435/140‑1 (MQ=255)
                                                                                                                tGCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGc                    >  2:868835/1‑140 (MQ=255)
                                                                                                                 gCGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGcc                   >  2:667730/1‑140 (MQ=255)
                                                                                                                  cGTCTGTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCt                  >  2:1320815/1‑140 (MQ=255)
                                                                                                                       gTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATa             <  2:714048/140‑1 (MQ=255)
                                                                                                                       gTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATa             <  1:1221087/140‑1 (MQ=255)
                                                                                                                       gTAACGGTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATa             <  2:880203/140‑1 (MQ=255)
                                                                                                                            ggTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCgtt        >  2:1111176/1‑140 (MQ=255)
                                                                                                                             gTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCgttg       >  1:447491/1‑140 (MQ=255)
                                                                                                                             gTCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCgttg       >  1:1038530/1‑140 (MQ=255)
                                                                                                                              tCATGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCgttgt      >  1:587004/1‑140 (MQ=255)
                                                                                                                                 tGTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGa   <  1:852007/140‑1 (MQ=255)
                                                                                                                                  gTAACTCCTGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGAt  >  2:832549/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                   
TAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTACCGTTACCTGAATTACTGCGTCTGTAACGGTCATGTAACTCCTGGGGGCCATCGATCGAGATACCAACCAGAAATTCATGTTCTTTGAGAAAGGCACACCATTCATTATTCAATAAAATGCCATTCGTTTGTAATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGAT  >  NC_000913/1580291‑1580561

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