Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,618,075 C→T *234* (TAG→TAA hda ← ATPase regulatory factor involved in DnaA inactivation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,618,0750CT97.5% 244.8 / ‑3.3 79*234* (TAG→TAAhdaATPase regulatory factor involved in DnaA inactivation
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (38/39);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TGTGCTCCATCGCTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACGCC  >  NC_000913/2617940‑2618208
                                                                                                                                       |                                                                                                                                     
tgtgCTCCATCGCTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTAc                                                                                                                                    <  1:680029/139‑1 (MQ=255)
       cATCGCTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCa                                                                                                                             >  1:739586/1‑139 (MQ=255)
               gTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAAtttct                                                                                                                     >  1:1906096/1‑139 (MQ=255)
                          ttGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACgg                                                                                                          >  1:834469/1‑139 (MQ=255)
                           tGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGa                                                                                                         >  2:308357/1‑139 (MQ=255)
                              cGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                     <  2:538819/139‑1 (MQ=255)
                               ggATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                     <  2:139715/139‑1 (MQ=255)
                                     gACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCt                                                                                               >  1:1206891/1‑139 (MQ=255)
                                                      aTCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATc                                                                              >  2:317926/1‑139 (MQ=255)
                                                           gCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGc                                                                         >  1:980122/1‑139 (MQ=255)
                                                           gCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGc                                                                         >  1:83713/1‑139 (MQ=255)
                                                            ccTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                        <  2:1206891/139‑1 (MQ=255)
                                                            ccTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                        <  2:673091/139‑1 (MQ=255)
                                                            ccTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                        <  2:834469/139‑1 (MQ=255)
                                                                cGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGAt                                                                    <  2:980122/139‑1 (MQ=255)
                                                                cGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGAt                                                                    <  2:1281502/139‑1 (MQ=255)
                                                                          ggATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAAc                                                                                                            <  1:919972/89‑1 (MQ=255)
                                                                          ggATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAAc                                                                                                            >  2:919972/1‑89 (MQ=255)
                                                                          ggATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATc                                                          >  1:921177/1‑139 (MQ=255)
                                                                          ggATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATc                                                          >  2:1714460/1‑139 (MQ=255)
                                                                            aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                        <  1:308357/139‑1 (MQ=255)
                                                                              tAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAAc                                                      <  2:83713/139‑1 (MQ=255)
                                                                               aaCTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                     >  2:398842/1‑91 (MQ=255)
                                                                               aaCTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                     <  1:398842/91‑1 (MQ=255)
                                                                                aCTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTc                                                                                                   <  2:1927930/92‑1 (MQ=255)
                                                                                aCTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTc                                                                                                   >  1:1927930/1‑92 (MQ=255)
                                                                                aCTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGt                                                    <  2:1906096/139‑1 (MQ=255)
                                                                                aCTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGt                                                    >  2:1820137/1‑139 (MQ=255)
                                                                                 cTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTc                                                   >  1:1376447/1‑139 (MQ=255)
                                                                                   tCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                        <  1:2013395/132‑1 (MQ=255)
                                                                                   tCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                        >  2:2013395/1‑132 (MQ=255)
                                                                                    cgcgTAGTTCGGATAAGGCGTTCGCGCCGCATTCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTGAGCTTTCGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                        <  2:1374627/131‑1 (MQ=255)
                                                                                    cgcgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                        >  1:1374627/1‑131 (MQ=255)
                                                                                     gcgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTgcg                                                                                      <  1:1513779/100‑1 (MQ=255)
                                                                                     gcgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTgcg                                                                                      >  2:1513779/1‑100 (MQ=255)
                                                                                      cgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                     <  1:1147942/84‑1 (MQ=255)
                                                                                      cgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                     >  2:1147942/1‑84 (MQ=255)
                                                                                          gTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                        <  1:1307129/109‑1 (MQ=255)
                                                                                          gTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                        >  2:1307129/1‑109 (MQ=255)
                                                                                           ttCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                        >  1:1136689/1‑108 (MQ=255)
                                                                                           ttCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                        <  2:1136689/108‑1 (MQ=255)
                                                                                               gATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTgcg                                     <  2:739586/139‑1 (MQ=255)
                                                                                               gATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACCCACGATCCAACTGATCCAACGTCATAAATAGCGTgcg                                     <  2:1376447/139‑1 (MQ=255)
                                                                                                aTAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                        >  2:624439/1‑119 (MQ=255)
                                                                                                aTAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                        <  1:624439/119‑1 (MQ=255)
                                                                                                aTAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGt                                                    <  2:86752/123‑1 (MQ=255)
                                                                                                aTAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGt                                                    >  1:86752/1‑123 (MQ=255)
                                                                                                    ggCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGt                                                                                  >  2:1412143/1‑89 (MQ=255)
                                                                                                    ggCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGt                                                                                  <  1:1412143/89‑1 (MQ=255)
                                                                                                    ggCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCAttt                                >  1:1107626/1‑139 (MQ=255)
                                                                                                     gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                     <  1:477244/69‑1 (MQ=255)
                                                                                                     gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                     >  2:477244/1‑69 (MQ=255)
                                                                                                     gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTc                                                                                                   <  1:1192205/71‑1 (MQ=255)
                                                                                                     gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTc                                                                                                   >  2:1192205/1‑71 (MQ=255)
                                                                                                     gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                        <  2:392645/114‑1 (MQ=255)
                                                                                                     gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                        >  1:392645/1‑114 (MQ=255)
                                                                                                     gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTtc                               >  1:1464118/1‑139 (MQ=255)
                                                                                                          cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                        <  2:1394791/93‑1 (MQ=255)
                                                                                                          cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                        >  1:1394791/1‑93 (MQ=255)
                                                                                                          cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTgcg                                     >  2:108827/1‑128 (MQ=255)
                                                                                                          cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTgcg                                     <  1:108827/128‑1 (MQ=255)
                                                                                                          cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGt                          >  2:1878019/1‑139 (MQ=255)
                                                                                                               cgcATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGc                     >  1:1721746/1‑139 (MQ=255)
                                                                                                               cgcATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGc                     >  2:1961360/1‑139 (MQ=255)
                                                                                                                  aTCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGc                  <  1:1820137/139‑1 (MQ=255)
                                                                                                                   tCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCt                 <  2:1721746/139‑1 (MQ=255)
                                                                                                                   tCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCt                 <  2:1464118/139‑1 (MQ=255)
                                                                                                                             aaCACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGaa       <  2:1107626/139‑1 (MQ=255)
                                                                                                                             aaCACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGaa       <  1:1714460/139‑1 (MQ=255)
                                                                                                                             aaCACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGaa       <  1:1961360/139‑1 (MQ=255)
                                                                                                                              acacCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACg                                                                          >  1:7458/1‑71 (MQ=255)
                                                                                                                              acacCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACg                                                                          <  2:7458/71‑1 (MQ=255)
                                                                                                                              acacCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGa                       >  1:2000237/1‑122 (MQ=255)
                                                                                                                              acacCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGa                       <  2:2000237/122‑1 (MQ=255)
                                                                                                                                  cTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                        >  2:1031633/1‑85 (MQ=255)
                                                                                                                                  cTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                        <  1:1031633/85‑1 (MQ=255)
                                                                                                                                  cTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACGcc  <  1:317926/139‑1 (MQ=255)
                                                                                                                                   ttATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAAc                                                               <  1:922324/77‑1 (MQ=255)
                                                                                                                                   ttATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAAc                                                               >  2:922324/1‑77 (MQ=255)
                                                                                                                                       |                                                                                                                                     
TGTGCTCCATCGCTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACGCC  >  NC_000913/2617940‑2618208

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