Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,291,901 C→A F130L (TTC→TTA nrfE → putative cytochrome c‑type biogenesis protein NrfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,291,9010CA98.7% 273.5 / ‑2.9 78F130L (TTC→TTAnrfEputative cytochrome c‑type biogenesis protein NrfE
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (32/45);  total (33/45)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

CTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGCGGCGAGCGTGGCGCTGGCGA  >  NC_000913/4291764‑4292037
                                                                                                                                         |                                                                                                                                        
ctgctTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGt                                                                                                                                        >  2:1245362/1‑140 (MQ=255)
 tgctTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTg                                                                                                                                       <  1:66417/140‑1 (MQ=255)
 tgctTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTg                                                                                                                                       <  1:155955/140‑1 (MQ=255)
    tttCCGCCTGGAGCGCGCTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCgg                                                                                                                                    <  2:832342/140‑1 (MQ=255)
        cGCCTGGAGCGCGCTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGatat                                                                                                                                >  1:126357/1‑140 (MQ=255)
        cGCCTGGAGCGCGCTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGatat                                                                                                                                >  1:300686/1‑140 (MQ=255)
             ggAGCGCGCTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTcca                                                                                                                           <  2:108543/140‑1 (MQ=255)
             ggAGCGCGCTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTcca                                                                                                                           >  1:61907/1‑140 (MQ=255)
              gAGCGCGCTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTctac                                                                                                                          <  2:1450869/140‑1 (MQ=255)
                 cgcgcTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAg                                                                                                                       <  1:1184235/140‑1 (MQ=255)
                 cgcgcTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAg                                                                                                                       <  1:862948/140‑1 (MQ=255)
                   cgcTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCa                                                                                                                     >  1:431087/1‑140 (MQ=255)
                     cTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAAt                                                                                                                   >  2:1147992/1‑140 (MQ=255)
                        tttGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGa                                                                                                                <  2:1391456/140‑1 (MQ=255)
                              tGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCcg                                                                                                          <  2:1278548/140‑1 (MQ=255)
                               ggCATTATCGGCATCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCcgc                                                                                                         <  2:758657/140‑1 (MQ=255)
                               ggCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCcgc                                                                                                         <  1:1542957/140‑1 (MQ=255)
                                gCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCcgcg                                                                                                        >  1:1596599/1‑140 (MQ=255)
                                gCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCcgcg                                                                                                        >  2:1041471/1‑140 (MQ=255)
                                  aTTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGAt                                                                                                      >  1:562315/1‑140 (MQ=255)
                                  aTTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGAt                                                                                                      <  1:1631488/140‑1 (MQ=255)
                                    tATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGAtct                                                                                                    <  1:1041471/140‑1 (MQ=255)
                                         gcagcaAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTACCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATc                                                                                               <  2:1752486/140‑1 (MQ=255)
                                           agcaAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCg                                                                                             >  1:1786848/1‑140 (MQ=255)
                                            gcaAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGa                                                                                            <  1:566964/140‑1 (MQ=255)
                                              aaaCCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGAtg                                                                                          <  2:742540/140‑1 (MQ=255)
                                              aaaCCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGAtg                                                                                          <  1:1147992/140‑1 (MQ=255)
                                                aCCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGAtgct                                                                                        >  1:248298/1‑140 (MQ=255)
                                                 ccGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGAtgctg                                                                                       <  2:562315/140‑1 (MQ=255)
                                                   gATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCa                                                                                     <  2:1763726/140‑1 (MQ=255)
                                                   gATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCa                                                                                     <  2:300686/140‑1 (MQ=255)
                                                    aTCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCa                                                                                     >  2:1100903/1‑139 (MQ=255)
                                                    aTCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCa                                                                                     <  1:1100903/139‑1 (MQ=255)
                                                        gCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATc                                                                                <  1:1272927/140‑1 (MQ=255)
                                                        gCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATc                                                                                <  2:431087/140‑1 (MQ=255)
                                                           aTTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCg                                                                             <  1:1245362/140‑1 (MQ=255)
                                                               ccGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCt                                                                         <  2:61907/140‑1 (MQ=255)
                                                                cGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCtt                                                                        <  2:1628766/140‑1 (MQ=255)
                                                                   tGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATc                                                                     >  1:231404/1‑140 (MQ=255)
                                                                    gACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATcc                                                                               <  1:1521328/129‑1 (MQ=255)
                                                                    gACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAAAATcc                                                                               >  1:148789/1‑129 (MQ=255)
                                                                    gACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAAAATcc                                                                               <  2:148789/129‑1 (MQ=255)
                                                                    gACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGATACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATcc                                                                               >  2:1521328/1‑129 (MQ=255)
                                                                        ttAGCCGTTTTATCTCTCATGCTCACCGCATTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGGTGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAGTCCGATGGTGCAACATCCCGGTCTTATCTTTCa                                                                <  2:1534190/140‑1 (MQ=255)
                                                                              gTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCAcc                                                          <  2:231404/140‑1 (MQ=255)
                                                                                tttATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCa                                                                >  2:1773807/1‑132 (MQ=255)
                                                                                tttATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCa                                                                <  1:1773807/132‑1 (MQ=255)
                                                                                    tctctcATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCgctgct                                                    <  2:320197/140‑1 (MQ=255)
                                                                                    tctctcATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCgctgct                                                    <  2:887281/140‑1 (MQ=255)
                                                                                       ctcATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTa                                                 <  1:1815806/140‑1 (MQ=255)
                                                                                         cATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTAcc                                               <  2:1704698/140‑1 (MQ=255)
                                                                                                  cgcACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATg                                      <  2:126357/140‑1 (MQ=255)
                                                                                                  cgcACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATg                                      >  2:1205061/1‑140 (MQ=255)
                                                                                                  cgcACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATg                                      >  2:671623/1‑140 (MQ=255)
                                                                                                  cgcACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATg                                      >  2:1772127/1‑140 (MQ=255)
                                                                                                   gcACTGCTACTGTTTGTGGTGCTGCGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATgg                                     <  2:873278/140‑1 (MQ=255)
                                                                                                    cACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGc                                    <  1:50296/140‑1 (MQ=255)
                                                                                                                   tggtgCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTTATGGTAgcg                     >  2:479321/1‑140 (MQ=255)
                                                                                                                   tggtgCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAgcg                     >  2:194473/1‑140 (MQ=255)
                                                                                                                      tgCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAgcggcg                  >  1:972956/1‑140 (MQ=255)
                                                                                                                      tgCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAgcggcg                  >  1:117047/1‑140 (MQ=255)
                                                                                                                      tgCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAgcggcg                  >  1:962921/1‑140 (MQ=255)
                                                                                                                       gCTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAgcggcgg                 >  2:1397442/1‑139 (MQ=255)
                                                                                                                        cTGTGGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGCGGCGAg                >  1:844838/1‑140 (MQ=255)
                                                                                                                         tgtgGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGCGGCGGGc               >  2:1189164/1‑140 (MQ=255)
                                                                                                                          gtgGTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGCGGCGAGCg              <  2:1786848/140‑1 (MQ=255)
                                                                                                                           tgGTCCGATCCCTTCGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTTATGGGAGCGGCGAGCgg             >  1:722488/1‑139 (MQ=255)
                                                                                                                            ggTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGCGGCGAGCGTg            <  2:952784/140‑1 (MQ=255)
                                                                                                                             gTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGCGGCGAGCGTgg           <  2:920157/140‑1 (MQ=255)
                                                                                                                             gTCCGATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGCGGCGAGCGTgg           >  1:626559/1‑140 (MQ=255)
                                                                                                                                 gATCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGCGGCGAGCGTGGCGCt       <  2:972956/140‑1 (MQ=255)
                                                                                                                                  aTCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGTTGGTAGCGGCGGGCGTGGGGCTg      >  2:510871/1‑140 (MQ=255)
                                                                                                                                  aTCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAgcgg                    >  1:949308/1‑126 (MQ=255)
                                                                                                                                  aTCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAgcgg                    <  2:949308/126‑1 (MQ=255)
                                                                                                                                  aTCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGCGGCGAGCGGGGGGCTg      >  2:1062366/1‑140 (MQ=255)
                                                                                                                                   tCCCTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGCGGCGAGCGTGGCGCTgg     <  1:1856669/140‑1 (MQ=255)
                                                                                                                                     ccTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGCGGCGAGCGTGGCGCTGGCg   >  1:1159229/1‑140 (MQ=255)
                                                                                                                                      cTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGCGGCGAGCGTGGCGCTGGCGa  <  2:117047/140‑1 (MQ=255)
                                                                                                                                      cTTAGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGCGGCGAGCGTGGCGCTGGCGa  <  2:962921/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
CTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTATCGGCAGCAAACCGATCCGCTATTTCCGCTGACGTTAGCCGTTTTATCTCTCATGCTCGCCGCACTGCTACTGTTTGTGGTGCTGTGGTCCGATCCCTTCGTGCGGATATTTCCACCAGCAATCGAAGGCCGCGATCTCAATCCGATGCTGCAACATCCCGGTCTTATCTTTCATCCACCGCTGCTTTACCTTGGCTATGGCGGTTTGATGGTAGCGGCGAGCGTGGCGCTGGCGA  >  NC_000913/4291764‑4292037

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: