Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,202,047 T→C 6.4% M590T (ATG→ACG)  yehI → DUF4132 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,202,0470TC 6.4% 188.6 / 2.8 62M590T (ATG→ACG) yehIDUF4132 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (26/32);  new base C (2/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TCCCCGGAAGCACAATTAGCAGAACAGGCGCTCCGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTA  >  NC_000913/2201899‑2202195
                                                                                                                                                    |                                                                                                                                                    
tCCCCGGAAGCACAATTAGCAGAACAGGCGCTCCGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGAt                                                                                                                                                      >  8:121294/1‑149 (MQ=255)
   ccGGAAGCACAATTAGCAGAACAGGCGCTCCGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGc                                                                                                                                                   >  3:316646/1‑149 (MQ=255)
   ccGGAAGCACAATTAGCAGAACAGGCGCTCCGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGc                                                                                                                                                   <  5:284750/149‑1 (MQ=255)
        aGCACAATTAGCAGAACAGGCGCTCCGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGa                                                                                                                                              <  3:315399/149‑1 (MQ=255)
        aGCACAATTAGCAGAACAGGCGCTCCGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGa                                                                                                                                              <  8:229692/149‑1 (MQ=255)
        aGCACAATTAGCAGAACAGGCGCTCCGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGa                                                                                                                                              <  7:195759/149‑1 (MQ=255)
        aGCACAATTAGCAGAACAGGCGCTCCGAAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGa                                                                                                                                              <  1:22332/149‑1 (MQ=255)
        aGCACAATTAGCAGAACAGGCCCTCCGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGAGAACTACGCACATTATCGTCCGGAAAAGAGTGTCAGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGa                                                                                                                                              >  5:146521/1‑149 (MQ=255)
                   cAGAACAGGCGCTCCGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTg                                                                                                                                   >  7:197049/1‑149 (MQ=255)
                     gAACAGGCGCTCCGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGATGCCGCGTGAGATGGCTTTGACCTTGTGGCTGCg                                                                                                                                 >  1:88850/1‑149 (MQ=255)
                      aaCAGGCGCTCCGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACACACATTATCGTCCCGAAAAGAGTGTCGGGAACCAGTATTTATTGGCACAGATGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCg                                                                                                                                 >  1:121375/1‑148 (MQ=255)
                       aCAGGCGCTCCGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTa                                                                                                                               >  8:244037/1‑149 (MQ=255)
                            cgcTCCGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAAt                                                                                                                          <  8:197049/149‑1 (MQ=255)
                                 cGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTGCGGAAAAGAGTGGCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATgag                                                                                                                       <  6:150636/147‑1 (MQ=255)
                                    aGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGaaaaa                                                                                                                  <  4:307925/149‑1 (MQ=255)
                                         tCACATAGAAGCGTTTATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCa                                                                                                               <  6:38940/147‑1 (MQ=255)
                                           acaTAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCt                                                                                                            <  4:131376/148‑1 (MQ=255)
                                                         attaaCATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAgg                                                                                              <  6:38366/148‑1 (MQ=255)
                                                          ttaaCATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGAc                                                                                            <  4:302297/149‑1 (MQ=255)
                                                           taaCATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGAc                                                                                            >  6:296240/1‑148 (MQ=255)
                                                             aCATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTAt                                                                                         <  1:45888/149‑1 (MQ=255)
                                                              cATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTAtt                                                                                        <  3:55669/149‑1 (MQ=255)
                                                               atatGTGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTAttt                                                                                       >  3:268616/1‑149 (MQ=255)
                                                                tatGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTAtttt                                                                                      <  4:316646/149‑1 (MQ=255)
                                                                     gAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAg                                                                                 >  2:108783/1‑149 (MQ=255)
                                                                              aaCTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCg                                                                         >  8:262359/1‑148 (MQ=255)
                                                                               actactCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGa                                                                       <  7:262359/149‑1 (MQ=255)
                                                                               actactCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGa                                                                       >  7:334587/1‑149 (MQ=255)
                                                                                  actCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACCGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTg                                                                    >  4:39966/1‑149 (MQ=255)
                                                                                    tCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGACGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGa                                                                  >  8:213590/1‑149 (MQ=255)
                                                                                      acaTTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTGTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGa                                                                  <  8:210280/147‑1 (MQ=255)
                                                                                        aTTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATgcg                                                              >  1:214503/1‑149 (MQ=255)
                                                                                        aTTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATgcg                                                              >  6:228619/1‑149 (MQ=255)
                                                                                              gTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGACGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTAc                                                          >  8:126676/1‑147 (MQ=255)
                                                                                                   gAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTAtt                                                                                        <  2:225775/112‑1 (MQ=255)
                                                                                                   gAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTAtt                                                                                        >  1:225775/1‑112 (MQ=255)
                                                                                                        gagTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGt                                              >  4:206558/1‑149 (MQ=255)
                                                                                                         agTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGtt                                             <  4:268616/149‑1 (MQ=255)
                                                                                                         agTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGtt                                             <  5:228619/149‑1 (MQ=255)
                                                                                                          gtgtCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTg                                            <  4:249930/149‑1 (MQ=255)
                                                                                                           tgtCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTgg                                           <  7:121294/149‑1 (MQ=255)
                                                                                                           ggtCTGGAACCTGTATTTATTGGCACAGTTGCCGCGGGAGATGGCTTTGACCTTCTGGCTGCGTAGCACTGAGAAAAAGCATCTGTTCGCGGGGGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTgg                                           <  2:88850/148‑1 (MQ=255)
                                                                                                             tCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGACGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGc                                          <  7:126676/148‑1 (MQ=255)
                                                                                                                ggAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACATTCTGGCTGCGTATCAATGAGAAATAGCATCTGTTCGCGGGTGAGGACTATTTTCTGAGTTTGCTCGGATTGGATGCGCTAGCAGGTCTGCTGTTGGCtttt                                      >  8:160440/1‑149 (MQ=255)
                                                                                                                 gAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCttttt                                     <  5:296240/149‑1 (MQ=255)
                                                                                                                  aaCCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTc                                    <  7:244037/149‑1 (MQ=255)
                                                                                                                     cTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTcaca                                 >  3:87534/1‑149 (MQ=255)
                                                                                                                     cTGTATTTATTGGCACAGGTGCCGCGTGAGATGGCGTTGACCTTCTGGCTGCGTATCAATGAGAAACAGCATCTGTGCGCGGGTGAGGACTAATTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTcaca                                 <  5:139894/149‑1 (MQ=255)
                                                                                                                      tGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGc                                                                                                                                  >  8:271352/1‑51 (MQ=255)
                                                                                                                      tGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGc                                                                                                                                  <  7:271352/51‑1 (MQ=255)
                                                                                                                        tatttattGGCACAGTTGCCGCGGGAGATGACTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATc                               <  2:105041/148‑1 (MQ=255)
                                                                                                                          tttattGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTc                            >  7:130185/1‑149 (MQ=255)
                                                                                                                            tattGGCACAGTTGCCGCGTGAGACGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTcc                           <  7:213590/148‑1 (MQ=255)
                                                                                                                                ggCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCTTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCaaaa                       <  1:108783/148‑1 (MQ=255)
                                                                                                                                 gCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCttt                                       <  2:11177/131‑1 (MQ=255)
                                                                                                                                 gCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCttt                                       >  1:11177/1‑131 (MQ=255)
                                                                                                                                 gCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAg                      >  7:307992/1‑148 (MQ=255)
                                                                                                                                 gCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAg                      >  3:1120/1‑148 (MQ=255)
                                                                                                                                 gCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGCATCAACGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAg                      >  1:65394/1‑148 (MQ=255)
                                                                                                                                     aGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCCTTTTCACATCGTCCAAAAGAAACa                 >  3:168050/1‑149 (MQ=255)
                                                                                                                                          ccGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTcc            <  1:134881/149‑1 (MQ=255)
                                                                                                                                          ccGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTc             <  3:170122/148‑1 (MQ=255)
                                                                                                                                              gTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGtt         >  6:62944/1‑148 (MQ=255)
                                                                                                                                                    tGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTa  <  4:87534/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCCCCGGAAGCACAATTAGCAGAACAGGCGCTCCGCAGCGGTCACATAGAAGCGTTAATTAACATATGGGAAAATGACAACTACTCACATTATCGTCCGGAAAAGAGTGTCTGGAACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTA  >  NC_000913/2201899‑2202195

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