Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,335,655 T→A 16.6% intergenic (‑65/+39) eptA ← / ← adiC lipid A phosphoethanolamine transferase/arginine:agmatine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,335,6550TA16.6% 443.5 / 70.6 173intergenic (‑65/+39)eptA/adiClipid A phosphoethanolamine transferase/arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (56/88);  new base A (14/15);  total (70/103)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 174 total aligned reads displayed.

TCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAT  >  NC_000913/4335557‑4335752
                                                                                                  |                                                                                                    
 cAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTc                                                                                                     >  2:1100946/1‑99 (MQ=255)
    agagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCAAACAAAGTGCGCATATGCTCGAAAAGCAAAATTAAGCCAACCTTAAGTTCTTaa                                                                                                 >  2:665323/1‑100 (MQ=255)
       gggTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGt                                                                                              >  2:1739319/1‑100 (MQ=255)
        ggTCTTTTTAGTAGGCGCTTCAACCTCACGGTGTTTCCATCGAACACAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGtt                                                                                             <  1:2508958/100‑1 (MQ=255)
         gTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTg                                                                                            <  2:2853499/100‑1 (MQ=255)
          tCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCGTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTgg                                                                                           >  2:2730371/1‑100 (MQ=255)
          tCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTgg                                                                                           <  2:1750018/100‑1 (MQ=255)
           cTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGc                                                                                          >  2:1240289/1‑100 (MQ=255)
           cTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTg                                                                                            <  2:3439428/98‑1 (MQ=255)
           cTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTGGc                                                                                          <  2:1932287/100‑1 (MQ=255)
            tttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCt                                                                                         <  2:3639181/100‑1 (MQ=255)
            tttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCt                                                                                         >  2:1027769/1‑100 (MQ=255)
            tttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTGGCt                                                                                         <  1:673531/100‑1 (MQ=255)
            tttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTGGCt                                                                                         <  2:2600546/100‑1 (MQ=255)
             ttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCt                                                                                         >  1:2180959/1‑99 (MQ=255)
             ttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTCGCt                                                                                         >  1:2760188/1‑99 (MQ=255)
             ttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTGGCtt                                                                                        <  2:3257507/100‑1 (MQ=255)
             gcttatggggCGCTTCAACATCACGGTGGTTCCATCGAACAAAGTGCGCAGATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTgg                                                                                           <  1:3371539/89‑1 (MQ=37)
              ttcggtacGCGTTTCAATATCACGGTGTGTCCATACAACAAAGTTCGTATATTTTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTGGCtta                                                                                       <  1:3789233/96‑2 (MQ=16)
               ttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCtt‑‑tt                                                                                    <  2:2840070/100‑1 (MQ=255)
                tagtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTa                                                                                   <  2:2992959/100‑1 (MQ=255)
                  gtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATg                                                                                 >  2:1874820/1‑100 (MQ=255)
                  gtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATg                                                                                 >  2:3297458/1‑100 (MQ=255)
                    agGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑tt                                                                              >  1:3051267/1‑100 (MQ=255)
                    agGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTGGCttaattttg                                                                                 <  2:1796773/100‑8 (MQ=37)
                      gcgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTg                                                                            <  2:1428504/100‑1 (MQ=255)
                       cgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGc                                                                           <  1:3407938/100‑1 (MQ=255)
                       cgcTTCAACATCACGGTGTTTCCATAGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑ttt                                                                             >  1:2985452/1‑98 (MQ=255)
                            cAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTgg                                                                        >  2:229837/1‑98 (MQ=255)
                            cAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGAt                                                                      <  2:373701/100‑1 (MQ=255)
                             aaCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGAtt                                                                     <  1:3677751/100‑1 (MQ=255)
                             aaCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGAtt                                                                     >  2:454832/1‑100 (MQ=255)
                              aCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTa                                                                    <  1:3111282/100‑1 (MQ=255)
                              aCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAAACTTAAGGTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTa                                                                    >  1:2721759/1‑100 (MQ=255)
                                  cACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAggg                                                                <  1:90544/100‑1 (MQ=255)
                                      gtgtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAAt                                                            <  1:3661879/100‑1 (MQ=255)
                                       tgtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAAtt                                                           <  2:705393/100‑1 (MQ=255)
                                        gtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTa                                                          >  2:192304/1‑100 (MQ=255)
                                        gtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTa                                                          <  1:1547844/100‑1 (MQ=255)
                                           tCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATc                                                       <  2:3142149/100‑1 (MQ=255)
                                             cATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCtt                                                     >  1:3301364/1‑100 (MQ=255)
                                              aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCttt                                                    >  2:1511473/1‑100 (MQ=255)
                                               tCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTg                                                   <  2:3220600/100‑1 (MQ=255)
                                                 gAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCt                                                 <  2:2193130/100‑1 (MQ=255)
                                                  aaCAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCtt                                                >  2:2648193/1‑100 (MQ=255)
                                                     aaaGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGAAATTAATCTTTGTTTAtt                                             >  1:1040004/1‑100 (MQ=255)
                                                      aaGTGCGCATATGCTCGCAAAGCGAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                            <  1:505314/100‑1 (MQ=255)
                                                      aaGTGCGCATATGCTCGCAAAGCAGAATTAAGCCAACCTTAAGTTCTTGAGGTTGGCGT‑‑TGATG‑TTTGCTGGGTTAAGGGTAATTAATCTTTGCTTATTg                                            >  1:768079/1‑100 (MQ=255)
                                                       aGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTgg                                           <  2:1720374/100‑1 (MQ=255)
                                                         tGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTg                                         <  1:1716397/100‑1 (MQ=255)
                                                         tGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTg                                         <  2:2721780/100‑1 (MQ=255)
                                                          gcgcATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGc                                        <  2:593600/100‑1 (MQ=255)
                                                          gcgcATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGc                                        <  1:3433332/100‑1 (MQ=255)
                                                          gcgcATATGCTCGCAAAGCAAAATTAAGACAACCTCAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGATGGATAAAGAGGAATTAATCTTTGCTTATTg                                            >  1:3559499/1‑96 (MQ=255)
                                                             cATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATc                                     <  1:3741516/100‑1 (MQ=255)
                                                              atatGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAATGTTAGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCt                                    <  1:1079382/100‑1 (MQ=255)
                                                               tatGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTa                                   <  1:1257440/100‑1 (MQ=255)
                                                                 ttcTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAg                                 <  1:2609130/98‑1 (MQ=255)
                                                                 tGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAg                                 <  2:1323088/100‑1 (MQ=255)
                                                                 tGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAg                                 >  2:2061839/1‑100 (MQ=255)
                                                                  gCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAgg                                >  1:1316314/1‑100 (MQ=255)
                                                                   cTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAggg                               <  1:200803/100‑1 (MQ=255)
                                                                    tCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGa                              >  2:1373155/1‑100 (MQ=255)
                                                                     cGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGAt                             >  2:3315164/1‑100 (MQ=255)
                                                                      gtaaagcaaaATTAAGCCAACCTTAAGTACTTAAGGTTGACTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTAGTGCAGGTAAGGGa                              <  1:1678486/96‑1 (MQ=255)
                                                                      gcaaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCGTATTGGTGCATCTAAGGGATa                            >  2:2572823/1‑100 (MQ=255)
                                                                        aaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACg                          <  2:385363/100‑1 (MQ=255)
                                                                         aagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACgg                         >  2:1224087/1‑100 (MQ=255)
                                                                         aagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACgg                         >  1:3515217/1‑100 (MQ=255)
                                                                            caaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGtt                      >  2:3678985/1‑100 (MQ=255)
                                                                              cataaaaGCCAACCTTAAGAACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGtttt                    >  1:3316449/6‑100 (MQ=38)
                                                                               ataaaaGCCAACCTTAAGAACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGttttt                   >  2:936910/5‑100 (MQ=38)
                                                                               aaTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGttttt                   <  1:1551271/100‑1 (MQ=255)
                                                                                aTTAAGCCAGCCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATGGGTGCATCTAAGGGATACGGGTTTTTa                  <  1:1099159/100‑1 (MQ=255)
                                                                                aTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATATGGGTTTTTa                  <  1:3870473/100‑1 (MQ=255)
                                                                                aTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                  <  1:2572805/100‑1 (MQ=255)
                                                                                 ttAAGCCAACCTTAAGTTCTTAAGTTTGTCTT‑‑TTATT‑TTTGCTGGATTAAGGAGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                  <  1:2746585/99‑1 (MQ=255)
                                                                                 ttAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTGTG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                 <  2:2783865/100‑1 (MQ=255)
                                                                                 ttAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                 <  2:1048552/100‑1 (MQ=255)
                                                                                  tAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATg                <  2:2312631/100‑1 (MQ=255)
                                                                                  tAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATg                <  2:2791778/100‑1 (MQ=255)
                                                                                   aaGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGTGGAATTTATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                 >  1:2771015/1‑98 (MQ=255)
                                                                                    aGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTa              <  1:988117/100‑1 (MQ=255)
                                                                                     gCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAg             >  1:2947033/1‑100 (MQ=255)
                                                                                     gCCAACCTTAAGAACTTAAGGTTGGCTTAATTTTGCTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATg                >  2:507642/1‑100 (MQ=21)
                                                                                     gCCAACCTTAAGAACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAg             >  2:3062755/1‑100 (MQ=255)
                                                                                     gCCAACCTTAAGAACTTAAGGTTGGC‑T‑‑TAAT‑‑TTTGCTGGATTAAGGGGAATTAATGTTTGATTATTGATGCATCAAAGGGATACGGGTATTTATGTa              >  1:424770/1‑97 (MQ=18)
                                                                                      ccAACCTTAAGAACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGc            >  2:3797057/1‑100 (MQ=255)
                                                                                      ccAACCTTAAGAACTTAAGGTTGGC‑T‑‑TAAT‑‑TTTGCTGGATTAAGGTCAATTAATCTTTGCGTATTGGTGCATCTAAGGGATACGGGGTTTTa                  >  1:3612990/1‑92 (MQ=18)
                                                                                       cAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGcc           >  1:3866018/1‑100 (MQ=255)
                                                                                       cAACCTTAAGAACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGTGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGCTTTTATGTAGc            >  1:662317/1‑99 (MQ=38)
                                                                                       cAACCTTAAGAACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGcc           >  1:2005454/1‑100 (MQ=255)
                                                                                          ccTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGt        <  2:3105572/100‑1 (MQ=255)
                                                                                           cTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGtt       >  1:1348618/1‑100 (MQ=255)
                                                                                            ttAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg      >  1:2759747/1‑100 (MQ=255)
                                                                                            ttAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg      <  2:3858173/100‑1 (MQ=255)
                                                                                            ttAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg      <  1:1781719/100‑1 (MQ=255)
                                                                                            ttAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg      <  2:235893/100‑1 (MQ=255)
                                                                                             tAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGt     >  2:2846713/1‑100 (MQ=255)
                                                                                               aGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa   <  2:3574871/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
TCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAT  >  NC_000913/4335557‑4335752

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