Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,302,455 G→T 8.9% 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 8.9% 311.3 / 16.5 90P100Q (CCG→CAG) napAnitrate reductase, periplasmic, large subunit
Reads supporting (aligned to +/- strand):  ref base G (41/41);  new base T (4/4);  total (45/45)
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.31e-01

CCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCC  >  NC_000913/2302308‑2302594
                                                                                                                                                   |                                                                                                                                           
ccAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCgg                                                                                                                                            <  4:1366385/149‑1 (MQ=255)
 cAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCTGc                                                                                                                                           >  2:760397/1‑149 (MQ=255)
 cAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGc                                                                                                                                           <  4:248076/149‑1 (MQ=255)
     aCCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCt                                                                                                                                          >  3:722542/1‑146 (MQ=255)
     aCCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCt                                                                                                                                          <  4:722542/146‑1 (MQ=255)
       cGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCTGCTGCGta                                                                                                                                     <  1:760397/149‑2 (MQ=255)
           cATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAAc                                                                                                                                 <  1:881787/149‑1 (MQ=255)
                                tttttCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGc                                                                                                                                        <  4:1418947/121‑1 (MQ=255)
                                tttttCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGGGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGc                                                                                                                                        >  3:1418947/1‑121 (MQ=255)
                                  tttctttcAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTg                                                                                                          <  4:562435/149‑1 (MQ=255)
                                      tttcAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGc                                                                                                                                        >  2:94423/1‑115 (MQ=255)
                                      tttcAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGc                                                                                                                                        <  1:94423/115‑1 (MQ=255)
                                                 ttttGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAAc                                                                                                                                 <  1:280901/111‑1 (MQ=255)
                                                 ttttGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAAc                                                                                                                                 >  2:280901/1‑111 (MQ=255)
                                                   ttGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCtt                                                                                         >  2:1236507/1‑149 (MQ=255)
                                                      aaTTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCa                                                                                                      <  3:570663/133‑1 (MQ=255)
                                                      aaTTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCa                                                                                                      >  4:570663/1‑133 (MQ=255)
                                                       aTTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTgg                                                                                                         <  4:1013162/129‑1 (MQ=255)
                                                       aTTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTgg                                                                                                         >  3:1013162/1‑129 (MQ=255)
                                                        tttCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTg                                                                                                          >  3:342735/1‑127 (MQ=255)
                                                        tttCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTg                                                                                                          <  4:342735/127‑1 (MQ=255)
                                                        tttCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGc                                                                                    >  1:21184/1‑149 (MQ=255)
                                                         ttCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCa                                                                                   >  4:1332169/1‑149 (MQ=255)
                                                              ttCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCTGCTGCGtatg                                                                                                                                   >  1:963362/1‑93 (MQ=255)
                                                              ttCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCTGCTGCGtatg                                                                                                                                   <  2:963362/96‑4 (MQ=255)
                                                              ttCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCa                                                                              >  4:687936/1‑149 (MQ=255)
                                                                   tCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGt                                                                                                                              >  1:458148/1‑96 (MQ=255)
                                                                   tCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGt                                                                                                                              <  2:458148/96‑1 (MQ=255)
                                                                   tCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGa                                                                                                   >  2:16504/1‑123 (MQ=255)
                                                                   tCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGa                                                                                                   <  1:16504/123‑1 (MQ=255)
                                                                   tCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCa                                                                         >  3:659122/1‑149 (MQ=255)
                                                                    cacaTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCa                                                                         >  2:604454/1‑148 (MQ=255)
                                                                    cacaTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCa                                                                         <  1:604454/148‑1 (MQ=255)
                                                                    cacaTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCAc                                                                        >  1:1274632/1‑149 (MQ=255)
                                                                    cacaTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCAc                                                                        <  2:21184/149‑1 (MQ=255)
                                                                     acaTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACg                                                                       >  2:940397/1‑149 (MQ=255)
                                                                      caTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTc                                                                                                                             <  3:696232/94‑1 (MQ=255)
                                                                      caTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTc                                                                                                                             >  4:696232/1‑94 (MQ=255)
                                                                       aTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGtt                                                                                >  1:995285/1‑138 (MQ=255)
                                                                       aTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGtt                                                                                <  2:995285/138‑1 (MQ=255)
                                                                          gAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTaa                                                                                       >  3:773126/1‑128 (MQ=255)
                                                                          gAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTaa                                                                                       <  4:773126/128‑1 (MQ=255)
                                                                               ccTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGt                                                             <  4:659122/149‑1 (MQ=255)
                                                                               ccTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGt                                                             <  2:1274632/149‑1 (MQ=255)
                                                                                cTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGcc                                                                                            >  3:625819/1‑117 (MQ=255)
                                                                                cTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGcc                                                                                            <  4:625819/117‑1 (MQ=255)
                                                                                 tGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTaa                                                                                       <  3:1019809/121‑1 (MQ=255)
                                                                                 tGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTaa                                                                                       >  4:1019809/1‑121 (MQ=255)
                                                                                    tCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCAc                                                                        >  2:976924/1‑133 (MQ=255)
                                                                                    tCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCAc                                                                        <  1:976924/133‑1 (MQ=255)
                                                                                      ccAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGc                                                                           <  4:724927/128‑1 (MQ=255)
                                                                                      ccAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGc                                                                           >  3:724927/1‑128 (MQ=255)
                                                                                      ccAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCAc                                                                        >  3:1306440/1‑131 (MQ=255)
                                                                                      ccAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCAc                                                                        <  4:1306440/131‑1 (MQ=255)
                                                                                        aGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGcc                                                                                            >  2:564073/1‑109 (MQ=255)
                                                                                        aGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGcc                                                                                            <  1:564073/109‑1 (MQ=255)
                                                                                        aGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAAc                                                                 <  3:448185/136‑1 (MQ=255)
                                                                                        aGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCggg                                                    <  3:1216651/149‑1 (MQ=255)
                                                                                        aGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCggg                                                    >  4:1216651/1‑149 (MQ=255)
                                                                                        aGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGAGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAAc                                                                 >  4:448185/1‑136 (MQ=255)
                                                                                                ggggTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTg                                            >  3:332355/1‑149 (MQ=255)
                                                                                                 gggTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGAtttt                                                                                                           >  4:166730/1‑85 (MQ=255)
                                                                                                 gggTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGAtttt                                                                                                           <  3:166730/85‑1 (MQ=255)
                                                                                                   gTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACa                                         >  1:1356627/1‑149 (MQ=255)
                                                                                                       aTTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATa                                                                                               <  2:516537/91‑1 (MQ=255)
                                                                                                       aTTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATa                                                                                               >  1:516537/1‑91 (MQ=255)
                                                                                                          cGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGccacc                                  <  3:687936/149‑1 (MQ=255)
                                                                                                           gCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATg                                                                                     >  2:1044749/1‑97 (MQ=255)
                                                                                                           gCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATg                                                                                     <  1:1044749/97‑1 (MQ=255)
                                                                                                              ttcttTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCacac                              >  3:879677/1‑149 (MQ=255)
                                                                                                                 tttGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTc                           >  1:993358/1‑149 (MQ=255)
                                                                                                                  ttGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTcc                          <  1:940397/149‑1 (MQ=255)
                                                                                                                     tCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGt                                                             <  2:599891/111‑1 (MQ=255)
                                                                                                                     tCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGt                                                             >  1:599891/1‑111 (MQ=255)
                                                                                                                     tCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGccc                                              >  4:874690/1‑126 (MQ=255)
                                                                                                                     tCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGccc                                              <  3:874690/126‑1 (MQ=255)
                                                                                                                           ttACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGt                 <  2:1356627/149‑1 (MQ=255)
                                                                                                                             aCCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGcc                                               >  2:1318449/1‑117 (MQ=255)
                                                                                                                             aCCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGcc                                               <  1:1318449/117‑1 (MQ=255)
                                                                                                                                gTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGa                                           >  4:386962/1‑118 (MQ=255)
                                                                                                                                gTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGa                                           <  3:386962/118‑1 (MQ=255)
                                                                                                                                   tttCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCa         <  1:1236507/149‑1 (MQ=255)
                                                                                                                                    ttCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGcc                                               >  4:25306/1‑110 (MQ=255)
                                                                                                                                    ttCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGcc                                               <  3:25306/110‑1 (MQ=255)
                                                                                                                                     tttgaCGCAGCCGCTGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGa       >  1:283536/5‑149 (MQ=255)
                                                                                                                                       tgaCGCAGCCGCTGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGt                                                             >  2:1218033/3‑93 (MQ=255)
                                                                                                                                       tgaCGCAGCCGCTGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGt                                                             <  1:1218033/91‑1 (MQ=255)
                                                                                                                                        tACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcg    <  4:332355/149‑1 (MQ=255)
                                                                                                                                          cGCAGCCGCTGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcgcc  <  2:283536/149‑1 (MQ=255)
                                                                                                                                           gcagcaGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGc                                      <  4:1292503/112‑1 (MQ=255)
                                                                                                                                           gcagcaGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGc                                      >  3:1292503/1‑112 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
CCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCC  >  NC_000913/2302308‑2302594

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