Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,302,455 G→T 6.0% P100Q (CCG→CAG)  napA ← nitrate reductase, periplasmic, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,302,4550GT 6.0% 285.6 / 7.0 83P100Q (CCG→CAG) napAnitrate reductase, periplasmic, large subunit
Reads supporting (aligned to +/- strand):  ref base G (40/38);  new base T (2/3);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCA  >  NC_000913/2302308‑2302597
                                                                                                                                                   |                                                                                                                                              
ccAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCgg                                                                                                                                               <  4:219460/149‑1 (MQ=255)
 cAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGc                                                                                                                                              <  3:465018/149‑1 (MQ=255)
  aGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCt                                                                                                                                             >  3:485909/1‑149 (MQ=255)
   gAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTg                                                                                                                                            <  2:412011/149‑1 (MQ=255)
    aaCCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGc                                                                                                                                           <  3:358611/149‑1 (MQ=255)
       cGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTTATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTc                                                                                                                                        <  1:530180/149‑1 (MQ=255)
       cGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCTGCTGCGta                                                                                                                                        <  1:194039/149‑2 (MQ=255)
       cGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCTGCTGCGta                                                                                                                                        <  2:407851/149‑2 (MQ=255)
          aCATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCaaa                                                                                                                                     >  4:21792/1‑149 (MQ=255)
          aCATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCaaa                                                                                                                                     <  4:304694/149‑1 (MQ=255)
                    cgatTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCACACGGTCTTTAc                                                                                                                           >  2:6053/1‑149 (MQ=255)
                     gatTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGGCTTTAcc                                                                                                                          >  1:366175/1‑149 (MQ=255)
                     gatTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTACTCTTCCATTACATCGATGGCCTGATCGCAAGTGCTTGGGGTAAATTCGCCTTCTCTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTAcc                                                                                                                          >  1:627642/1‑149 (MQ=255)
                      atTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCg                                                                                                                         <  2:445150/149‑1 (MQ=255)
                      atTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCg                                                                                                                         >  2:197634/1‑149 (MQ=255)
                        tCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTa                                                                                                                       <  4:485909/149‑1 (MQ=255)
                            gCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATg                                                                                                                   >  4:273197/1‑149 (MQ=255)
                               cTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGAtt                                                                                                                >  1:561775/1‑149 (MQ=255)
                               cTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGAtt                                                                                                                <  2:561775/149‑1 (MQ=255)
                                tttttCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTa                                                                                                                       >  4:236868/1‑141 (MQ=255)
                                tttttCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTa                                                                                                                       <  3:236868/141‑1 (MQ=255)
                                   ttctttcAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTgg                                                                                                            >  2:656454/1‑149 (MQ=255)
                                    tctttcAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTggg                                                                                                           <  2:415825/149‑1 (MQ=255)
                                      tttcAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCt                                                                                                                               <  3:556290/127‑1 (MQ=255)
                                      tttcAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCt                                                                                                                               >  4:556290/1‑127 (MQ=255)
                                      tttcAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCa                                                                                                         >  2:301218/1‑149 (MQ=255)
                                      tttcAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCa                                                                                                         >  4:354853/1‑149 (MQ=255)
                                             gCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATa                                                                                                  <  3:354853/149‑1 (MQ=255)
                                              cGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAg                                                                                                 <  1:656454/149‑1 (MQ=255)
                                                gTTTTGTATTTCTCTTCCATCACATCGAAGGCCTGATACCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTAc                                                                                                                           >  2:235572/1‑121 (MQ=255)
                                                gTTTTGAATTTCTCTTCCATCTCATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTAc                                                                                                                           <  1:235572/121‑1 (MQ=255)
                                                 ttttGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGccc                                                                                              >  4:101530/1‑149 (MQ=255)
                                                 ttttGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGccc                                                                                              >  3:324580/1‑149 (MQ=255)
                                                 ttttGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGccc                                                                                              >  1:27620/1‑149 (MQ=255)
                                                  tttGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCt                                                                                             <  2:27620/149‑1 (MQ=255)
                                                        tttCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGc                                                                                       >  3:206608/1‑149 (MQ=255)
                                                              ttCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCa                                                                                 >  4:60049/1‑149 (MQ=255)
                                                              ttCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCa                                                                                 <  3:21792/149‑1 (MQ=255)
                                                              ttCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAAGAGCGCTTGATGGAGTTCa                                                                                 <  2:627642/149‑1 (MQ=255)
                                                               tCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAg                                                                                >  4:181397/1‑149 (MQ=255)
                                                                 cATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGc                                                                              <  4:582236/149‑1 (MQ=255)
                                                                   tCACATCGAAGGCCTGATCGCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATCTTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCa                                                                            <  1:197634/149‑1 (MQ=255)
                                                                      caTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACgg                                                                         <  2:677061/149‑1 (MQ=255)
                                                                                cTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCTGCTGCGta                                                                                                                                        >  2:683031/1‑75 (MQ=255)
                                                                                cTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCTGCTGCGta                                                                                                                                        <  1:683031/76‑2 (MQ=255)
                                                                                cTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTg                                                               >  3:10887/1‑149 (MQ=255)
                                                                                cTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTg                                                               >  4:567573/1‑149 (MQ=255)
                                                                                cTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTg                                                               >  2:313669/1‑149 (MQ=255)
                                                                                 tGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGc                                                              <  4:558848/149‑1 (MQ=255)
                                                                                   aTCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGt                                                            >  3:44910/1‑149 (MQ=255)
                                                                                      ccAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCg                                                         >  1:495451/1‑149 (MQ=255)
                                                                                       cAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCgg                                                        <  3:60049/149‑1 (MQ=255)
                                                                                        aGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCggg                                                       >  1:185323/1‑149 (MQ=255)
                                                                                          gTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTc                                                     >  3:212838/1‑149 (MQ=255)
                                                                                           tGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCg                                                    <  1:301218/149‑1 (MQ=255)
                                                                                            gATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGc                                                                                       >  2:207590/1‑113 (MQ=255)
                                                                                            gATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGc                                                                                       <  1:207590/113‑1 (MQ=255)
                                                                                             aTTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCAc                                                                           >  3:321715/1‑124 (MQ=255)
                                                                                             aTTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCAc                                                                           <  4:321715/124‑1 (MQ=255)
                                                                                                 gggTAGATTCGCCTTCTTTGTCATATTTACCGTATTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGACATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGa                                              >  1:196969/1‑149 (MQ=255)
                                                                                                 gggTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGa                                              <  2:185323/149‑1 (MQ=255)
                                                                                                    tAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAg                                           <  2:152762/149‑1 (MQ=255)
                                                                                                       aTTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGcc                                        >  1:344100/1‑149 (MQ=255)
                                                                                                        ttCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGtt                                                                       >  3:623614/1‑117 (MQ=255)
                                                                                                        ttCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGtt                                                                       <  4:623614/117‑1 (MQ=255)
                                                                                                         tCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGccac                                      <  4:10887/149‑1 (MQ=255)
                                                                                                         tCGCCTTCTTAGACGTATTTGCCGTATTTCATGCGCGGCAGCGACTGCGTCAAACGGGCTTTTCCGTACACGATTATCGGCGGGAAATAGCCCTTATTGCAGTTCAGGCCACGCTTAACCGGTGCGTCCGGGTCGCCCTGACGGGccac                                      <  2:344100/149‑1 (MQ=255)
                                                                                                           gCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGccacca                                    >  3:490696/1‑149 (MQ=255)
                                                                                                              ttcttTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCacac                                 <  2:620559/149‑1 (MQ=255)
                                                                                                                     tCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCctg                          >  1:609820/1‑149 (MQ=255)
                                                                                                                     gcATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGCCAGGCCACCACACGTCCctg                          >  1:674109/2‑149 (MQ=255)
                                                                                                                      cATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCctgc                         <  3:567573/149‑1 (MQ=255)
                                                                                                                      cATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCctgc                         <  1:313669/149‑1 (MQ=255)
                                                                                                                       atatTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCctgct                        <  2:674109/149‑1 (MQ=255)
                                                                                                                                 tttttCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGAc              >  2:186004/1‑149 (MQ=255)
                                                                                                                                 tttttCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGAc              <  2:366175/149‑1 (MQ=255)
                                                                                                                                   tttCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCa                                                                                      <  1:64331/75‑1 (MQ=255)
                                                                                                                                   tttCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCa                                                                                      >  2:64331/1‑75 (MQ=255)
                                                                                                                                    gtttgaCGCAGCCGCTGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAg           >  4:632479/6‑149 (MQ=255)
                                                                                                                                     tCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGa          <  3:101530/149‑1 (MQ=255)
                                                                                                                                      cATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGaa         <  2:495451/149‑1 (MQ=255)
                                                                                                                                       aTACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAc        >  4:341456/1‑149 (MQ=255)
                                                                                                                                         aCGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcgc      >  2:349016/1‑149 (MQ=255)
                                                                                                                                             agcagcGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCa  >  3:599015/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCA  >  NC_000913/2302308‑2302597

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