Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,201,128 T→C intergenic (‑96/‑369) xisE ← / → ymfI e14 prophage, putative excisionase/e14 prophage, uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,201,1280TC91.1% 392.1 / 18.7 123intergenic (‑96/‑369)xisE/ymfIe14 prophage, putative excisionase/e14 prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (6/5);  new base C (58/54);  total (64/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

Only 100 of 123 total aligned reads displayed.

TCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGTGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTGGTACCGCCAAGACTACACACAGCATAAA  >  NC_000913/1200995‑1201260
                                                                                                                                     |                                                                                                                                    
tCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGTGGCca                                                                                                                                 <  1:827468/139‑1 (MQ=255)
tCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCca                                                                                                                                 <  1:217512/139‑1 (MQ=255)
         cATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCa                                                                                                                        >  2:189543/1‑139 (MQ=255)
                      ttAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGt                                                                                                                >  2:813570/1‑134 (MQ=255)
                      ttAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTc                                                                                                               >  2:377793/1‑135 (MQ=255)
                      ttAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGt                                                                                                           >  1:225519/1‑139 (MQ=255)
                        aaCATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTa                                                                                                                           >  2:1479776/1‑121 (MQ=255)
                        aaCATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGtta                                                                                                         >  1:1353303/1‑139 (MQ=255)
                             ttGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAAc                                                                                                    <  2:1248759/139‑1 (MQ=255)
                                     tatGCGATCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGAGTTATTAACCCGTCATc                                                                                            <  2:1439495/139‑1 (MQ=255)
                                        gCGATCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCgg                                                                                         <  2:358014/139‑1 (MQ=255)
                                        gCGATCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCgg                                                                                         <  1:626484/139‑1 (MQ=255)
                                           aTCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGTGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGa                                                                                      <  1:996722/139‑1 (MQ=255)
                                            tCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCgg                                                                                                             >  2:838507/1‑115 (MQ=255)
                                            tCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCgg                                                                                                             <  1:838507/115‑1 (MQ=255)
                                             cACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGTGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCg                                                                                                 >  1:1251762/1‑126 (MQ=255)
                                             cACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGttatt                                                                                                       <  1:585878/120‑1 (MQ=255)
                                             cACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGttatt                                                                                                       >  2:585878/1‑120 (MQ=255)
                                              aCCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATAc                                                                                   >  1:94896/1‑139 (MQ=255)
                                               ccTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATAct                                                                                  <  1:189543/139‑1 (MQ=255)
                                                    ttaCTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTATCCCGTCATCCGGGGATACtcttct                                                                             <  2:94896/139‑1 (MQ=255)
                                                        tacACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGt                                                                         >  2:497751/1‑139 (MQ=255)
                                                        tacACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGt                                                                         <  1:920025/139‑1 (MQ=255)
                                                          cacTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGttt                                                                       >  2:1180105/1‑139 (MQ=255)
                                                          cacTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGttt                                                                       >  2:467726/1‑139 (MQ=255)
                                                             tAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGc                                                                                                                                   <  1:147796/76‑1 (MQ=255)
                                                              aaCTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGcc                                                                                                                                  <  1:814859/76‑1 (MQ=255)
                                                              aaCTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGcc                                                                                                                                  >  2:814859/1‑76 (MQ=255)
                                                              aaCTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCcac                                                                                                                                <  1:680380/78‑1 (MQ=255)
                                                              aaCTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTa                                                                   >  1:331519/1‑139 (MQ=255)
                                                               aCTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCcaca                                                                                                                               <  2:613923/78‑1 (MQ=255)
                                                               aCTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCca                                                                                                                                 >  2:1096826/1‑76 (MQ=255)
                                                               aCTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCca                                                                                                                                 <  1:1096826/76‑1 (MQ=255)
                                                               aCTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAg                                                                                                                              <  1:1361260/79‑1 (MQ=255)
                                                               aCTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAg                                                                                                                              >  2:1361260/1‑79 (MQ=255)
                                                               aCTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGt                                                                                                >  1:1364761/1‑109 (MQ=255)
                                                               aCTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACtcttct                                                                             >  1:1402458/1‑128 (MQ=255)
                                                               aCTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACtcttct                                                                             <  2:1402458/128‑1 (MQ=255)
                                                               aCTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTaa                                                                  >  1:263126/1‑139 (MQ=255)
                                                                  gCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTaaaaa                                                               >  2:1564751/1‑139 (MQ=255)
                                                                  gCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTaaaaa                                                               >  2:915154/1‑139 (MQ=255)
                                                                     tAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAggg                                                            <  2:1353303/139‑1 (MQ=255)
                                                                      aGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGc                                                           <  1:1424083/139‑1 (MQ=255)
                                                                       gTCTCAGCATGTCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAggg                                                            <  2:1341829/137‑1 (MQ=255)
                                                                       gTCTCAGCATGTCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAggg                                                            >  1:1341829/1‑137 (MQ=255)
                                                                        tctcAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCgg                                                         <  2:916130/139‑1 (MQ=255)
                                                                         ctcAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGTGGCcac                                                                                                                                >  2:1160291/1‑67 (MQ=255)
                                                                          tcAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTa                                                       <  2:225519/139‑1 (MQ=255)
                                                                          tcAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTa                                                       <  2:331519/139‑1 (MQ=255)
                                                                           cAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCca                                                                                                                                 <  2:1065367/64‑1 (MQ=255)
                                                                           cAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCca                                                                                                                                 >  1:1065367/1‑64 (MQ=255)
                                                                           cAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTAc                                                      <  2:1279411/139‑1 (MQ=255)
                                                                             gCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTAt                                                                                                                          <  1:639442/69‑1 (MQ=255)
                                                                             gCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTAt                                                                                                                          >  2:639442/1‑69 (MQ=255)
                                                                             gCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGCTGAGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCa                                                    >  2:1010059/1‑139 (MQ=255)
                                                                             gCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGttatt                                                                                                       >  2:1531228/1‑88 (MQ=255)
                                                                             gCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGttatt                                                                                                       <  1:1531228/88‑1 (MQ=255)
                                                                               atatCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGa                                                  <  2:263126/139‑1 (MQ=255)
                                                                                 atCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCcac                                                                                                                                <  2:1460973/59‑1 (MQ=255)
                                                                                 atCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCcac                                                                                                                                >  1:1460973/1‑59 (MQ=255)
                                                                                 atCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGaaa                                                >  1:164463/1‑139 (MQ=255)
                                                                                  tCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTc                                                                                               <  1:1553438/91‑1 (MQ=255)
                                                                                  tCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTc                                                                                               >  2:1553438/1‑91 (MQ=255)
                                                                                  tCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATa                                                                                    >  1:1571552/1‑102 (MQ=255)
                                                                                        cTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATa                                                                                    <  1:309775/96‑1 (MQ=255)
                                                                                        cTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTaa                                         >  1:484539/1‑139 (MQ=255)
                                                                                                tcggtTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTGCCAGAAAGGACTAAGGaaaaaa                                 >  2:394718/1‑139 (MQ=255)
                                                                                                tcggtTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGaaaaaa                                 >  1:344997/1‑139 (MQ=255)
                                                                                                tcagtTAGTTTCTCCACAAGACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGaaaaaa                                 >  2:1096437/4‑139 (MQ=255)
                                                                                                   gtTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCcac                                                                                                                                >  1:458556/1‑41 (MQ=255)
                                                                                                   gtTAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCcac                                                                                                                                <  2:458556/41‑1 (MQ=255)
                                                                                                    ttAGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAggg                                                            <  1:1054333/108‑1 (MQ=255)
                                                                                                      aGTTTCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTGCCAGAAAGGACTAAGGAAAAAACTGGTa                           <  1:394718/139‑1 (MQ=255)
                                                                                                        tttCTCCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTGGTAcc                         <  1:497751/139‑1 (MQ=255)
                                                                                                          tctcCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTGGt                            <  1:865205/134‑1 (MQ=255)
                                                                                                          tctcCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTGGTACCGc                       <  1:1180105/139‑1 (MQ=255)
                                                                                                            tcCACAAAACAGAGAAGAGCACCTGTGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTGGTACCGCCa                     >  2:1648357/1‑139 (MQ=255)
                                                                                                            tcCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTaaaa                                                                <  1:1169257/96‑1 (MQ=255)
                                                                                                            tcCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTaaaa                                                                >  2:1169257/1‑96 (MQ=255)
                                                                                                            tcCACAAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTGGTACCGCCa                     <  1:1564751/139‑1 (MQ=255)
                                                                                                              cacaAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTctct                                                                           <  1:518464/83‑1 (MQ=255)
                                                                                                              cacaAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTctct                                                                           >  2:518464/1‑83 (MQ=255)
                                                                                                              cacaAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTgg                             <  1:127254/129‑1 (MQ=255)
                                                                                                              cacaAAACAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTgg                             >  2:127254/1‑129 (MQ=255)
                                                                                                                  aaaCAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTGGTACCGCCAAGACTa               >  1:555703/1‑139 (MQ=255)
                                                                                                                     cAGAGAAGAGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGa                                                  <  1:416811/101‑1 (MQ=255)
                                                                                                                        agaagaGCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTGGTACCGCCAAGACTAcacaca         >  1:1440330/1‑139 (MQ=255)
                                                                                                                             agCACCTGCGGCCACAGCTATCAGGATGGGTCGGGttatt                                                                                                       >  2:918540/1‑40 (MQ=255)
                                                                                                                             agCACCTGCGGCCACAGCTATCAGGATGGGTCGGGttatt                                                                                                       <  1:918540/40‑1 (MQ=255)
                                                                                                                             agCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATc                                                                                            <  2:309452/51‑1 (MQ=255)
                                                                                                                             agCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATc                                                                                            >  1:309452/1‑51 (MQ=255)
                                                                                                                             agCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATa                                                                                    >  2:1409120/1‑59 (MQ=255)
                                                                                                                             agCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTGGTACCGCCAAGACTACACACAGCATa    >  1:1153686/1‑139 (MQ=255)
                                                                                                                             agCACCTGCGGCCACAGCTATCAGGATGGGTCGGTTATTAACCCGTCATCCgg                                                                                         <  1:1055801/53‑1 (MQ=255)
                                                                                                                             agCACCTGCGGCCACAGCTATCAGGATGGGTCGGTTATTAACCCGTCATCCgg                                                                                         >  2:1055801/1‑53 (MQ=38)
                                                                                                                              gCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTGGTACCGCCa                     >  1:429151/1‑121 (MQ=255)
                                                                                                                              gCACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTGGTACCGCCa                     <  2:429151/121‑1 (MQ=255)
                                                                                                                               cACCTGTGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCgggg                                                                                       >  2:729283/1‑54 (MQ=255)
                                                                                                                               cACCTGTGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCgggg                                                                                       <  1:729283/54‑1 (MQ=255)
                                                                                                                               cACCTGCGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTGGTACCGCCAAGACTACACACAGCATaaa  >  2:449682/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                    
TCAGCAGCCCATTCTTCAAGAGTTAACATTTGAAGCATATGCGATCACCTTATTACTACACTAACTGCTTAGTCTCAGCATATCGACCCTGCACGGTCGGTTAGTTTCTCCACAAAACAGAGAAGAGCACCTGTGGCCACAGCTATCAGGATGGGTCGGGTTATTAACCCGTCATCCGGGGATACTCTTCTCTGTTTTGTAAAAAGGGCGGTACCAGAAAGGACTAAGGAAAAAACTGGTACCGCCAAGACTACACACAGCATAAA  >  NC_000913/1200995‑1201260

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