Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,197,739 A→G 100% S167P (TCC→CCC)  ymfE ← e14 prophage; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,7390AG98.2% 430.3 / ‑4.2 131S167P (TCC→CCC) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/5);  new base G (69/56);  total (69/62)
Fisher's exact test for biased strand distribution p-value = 2.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

TGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTA  >  NC_000913/1197644‑1197834
                                                                                               |                                                                                               
tGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAtt                                                                                             <  2:2283225/100‑1 (MQ=255)
tGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAtt                                                                                             <  2:4130521/100‑1 (MQ=255)
  agaATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCAATTATTAAAAAAAAAAGATTATTAGGAAGGGAAt                                                                                              <  1:1634815/97‑1 (MQ=255)
  agaATGATAAAATGACAAAAAAGAAGAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTaa                                                                                           <  1:230666/100‑1 (MQ=255)
   gaATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAAt                                                                                          >  2:1188589/1‑100 (MQ=255)
    aaTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATg                                                                                         >  1:581213/1‑100 (MQ=255)
    aaTGATAAAATTACAAAAAAGAATAAAGACAGAAAAGAAACTTAAACATGAAAAATATACATTTGTTAGAAAAAAAAGATTATTAGGAAGGGAATTAATg                                                                                         <  1:865353/100‑1 (MQ=255)
     aTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGc                                                                                        <  2:1649446/100‑1 (MQ=255)
     aTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGc                                                                                        <  2:1301511/100‑1 (MQ=255)
       gATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCaa                                                                                      >  2:2508183/1‑100 (MQ=255)
        aTAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAAt                                                                                     <  2:448691/100‑1 (MQ=255)
         tAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAAtt                                                                                    >  2:253229/1‑100 (MQ=255)
           aaaTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTAt                                                                                  >  2:3052326/1‑100 (MQ=255)
             aTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCg                                                                                >  1:2436088/1‑100 (MQ=255)
               tACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaa                                                                              <  2:181474/100‑1 (MQ=255)
                aCAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaaa                                                                             >  1:722800/1‑100 (MQ=255)
                 cAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaaaa                                                                            >  1:1689440/1‑100 (MQ=255)
                  aaaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAt                                                                           >  2:116391/1‑100 (MQ=255)
                   aaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAtt                                                                          >  2:33907/1‑100 (MQ=255)
                    aaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATAACCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaaaa                                                                            >  1:490262/1‑97 (MQ=255)
                    aaaaGAATAAAGACAGCAAAGAAAAAAAAAAAAAAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTc                                                                         <  1:1583728/100‑1 (MQ=255)
                            aaaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAata                                                                 <  2:3759931/100‑1 (MQ=255)
                             aaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAatat                                                                >  1:1233714/1‑100 (MQ=255)
                              aGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATc                                                               >  1:2019757/1‑100 (MQ=255)
                                ataataaaaaaCCATCAATATAAAAAATATACATTTGTTTGCAAAAAAAAATTATTAGGAAAGGAATTAATTCAATTATAAAAAATAAAAAAAATATCCa                                                             <  1:3421464/95‑1 (MQ=255)
                                 cAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaa                                                            >  2:4102074/1‑100 (MQ=255)
                                   gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                          >  2:1347354/1‑100 (MQ=255)
                                   gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                          >  2:634516/1‑100 (MQ=255)
                                      aaGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATa                                                       <  2:2311121/100‑1 (MQ=255)
                                       aGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAg                                                      <  2:3317901/100‑1 (MQ=255)
                                        naaCCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGt                                                     >  1:508167/2‑100 (MQ=255)
                                           ccTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGtata                                                  >  2:292706/1‑100 (MQ=255)
                                            cTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGtata                                                  >  2:1646076/1‑99 (MQ=255)
                                              tCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACtt                                               <  2:3175589/100‑1 (MQ=255)
                                               cAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACttt                                              >  2:2255642/1‑100 (MQ=255)
                                                aagaTGAAAAATATGCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACtt                                               <  2:4175471/95‑1 (MQ=255)
                                                  cATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTAtt                                           <  2:292352/100‑1 (MQ=255)
                                                  cATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTAtt                                           >  2:2757693/1‑100 (MQ=255)
                                                   aTGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTc                                          >  1:2753952/1‑100 (MQ=255)
                                                    tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTcc                                         <  2:3157175/100‑1 (MQ=255)
                                                    tGAAAAATATACATTTGTTTACAAAAAAAGATTATTAAGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTcc                                         <  1:1386434/100‑1 (MQ=255)
                                                     gAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCa                                        >  2:3131149/1‑100 (MQ=255)
                                                     gAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCa                                        >  2:2951917/1‑100 (MQ=255)
                                                      aaaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCag                                       >  1:3235109/1‑100 (MQ=255)
                                                      aaaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCa                                        <  1:1391475/99‑1 (MQ=255)
                                                      aaaaaTATCCATTTGAATGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCag                                       <  2:3155272/100‑1 (MQ=255)
                                                       aaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCaga                                      >  2:190873/1‑100 (MQ=255)
                                                       aaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCaga                                      <  2:2973919/100‑1 (MQ=255)
                                                         aaaaaccATTTGATAGCAAAAAAAGATTATTAGGAAAAGAATTAATGCAATTAAAGAAAAAAAAAAAAATATCCAAAAATAGTATACTTTATTCCagaag                                    <  1:253228/95‑1 (MQ=255)
                                                         aaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCagaag                                    <  2:4083504/100‑1 (MQ=255)
                                                         aaTATCCATTTGTTAGAAAAAAAAGATTAATAGGAAGGGAATTAATGCAATTATCGAAAATTCAAACAATATCCAAAAATAGTATACTTTATTCCagaag                                    <  1:2571470/100‑1 (MQ=255)
                                                           tatCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAagag                                  <  2:3908765/100‑1 (MQ=255)
                                                             tCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGtt                                >  2:3260751/1‑100 (MQ=255)
                                                              ccATTTGTTTGCAAAATAAGATTCTTAGGAAGGGAATTAATGCAATTATCGCACATTCAAAAAATATCCAATAATAGTATACTATATTCCagaag                                    >  1:1109747/1‑95 (MQ=255)
                                                              ccATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTc                               >  2:1474419/1‑100 (MQ=255)
                                                              ccATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTc                               >  2:2168446/1‑100 (MQ=255)
                                                               cATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCa                              >  2:782768/1‑100 (MQ=255)
                                                               cATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCa                              >  1:3897527/1‑100 (MQ=255)
                                                               cATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCa                              >  1:2875642/1‑100 (MQ=255)
                                                                 tttgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAat                            >  1:3205783/1‑100 (MQ=255)
                                                                 tttgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAat                            <  1:1795194/100‑1 (MQ=255)
                                                                   tgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAatat                          >  2:1762270/1‑100 (MQ=255)
                                                                   tgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAatat                          >  2:2499801/1‑100 (MQ=255)
                                                                     tttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATaa                        >  1:2054648/1‑100 (MQ=255)
                                                                       tgCAAACAAAGATTATTAGGAAGGCAATTAATGCAATTATCGAAAATTAAACAACTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATCTAAtt                      <  1:1916909/100‑2 (MQ=255)
                                                                       tgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATg                      >  1:100166/1‑100 (MQ=255)
                                                                        gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                     >  2:2464083/1‑100 (MQ=255)
                                                                        gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                     >  2:1050162/1‑100 (MQ=255)
                                                                        gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATACAAAAATAGTATACTTTATTCCAGAAGAGATCAATATAATGt                     >  1:712110/1‑100 (MQ=255)
                                                                         cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                    <  2:2689299/100‑1 (MQ=255)
                                                                         cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                    >  2:3096964/1‑100 (MQ=255)
                                                                         cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                    >  1:2082805/1‑100 (MQ=255)
                                                                         cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCAAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                    >  1:623924/1‑100 (MQ=255)
                                                                         cAAAAAAAGATTATTAGGAAGGGAATTAATCCAATTATCGAAAATTCAAAAAATATCCAAAACTAGTATACTTTATTCCAAAAGAGTTCAATATAATGtt                    >  1:2766485/1‑100 (MQ=255)
                                                                         cAAAAAAAGATTATGAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAat                            <  1:455144/92‑1 (MQ=255)
                                                                          aaaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGttt                   >  1:1128057/1‑100 (MQ=255)
                                                                          aaaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGttt                   >  2:2457499/1‑100 (MQ=255)
                                                                          aaaaaaaGATTATTAAGAAAAAAATTAATGCAATTATCTAAAATTAAAAAAAAAAAAAAAAATAGTATACTTTATTCCAGAagag                                  <  1:2514615/85‑1 (MQ=255)
                                                                           aaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                  <  2:2928816/100‑1 (MQ=255)
                                                                           aaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                  <  2:785329/100‑1 (MQ=255)
                                                                              aaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCt               <  1:399848/100‑1 (MQ=255)
                                                                              aaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCt               >  1:23633/1‑100 (MQ=255)
                                                                               aaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCtt              >  1:349254/1‑100 (MQ=255)
                                                                                aGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTc             <  1:2132233/100‑1 (MQ=255)
                                                                                    tattaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAtt         >  1:3040032/1‑100 (MQ=255)
                                                                                    tattaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAtt         <  2:1053280/100‑1 (MQ=255)
                                                                                     attaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttt        <  2:1833053/100‑1 (MQ=255)
                                                                                      ttaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAtttt       <  2:3375433/100‑1 (MQ=255)
                                                                                       taGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttttt      <  2:1512725/100‑1 (MQ=255)
                                                                                       taGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttttt      <  1:1280431/100‑1 (MQ=255)
                                                                                       taGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttttt      <  2:1512351/100‑1 (MQ=255)
                                                                                        aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc     <  2:3765578/100‑1 (MQ=255)
                                                                                        aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc     >  1:3567351/1‑100 (MQ=255)
                                                                                        aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc     <  2:2433601/100‑1 (MQ=255)
                                                                                        aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc     <  2:2402482/100‑1 (MQ=255)
                                                                                         ggaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCt    >  1:205473/1‑100 (MQ=255)
                                                                                          gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt   >  1:2083370/1‑100 (MQ=255)
                                                                                          gaagggaATTAATGCAATTAAAAAAAAAAAAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt   <  1:2891432/100‑1 (MQ=255)
                                                                                           aagggaATTAATGCAATTATTGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTa  >  1:2611888/1‑100 (MQ=255)
                                                                                           aagggaATTAATGCAATTATCGAAAATTAAAAAAATATACAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCtt   <  1:2951905/99‑1 (MQ=255)
                                                                                               |                                                                                               
TGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTA  >  NC_000913/1197644‑1197834

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