Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,300,930 T→C E608E (GAA→GAG napA ← nitrate reductase, periplasmic, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,300,9300TC97.1% 339.1 / ‑0.5 102E608E (GAA→GAGnapAnitrate reductase, periplasmic, large subunit
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (44/55);  total (45/57)
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.42e-01

Only 100 of 103 total aligned reads displayed.

CGTGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGTTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCCAGCAGATCTTCCGGC  >  NC_000913/2300833‑2301029
                                                                                                 |                                                                                                   
cgtgaccgCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTc                                                                                                   <  1:2649236/100‑1 (MQ=255)
 gtgaccgCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCg                                                                                                  <  2:1451340/100‑1 (MQ=255)
  tgaccgCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCgg                                                                                                 <  2:919672/100‑1 (MQ=255)
  tgaccgCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCgg                                                                                                 <  2:1562436/100‑1 (MQ=255)
   gaccgCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGa                                                                                                <  2:2251046/100‑1 (MQ=255)
    accgCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGAt                                                                                               <  1:470571/100‑1 (MQ=255)
     ccgCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATa                                                                                              >  1:865934/1‑100 (MQ=255)
        cgACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCg                                                                                           >  1:1108324/1‑100 (MQ=255)
         gACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCgg                                                                                          <  2:678975/100‑1 (MQ=255)
         gACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCgg                                                                                          <  2:501427/100‑1 (MQ=255)
          aCCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCggg                                                                                         <  1:988498/100‑1 (MQ=255)
           ccAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGa                                                                                        >  1:1869126/1‑100 (MQ=255)
           ccAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGa                                                                                        <  2:2277869/100‑1 (MQ=255)
           ccAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGa                                                                                        <  2:2201265/100‑1 (MQ=255)
             aaaCCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAAt                                                                                      >  2:1731859/1‑100 (MQ=255)
             aaaCCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAAt                                                                                      >  1:2104487/1‑100 (MQ=255)
             aaaCCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAAt                                                                                      >  2:900768/1‑100 (MQ=255)
              aaCCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGaa                                                                                       <  2:2409817/98‑1 (MQ=255)
               aCCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAAttt                                                                                    <  2:918563/100‑1 (MQ=255)
                 cATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGc                                                                                  <  1:564468/100‑1 (MQ=255)
                  aTGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCt                                                                                 <  1:933003/100‑1 (MQ=255)
                   tGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTc                                                                                <  1:50917/100‑1 (MQ=255)
                    gCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCa                                                                               <  1:1975177/100‑1 (MQ=255)
                     cGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCAc                                                                              <  1:660585/100‑1 (MQ=255)
                          tctTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCg                                                                         <  1:1214932/100‑1 (MQ=255)
                          tctTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCg                                                                         <  1:1002138/100‑1 (MQ=255)
                           ctTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCgg                                                                        <  2:2637891/100‑1 (MQ=255)
                               gAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGgtgt                                                                     >  2:1465084/1‑99 (MQ=255)
                               gAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGtgtg                                                                    <  1:2727354/100‑1 (MQ=255)
                                aaCAGCCCTTTTTGCGGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTgg                                                                   <  1:2842962/100‑1 (MQ=255)
                                aaCAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTgg                                                                   <  1:2352653/100‑1 (MQ=255)
                                 aCAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGc                                                                  <  2:922447/100‑1 (MQ=255)
                                 aCAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGc                                                                  <  1:1302521/100‑1 (MQ=255)
                                   aGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCAt                                                                >  1:139997/1‑100 (MQ=255)
                                    gCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATa                                                               >  2:1430582/1‑100 (MQ=255)
                                    gCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATa                                                               >  1:27304/1‑100 (MQ=255)
                                      ccTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGTTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACa                                                             >  1:1182728/1‑100 (MQ=255)
                                      ccTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACa                                                             >  1:2086846/1‑100 (MQ=255)
                                       cTTTTTGCAGATAGAAGCCCAGCTCGCGGGCTTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAg                                                            >  1:1254566/1‑100 (MQ=255)
                                       cTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAg                                                            <  1:428624/100‑1 (MQ=255)
                                       cTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAg                                                            >  2:2667757/1‑100 (MQ=255)
                                        tttttGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCACACTTCCGGTGTGGCATACAGc                                                           <  2:508549/100‑1 (MQ=255)
                                        tttttGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCCGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGc                                                           <  2:1849655/100‑1 (MQ=255)
                                          tttGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCAc                                                         <  2:1781067/100‑1 (MQ=255)
                                            tGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACtt                                                       >  1:895460/1‑100 (MQ=255)
                                             gCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGTTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTc                                                      >  1:1234137/1‑100 (MQ=255)
                                             gCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTc                                                      >  1:420830/1‑100 (MQ=255)
                                             gCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTc                                                      >  2:1888658/1‑100 (MQ=255)
                                              cAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACTGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCg                                                     >  2:565615/1‑100 (MQ=255)
                                              cAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCg                                                     >  1:98175/1‑100 (MQ=255)
                                              cAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCg                                                     <  1:2283158/100‑1 (MQ=255)
                                              cAGATAGAAGCCCAGCTCGCGGGATTCATCGTTAAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCg                                                     >  1:630115/1‑100 (MQ=255)
                                               aGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGt                                                    <  1:1430588/100‑1 (MQ=255)
                                               aGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGt                                                    >  2:1111434/1‑100 (MQ=255)
                                                gATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTa                                                   <  2:2052515/100‑1 (MQ=255)
                                                 aTAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTAc                                                  <  1:587079/100‑1 (MQ=255)
                                                 aTAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTAc                                                  <  1:1112951/100‑1 (MQ=255)
                                                  tAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCAGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACa                                                 >  2:1739644/1‑100 (MQ=255)
                                                   aGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAg                                                <  1:996051/100‑1 (MQ=255)
                                                    gAAGCCCCGATCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGc                                               >  2:2700035/1‑100 (MQ=255)
                                                    gAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGc                                               >  2:2781751/1‑100 (MQ=255)
                                                       gCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGtt                                            >  1:2548756/1‑100 (MQ=255)
                                                       gCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGtt                                            >  1:1264633/1‑100 (MQ=255)
                                                         ccAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGtttt                                          <  1:2667777/100‑1 (MQ=255)
                                                         ccAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGtttt                                          >  2:1876224/1‑100 (MQ=255)
                                                          cAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAgatcggaa                                         >  1:962916/1‑93 (MQ=255)
                                                           aGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGc                                        <  1:1731872/100‑1 (MQ=255)
                                                           aGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGc                                        <  1:2117737/100‑1 (MQ=255)
                                                           aGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGc                                        <  2:490793/100‑1 (MQ=255)
                                                             cTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCa                                      <  1:37560/100‑1 (MQ=255)
                                                               cgcgGGATTCATCGTTCAGCTGATCTTCCGCCAGTTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCAc                                     <  2:2626891/99‑1 (MQ=255)
                                                                      ttCATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTcc                             <  2:408263/100‑1 (MQ=255)
                                                                        cATCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCgg                           >  2:516885/1‑100 (MQ=255)
                                                                         aTCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGTCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGt                          <  2:1588528/100‑1 (MQ=255)
                                                                         aTCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGt                          <  2:118984/100‑1 (MQ=255)
                                                                          tCGTTCAGCTGATCTTCCGCCAGTTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGtt                         <  2:1234133/100‑1 (MQ=255)
                                                                          tCGTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGtt                         <  2:1148035/100‑1 (MQ=255)
                                                                            gTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTttc                       >  2:603080/1‑100 (MQ=255)
                                                                            gTTCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTttc                       >  1:1712426/1‑100 (MQ=255)
                                                                             ttCAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTttct                      <  2:14300/100‑1 (MQ=255)
                                                                               cAGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCTCTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTttcttn                    <  2:406018/100‑2 (MQ=255)
                                                                                aGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCg                   <  2:2623672/100‑1 (MQ=255)
                                                                                aGCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCg                   <  2:887290/100‑1 (MQ=255)
                                                                                 gCTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGc                  <  1:1555031/100‑1 (MQ=255)
                                                                                  cTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGcc                 >  1:1964774/1‑100 (MQ=255)
                                                                                  cTGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGcc                 >  1:1142805/1‑100 (MQ=255)
                                                                                   tGATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCca                >  1:2496680/1‑100 (MQ=255)
                                                                                    gATCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGcc                 >  2:2167661/1‑98 (MQ=255)
                                                                                     aTCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCcagc              <  2:1108320/100‑1 (MQ=255)
                                                                                     aTCTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCcagc              <  1:2423578/100‑1 (MQ=255)
                                                                                       cTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCcagcag            >  1:1601301/1‑100 (MQ=255)
                                                                                       cTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCcagcag            <  2:420829/100‑1 (MQ=255)
                                                                                       cTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCcagcag            >  2:1024982/1‑100 (MQ=255)
                                                                                       cTTCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCcagcag            >  1:379445/1‑100 (MQ=255)
                                                                                         tCCGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCCAGCAGAt          >  1:1837416/1‑100 (MQ=255)
                                                                                          ccGCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCCAGCAGATc         >  2:2152736/1‑100 (MQ=255)
                                                                                            gCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCCAGCAGATCtt       >  1:2590447/1‑100 (MQ=255)
                                                                                            gCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCCAGCAGATCtt       >  1:2852883/1‑100 (MQ=255)
                                                                                            gCCAGCTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCCAGCAGATCtt       >  2:2785603/1‑100 (MQ=255)
                                                                                                 ttCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCCAGCAGATCTTCCGGc  <  1:2801074/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
CGTGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGTTCAGCTGATCTTCCGCCAGTTCGGATACCGGGAATTTGCTCACTTCCGGTGTGGCATACAGCACTTCGTACAGCGTTTTGCCACGCAGTTCCGGTTTCTTCGCCAGCAGATCTTCCGGC  >  NC_000913/2300833‑2301029

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