Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,907,713 (T)7→8 intergenic (‑122/+548) yobF ← / ← yebO DUF2527 family heat‑induced protein/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,907,7061.T95.6% 232.6 / 2.8 68intergenic (‑115/+555)yobF/yebODUF2527 family heat‑induced protein/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base T (35/30);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

AAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGT  >  NC_000913/1907578‑1907822
                                                                                                                                 |                                                                                                                     
aaaaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACttttttttgc                                                                                                            >  1:565770/1‑136 (MQ=255)
 aaaTGCCACACCTAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACttttttttgca                                                                                                           <  2:395147/139‑5 (MQ=255)
 aaaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACttttttttgca                                                                                                           <  2:211272/139‑5 (MQ=255)
 aaaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACttttttttgca                                                                                                           <  2:610781/139‑5 (MQ=255)
 aaaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACttttttttgca                                                                                                           <  2:405759/139‑5 (MQ=255)
          acaTAAACAGAACTGTACCTCGTTTAAACCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGc                                                                                                   >  1:363651/1‑139 (MQ=255)
                 cAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGaa                                                                                            <  1:416728/139‑1 (MQ=255)
                    aaCTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATAc                                                                                         >  1:109086/1‑139 (MQ=255)
                     aCTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATAcg                                                                                        >  1:391434/1‑139 (MQ=255)
                      cTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATAcgc                                                                                       >  1:122792/1‑139 (MQ=255)
                      cTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATAcgc                                                                                       <  2:140515/139‑1 (MQ=255)
                        gttccTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCCCAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACgcgc                                                                                     <  1:403475/136‑1 (MQ=255)
                         tACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGcc                                                                                    >  1:160978/1‑139 (MQ=255)
                          aCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCg                                                                                   <  1:199320/139‑1 (MQ=255)
                            cTCGTTTATCCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGat                                                                                 <  2:363651/139‑1 (MQ=255)
                            cTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGat                                                                                 >  2:564604/1‑139 (MQ=255)
                             tCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGata                                                                                <  2:391434/139‑1 (MQ=255)
                                tttAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATAcgc                                                                                       <  1:15240/129‑1 (MQ=255)
                                tttAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATAcgc                                                                                       >  2:15240/1‑129 (MQ=255)
                                  tAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTGTGCACACTTACGCAACGGAATACGCGCCGATATATCa                                                                           >  1:388830/1‑139 (MQ=255)
                                        gAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGc                                                                     <  2:122792/139‑1 (MQ=255)
                                                       aTACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGtt                                                      <  2:109086/139‑1 (MQ=255)
                                                        tACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGc                                                                     >  2:265293/1‑123 (MQ=255)
                                                        tACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGc                                                                     <  1:265293/123‑1 (MQ=255)
                                                        tACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACttttttttg                                                                                                             <  2:141405/82‑3 (MQ=255)
                                                        tACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACttttttttg                                                                                                             >  1:141405/1‑80 (MQ=255)
                                                          cAACGTTAATTATGGCAAGCCATTTTTAGACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCAACTTTTTTTTGCACACTTATGCAACGGAATAcg                                                                                        >  2:378741/1‑102 (MQ=255)
                                                          cAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATAcg                                                                                        <  1:378741/102‑1 (MQ=255)
                                                          cAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACttttttttgca                                                                                                           <  2:586113/82‑5 (MQ=255)
                                                          cAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACttttttttgca                                                                                                           >  1:586113/1‑78 (MQ=255)
                                                             cGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCa                                                >  2:423780/1‑139 (MQ=255)
                                                             cGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCa                                                >  2:441122/1‑139 (MQ=255)
                                                             cGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCa                                                >  1:52326/1‑139 (MQ=255)
                                                             cGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCa                                                <  1:310421/139‑1 (MQ=255)
                                                               ttAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATAcg                                                                                        <  1:551631/97‑1 (MQ=255)
                                                               ttAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATAcg                                                                                        >  2:551631/1‑97 (MQ=255)
                                                                tAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGAt                                             >  1:303655/1‑139 (MQ=255)
                                                                 aaTTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACgcgc                                                                                     <  1:345119/98‑1 (MQ=255)
                                                                 aaTTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACgcgc                                                                                     >  2:345119/1‑98 (MQ=255)
                                                                 aaTTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATa                                            >  1:300915/1‑139 (MQ=255)
                                                                       ggCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAg                                               >  1:232430/1‑130 (MQ=255)
                                                                       ggCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAg                                               <  2:232430/130‑1 (MQ=255)
                                                                       ggCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTaaa                                      >  1:246441/1‑139 (MQ=255)
                                                                       ggCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTaaa                                      >  1:307309/1‑139 (MQ=255)
                                                                       ggCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTaaa                                      >  1:52192/1‑139 (MQ=255)
                                                                         cAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACt                                    >  2:194677/1‑139 (MQ=255)
                                                                           aGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTgg                                  >  1:73003/1‑139 (MQ=255)
                                                                            gCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGc                                 >  1:336762/1‑139 (MQ=255)
                                                                            gCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGc                                 >  1:318385/1‑139 (MQ=255)
                                                                             ccATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCg                                <  2:565770/139‑1 (MQ=255)
                                                                               aTTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTa                              <  2:73003/139‑1 (MQ=255)
                                                                                        cATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAAc                     <  2:246441/139‑1 (MQ=255)
                                                                                        cATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAAc                     <  2:300915/139‑1 (MQ=255)
                                                                                          tGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATAcgc                                                                                       >  2:321312/1‑71 (MQ=255)
                                                                                          tGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATAcgc                                                                                       <  1:321312/71‑1 (MQ=255)
                                                                                           gTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACgcgc                                                                                     >  1:472824/1‑72 (MQ=255)
                                                                                           gTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACgcgc                                                                                     <  2:472824/72‑1 (MQ=255)
                                                                                           gTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCg                                <  2:499958/125‑1 (MQ=255)
                                                                                           gTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCg                                >  1:499958/1‑125 (MQ=255)
                                                                                           gTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATa                  >  1:260980/1‑139 (MQ=255)
                                                                                           gTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACAt                   >  1:482753/1‑137 (MQ=255)
                                                                                           gTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACAt                   <  2:482753/137‑1 (MQ=255)
                                                                                                     gACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTa                              >  2:601961/1‑117 (MQ=255)
                                                                                                     gACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTa                              <  1:601961/117‑1 (MQ=255)
                                                                                                     gACGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAg        >  1:441557/1‑139 (MQ=255)
                                                                                                      aCGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCAt                                                                          >  2:89380/1‑72 (MQ=255)
                                                                                                      aCGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCAt                                                                          <  1:89380/72‑1 (MQ=255)
                                                                                                      aCGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAAc                                                                <  2:113730/82‑1 (MQ=255)
                                                                                                      aCGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAAc                                                                >  1:113730/1‑82 (MQ=255)
                                                                                                      aCGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAAt           >  2:282406/1‑135 (MQ=255)
                                                                                                      aCGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAAt           <  1:282406/135‑1 (MQ=255)
                                                                                                      aCGCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGt       >  1:154790/1‑139 (MQ=255)
                                                                                                        gCACAAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTGAATAAGTTc     >  2:390338/1‑139 (MQ=255)
                                                                                                          acaAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAg   <  2:52192/139‑1 (MQ=255)
                                                                                                           caAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAgt  <  2:441557/139‑1 (MQ=255)
                                                                                                           caAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAgt  <  2:318385/139‑1 (MQ=255)
                                                                                                           caAATCTCGATACGTCTTCCACTTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAgt  <  1:441122/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                     
AAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGT  >  NC_000913/1907578‑1907822

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