Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 472,811 G→T 6.7% E72* (GAA→TAA)  glnK → nitrogen assimilation regulatory protein for GlnL, GlnE, and AmtB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913472,8110GT 6.7% 211.6 / 5.7 60E72* (GAA→TAA) glnKnitrogen assimilation regulatory protein for GlnL, GlnE, and AmtB
Reads supporting (aligned to +/- strand):  ref base G (31/25);  new base T (2/2);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

CATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAA  >  NC_000913/472663‑472959
                                                                                                                                                    |                                                                                                                                                    
cATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATg                                                                                                                                                      <  3:546410/149‑1 (MQ=255)
cATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATg                                                                                                                                                      <  1:396218/149‑1 (MQ=255)
 aTTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGa                                                                                                                                                     >  1:352731/1‑149 (MQ=255)
  ttGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGaa                                                                                                                                                    <  1:88823/149‑1 (MQ=255)
  ttGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGaa                                                                                                                                                    >  1:74803/1‑149 (MQ=255)
  ttGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGaa                                                                                                                                                    <  1:337300/149‑1 (MQ=255)
     gTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTg                                                                                                                                                 >  4:337306/1‑149 (MQ=255)
      tATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGa                                                                                                                                                <  4:392411/149‑1 (MQ=255)
       aTTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGAt                                                                                                                                               <  2:27215/149‑1 (MQ=255)
                      gTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTaa                                                                                                                                <  3:289957/149‑1 (MQ=255)
                          ccGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAggcg                                                                                                                            >  3:420401/1‑149 (MQ=255)
                            gAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAggcggc                                                                                                                          <  2:74803/149‑1 (MQ=255)
                            gAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAggcggc                                                                                                                          >  4:240781/1‑149 (MQ=255)
                                  aaaGGTTTCGGGCGTCAGAGAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTACTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTacac                                                                                                                    >  4:390335/1‑149 (MQ=255)
                                  aaaGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTacac                                                                                                                    >  3:538700/1‑149 (MQ=255)
                                  aaaGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTacac                                                                                                                    >  1:329212/1‑149 (MQ=255)
                                            ggCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAAtt                                                                                                          <  3:615222/149‑1 (MQ=255)
                                              cGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTgg                                                                                                        <  4:7075/149‑1 (MQ=255)
                                                   gAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACg                                                                                                   >  4:196145/1‑149 (MQ=255)
                                                    aaaGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAg                                                                                                                                                   <  4:479842/100‑1 (MQ=255)
                                                    aaaGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAg                                                                                                                                                   >  3:479842/1‑100 (MQ=255)
                                                              ccGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTtc                                                                                        <  3:144835/149‑1 (MQ=255)
                                                              ccGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTtc                                                                                        >  4:93689/1‑149 (MQ=255)
                                                              ccGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTtc                                                                                        <  3:613935/149‑1 (MQ=255)
                                                              ccGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTtc                                                                                        <  1:237871/149‑1 (MQ=255)
                                                                 aGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTtcgtc                                                                                     <  4:318614/149‑1 (MQ=255)
                                                                        ccGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATTAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAAt                                                                              >  4:13310/1‑149 (MQ=255)
                                                                        ccGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATTAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAAt                                                                              >  2:204846/1‑149 (MQ=255)
                                                                             gggCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATTAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCaa                                                                         <  3:13310/149‑1 (MQ=255)
                                                                             gggCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATTAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCaa                                                                         <  1:204846/149‑1 (MQ=255)
                                                                             gggCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCaa                                                                         >  1:588711/1‑149 (MQ=255)
                                                                                   aaTACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCg                                                                                                                  >  2:61359/1‑102 (MQ=255)
                                                                                   aaTACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCg                                                                                                                  <  1:61359/102‑1 (MQ=255)
                                                                                   aaTACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTc                                                                   <  1:20170/149‑1 (MQ=255)
                                                                                        aGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCattcg                                                              >  3:57180/1‑149 (MQ=255)
                                                                                          cGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCattcgta                                                            <  2:329212/149‑1 (MQ=255)
                                                                                             cAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAgg                                                                                                                              <  4:264772/80‑1 (MQ=255)
                                                                                             cAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAgg                                                                                                                              >  3:264772/1‑80 (MQ=255)
                                                                                             cAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCattcgtattc                                                         <  2:352731/149‑1 (MQ=255)
                                                                                                  tCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTAcc                                                    <  3:240781/149‑1 (MQ=255)
                                                                                                          aaaGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGc                                            <  1:491876/149‑1 (MQ=255)
                                                                                                              tAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCcgac                                        <  3:442502/149‑1 (MQ=255)
                                                                                                              tAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCcgac                                        >  4:23903/1‑149 (MQ=255)
                                                                                                                   aTTGATGTGGCGATTGCTGATGACCAACTCTATGAAGTGATCGATATTGTCGGTAAGGCGGCTTAAACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAgc                                   >  1:315814/1‑149 (MQ=255)
                                                                                                                         gtgGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCt                             >  4:372071/1‑149 (MQ=255)
                                                                                                                           ggCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCt                             >  3:585053/1‑147 (MQ=255)
                                                                                                                           ggCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCt                             <  4:585053/147‑1 (MQ=255)
                                                                                                                           ggCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGt                           >  1:85825/1‑149 (MQ=255)
                                                                                                                              gATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAAt                        >  4:75112/1‑149 (MQ=255)
                                                                                                                                  gCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATctct                    <  3:23903/149‑1 (MQ=255)
                                                                                                                                    tgatgaCCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTgg                  >  4:109778/1‑149 (MQ=255)
                                                                                                                                          ccAACTCGATGAAGTGATCGATATCGTCAGTCAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCACTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGcacaca            >  4:379767/1‑149 (MQ=255)
                                                                                                                                           cAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCattcgtattcg                                                        >  2:537032/1‑104 (MQ=255)
                                                                                                                                           cAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCattcgtactcg                                                        <  1:537032/104‑1 (MQ=255)
                                                                                                                                            aaCTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGc          >  3:197698/1‑149 (MQ=255)
                                                                                                                                            aaCTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGc          >  3:375265/1‑149 (MQ=255)
                                                                                                                                             aCTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCa         >  2:116729/1‑149 (MQ=255)
                                                                                                                                             aCTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTAAACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCa         >  4:146497/1‑149 (MQ=255)
                                                                                                                                                 gATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAg     >  1:11742/1‑149 (MQ=255)
                                                                                                                                                 gATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAg     >  3:389639/1‑149 (MQ=255)
                                                                                                                                                  aggAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAgg    >  2:426353/3‑149 (MQ=255)
                                                                                                                                                  aTGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAgg    >  2:341144/1‑149 (MQ=255)
                                                                                                                                                  aTGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAgg    <  3:337306/149‑1 (MQ=255)
                                                                                                                                                    gAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGaa  <  2:315814/149‑1 (MQ=255)
                                                                                                                                                    gAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTGCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGaa  <  4:538700/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAA  >  NC_000913/472663‑472959

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