Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 399,801 A→G 5.3% Q70R (CAA→CGA)  yaiZ → DUF2754 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913399,8010AG 5.3% 262.0 / 3.5 76Q70R (CAA→CGA) yaiZDUF2754 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (38/34);  new base G (2/2);  total (40/36)
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.87e-01

TTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTG  >  NC_000913/399656‑399949
                                                                                                                                                 |                                                                                                                                                    
ttcattcTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATa                                                                                                                                                   >  4:1464833/1‑149 (MQ=255)
 tcattcTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATaa                                                                                                                                                  <  2:389412/149‑1 (MQ=255)
           tGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCGATAAGAATATTATc                                                                                                                                        <  2:246492/149‑1 (MQ=255)
             gCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCGATAATATTCTTATTgg                                                                                                                                      >  2:548170/1‑149 (MQ=255)
               gtggtgGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCGATAATATTCTTATTGGGc                                                                                                                                    <  1:548170/149‑1 (MQ=255)
                      ggTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTaaaa                                                                                                                             <  3:1464833/149‑1 (MQ=255)
                        ttACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAata                                                                                                                           >  4:285036/1‑149 (MQ=255)
                        ttACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAata                                                                                                                           <  3:285036/149‑1 (MQ=255)
                           cTGGGATTTGATGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTa                                                                                                                        >  3:945821/1‑149 (MQ=255)
                                     aaGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGtt                                                                                                              >  3:426507/1‑149 (MQ=255)
                                       gCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGtttt                                                                                                            >  3:855017/1‑149 (MQ=255)
                                        cAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAata                                                                                                                                             <  2:730301/115‑1 (MQ=255)
                                        cAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAata                                                                                                                                             >  1:730301/1‑115 (MQ=255)
                                            aGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCTGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATCGTGGTTTTCAAtc                                                                                                       <  4:945821/149‑2 (MQ=255)
                                            aGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAatat                                                                                                                          >  4:1104305/1‑130 (MQ=255)
                                            aGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAatat                                                                                                                          <  3:1104305/130‑1 (MQ=255)
                                            aGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAatat                                                                                                                          >  2:368552/1‑130 (MQ=255)
                                            aGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAatat                                                                                                                          <  1:368552/130‑1 (MQ=255)
                                              gTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATtgtg                                                                                                                 >  2:721498/1‑137 (MQ=255)
                                              gTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATtgtg                                                                                                                 <  1:721498/137‑1 (MQ=255)
                                                    aGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATcc                                                                                               <  4:855017/149‑1 (MQ=255)
                                                        cTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAatat                                                                                                                                            <  3:745308/100‑1 (MQ=255)
                                                        cTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAatat                                                                                                                                            >  4:745308/1‑100 (MQ=255)
                                                        cTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAg                                                                                           >  3:856791/1‑149 (MQ=255)
                                                              cGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGAc                                                                                     >  1:790673/1‑149 (MQ=255)
                                                              cGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGAc                                                                                     >  2:1064880/1‑149 (MQ=255)
                                                               gTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAatat                                                                                                                          <  1:1416755/111‑1 (MQ=255)
                                                               gTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAatat                                                                                                                          >  2:1416755/1‑111 (MQ=255)
                                                                tCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGc                                                                                   <  4:856791/149‑1 (MQ=255)
                                                                   gTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTcc                                                                                >  4:1380077/1‑149 (MQ=255)
                                                                   gTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTcc                                                                                <  4:426507/149‑1 (MQ=255)
                                                                   gTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTcc                                                                                >  2:674224/1‑149 (MQ=255)
                                                                        gTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAgcgc                                                                           >  3:337926/1‑149 (MQ=255)
                                                                              tGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCg                                                                                                                                       >  2:144913/1‑83 (MQ=255)
                                                                              tGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCg                                                                                                                                       <  1:144913/83‑1 (MQ=255)
                                                                              tGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAata                                                                                                                           <  3:937492/95‑1 (MQ=255)
                                                                              tGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAata                                                                                                                           >  4:937492/1‑95 (MQ=255)
                                                                               gggTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGtt                                                                      >  3:1051923/1‑147 (MQ=255)
                                                                               gggTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGtt                                                                      <  4:1051923/147‑1 (MQ=255)
                                                                                ggTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTtcaa                                                                   <  1:1064880/149‑1 (MQ=255)
                                                                                    aTTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCACAATAAGAATATTATCGGGCGTTAAAATATTACATtgt                                                                                                                  <  3:919740/98‑1 (MQ=255)
                                                                                    aTTAGTTGCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACAttct                                                                                                                  >  4:919740/1‑96 (MQ=255)
                                                                                     ttAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCaaa                                                              <  4:337926/149‑1 (MQ=255)
                                                                                       aGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATc                                                                                                >  1:296434/1‑113 (MQ=255)
                                                                                       aGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATc                                                                                                <  2:296434/113‑1 (MQ=255)
                                                                                        gTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACgt                                                           <  3:1380077/149‑1 (MQ=255)
                                                                                         tttCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCa                                                                                                          >  3:1114578/1‑101 (MQ=255)
                                                                                         tttCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCa                                                                                                          <  4:1114578/101‑1 (MQ=255)
                                                                                         tttCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACgtg                                                          >  4:873229/1‑149 (MQ=255)
                                                                                              ggctggcGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCaa                                                                                                         >  4:726656/1‑97 (MQ=255)
                                                                                              ggctggcGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCaa                                                                                                         <  3:726656/97‑1 (MQ=255)
                                                                                                  ggcGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTCATCAGAt                                                 <  2:790673/149‑1 (MQ=255)
                                                                                                         ttGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGgt                                             <  1:566416/146‑1 (MQ=255)
                                                                                                         ttGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGgt                                             >  2:566416/1‑146 (MQ=255)
                                                                                                               attCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAg                                    >  1:854620/1‑149 (MQ=255)
                                                                                                                ttCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTGACCCAGa                                   >  1:1193873/1‑149 (MQ=255)
                                                                                                                 tcgtcgtcgCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGt                                                                                                    >  3:838039/1‑83 (MQ=255)
                                                                                                                 tcgtcgtcgCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGt                                                                                                    <  4:838039/83‑1 (MQ=255)
                                                                                                                    tcgtcgCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGc                               >  3:962937/1‑149 (MQ=255)
                                                                                                                          cgATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCtt                         <  4:962937/149‑1 (MQ=255)
                                                                                                                            aTGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGc                       <  3:873229/149‑1 (MQ=255)
                                                                                                                               aagaaACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAgcgc                                                                           >  4:120938/1‑94 (MQ=255)
                                                                                                                               aagaaACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAgcgc                                                                           <  3:120938/94‑1 (MQ=255)
                                                                                                                                agaaACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTtcaat                                                                  >  1:396484/1‑102 (MQ=255)
                                                                                                                                agaaACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTtcaat                                                                  <  2:396484/102‑1 (MQ=255)
                                                                                                                                agaaACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCcaca                   >  4:241396/1‑149 (MQ=255)
                                                                                                                                 gaaACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGAc                                                                                     >  2:997552/1‑82 (MQ=255)
                                                                                                                                 gaaACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGAc                                                                                     <  1:997552/82‑1 (MQ=255)
                                                                                                                                 gaaACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAg                  >  3:739720/1‑149 (MQ=255)
                                                                                                                                 gaaACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAg                  >  4:83088/1‑149 (MQ=255)
                                                                                                                                        aaaaCGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGa           <  4:739720/149‑1 (MQ=255)
                                                                                                                                        aaaaCGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGa           <  1:674224/149‑1 (MQ=255)
                                                                                                                                              cccGATAATATTCTTATTGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAg                                                                        >  2:1389740/1‑82 (MQ=37)
                                                                                                                                               ccAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGatacatac    >  1:715790/1‑149 (MQ=255)
                                                                                                                                                cAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTaa                                         <  1:155563/111‑1 (MQ=255)
                                                                                                                                                cAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTCTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTaa                                         >  2:155563/1‑111 (MQ=255)
                                                                                                                                                 aaTAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCca                     <  2:1034119/130‑1 (MQ=255)
                                                                                                                                                 aaTAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCca                     >  1:1034119/1‑130 (MQ=255)
                                                                                                                                                 aaTAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCGCAGCTTCGGatacatac    >  4:405986/1‑147 (MQ=255)
                                                                                                                                                 aaTAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCGCAGCTTCGGatacatac    <  3:405986/147‑1 (MQ=255)
                                                                                                                                                 aaTAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTg  >  1:620396/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTCATTCTTAATGGCGTGGTGGGGTTACTGGGATTTGAAGCAAAAGGTTGGCAGACCTATGCCGTCGGTCTGGTGACGTGGGTGATTAGTTTCTGGCTGGCGGGGTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTATCCGCAGCGTGACGCTCCAGCGCCAGTTCAATCAAACGTGTGATCAGATCGGTGTAACCCAGACCGCTGGCTTGCCACAGCTTCGGATACATACTG  >  NC_000913/399656‑399949

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