Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 585,077 C→T 23.5% G186E (GGA→GAA)  ompT ← DLP12 prophage; outer membrane protease VII; outer membrane protein 3b

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913585,0770CT23.5% 142.3 / 44.4 81G186E (GGA→GAA) ompTDLP12 prophage; outer membrane protease VII; outer membrane protein 3b
Reads supporting (aligned to +/- strand):  ref base C (28/34);  new base T (10/9);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-02

TTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATC  >  NC_000913/584978‑585176
                                                                                                   |                                                                                                   
ttAAATGTGCCTCCGAGTTCAAAATCTTGATAACGATAGCTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTtctc                                                                                                     <  2:1019887/100‑1 (MQ=255)
 tAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTtctc                                                                                                     >  1:1377382/1‑99 (MQ=255)
 tAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTcc                                                                                                    <  1:707717/100‑1 (MQ=255)
  aaaTGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATCTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCa                                                                                                   <  1:206916/100‑1 (MQ=255)
     tgtgCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTc                                                                                                >  1:307127/1‑100 (MQ=255)
     tgtgCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTc                                                                                                <  1:1848573/100‑1 (MQ=255)
      gtgCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCAATTGCTCTTTCTTCATTCg                                                                                               <  1:330349/100‑1 (MQ=255)
       tgCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTTCATTCgg                                                                                              <  1:1442271/100‑1 (MQ=255)
        gCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCggg                                                                                             <  1:9375/100‑1 (MQ=255)
        gCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCggg                                                                                             <  2:1989107/100‑1 (MQ=255)
         ccACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCggg                                                                                             >  2:1635129/1‑99 (MQ=255)
          cACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGaa                                                                                           >  1:511567/1‑100 (MQ=255)
          cACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGa                                                                                            >  1:2230428/1‑99 (MQ=255)
           aacaaattCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAg                                                                                          <  2:1941034/94‑1 (MQ=255)
           aCCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAg                                                                                          <  1:1645586/100‑1 (MQ=255)
              gAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTTCATTCGGGAAGGAg                                                                                       <  1:1451890/100‑1 (MQ=255)
                gTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAg                                                                                       >  2:311133/1‑98 (MQ=255)
                 ttCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCg                                                                                    >  1:1815511/1‑100 (MQ=255)
                  tCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTTCATTCGGGAAGGAGCCGa                                                                                   <  1:161557/100‑1 (MQ=255)
                   cAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTTCATTCGGGAAGGAGCCGat                                                                                  >  2:1200962/1‑100 (MQ=255)
                   cAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGat                                                                                  <  1:1193825/100‑1 (MQ=255)
                      aaTCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGatat                                                                                >  2:2390079/1‑99 (MQ=255)
                         cTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATatcatc                                                                            >  1:1116250/1‑100 (MQ=255)
                          ttCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTTCATTCGGGAAGGAGCCGATATCAtct                                                                           >  1:1418899/1‑100 (MQ=255)
                          ttCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCAtct                                                                           >  1:700246/1‑100 (MQ=255)
                           tCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGTAGCCGATATCAtctc                                                                          >  2:247008/1‑100 (MQ=255)
                           tCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCAtctc                                                                          <  1:342752/100‑1 (MQ=255)
                           tCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAGTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCAtctc                                                                          <  2:232846/100‑1 (MQ=255)
                             aTAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTTCATTCGGGAAGGAGCCGATATCATCTCTg                                                                        <  2:67682/100‑1 (MQ=255)
                             aTAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTg                                                                        <  1:1872565/100‑1 (MQ=255)
                             aTAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTg                                                                        >  1:146989/1‑100 (MQ=255)
                              tAACGATAATTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGa                                                                       <  1:605497/100‑1 (MQ=255)
                              tAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGa                                                                       <  1:1549025/100‑1 (MQ=255)
                               aaCGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGaa                                                                      <  1:2349071/100‑1 (MQ=255)
                                 cGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAAt                                                                     <  1:124449/99‑1 (MQ=255)
                                 cGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATc                                                                    <  2:1963515/100‑1 (MQ=255)
                                      aCTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTTCATTCGGGAAGGAGCCGATATCATCTCTGAATCcctcc                                                               <  1:556117/100‑1 (MQ=255)
                                      aCTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCcctcc                                                               >  2:935105/1‑100 (MQ=255)
                                       cTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCcctcct                                                              >  1:2077947/1‑100 (MQ=255)
                                        ttacAGTCAAGCCAATGTAGGGTATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCctcctc                                                             <  1:438014/97‑1 (MQ=255)
                                        ttCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCctcctc                                                             <  1:1741705/100‑1 (MQ=255)
                                          ccAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCctcctc                                                             >  1:955029/1‑98 (MQ=255)
                                          ccAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAg                                                           <  2:909033/100‑1 (MQ=255)
                                           cAGTCGAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAACCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGa                                                          <  2:1725639/100‑1 (MQ=255)
                                             gTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAAc                                                        >  1:969722/1‑100 (MQ=255)
                                             gTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAAc                                                        <  2:494841/100‑1 (MQ=255)
                                               cAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTg                                                      >  2:605981/1‑100 (MQ=255)
                                                aaGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGt                                                     <  2:2230402/100‑1 (MQ=255)
                                                 aGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTTCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTa                                                    >  2:1403239/1‑100 (MQ=255)
                                                 aGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTTCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTa                                                    <  1:1403254/100‑1 (MQ=255)
                                                 aGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTa                                                    <  2:1864259/100‑1 (MQ=255)
                                                 aGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTa                                                    <  2:655259/100‑1 (MQ=255)
                                                  gCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAg                                                   >  2:524229/1‑100 (MQ=255)
                                                    cAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGat                                                 <  1:524232/100‑1 (MQ=255)
                                                       tGTATGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTTTTGCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGat                                                 <  2:2254255/97‑1 (MQ=255)
                                                       tGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGAtata                                              <  1:727884/100‑1 (MQ=255)
                                                         tAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAgg                                            >  2:2326449/1‑100 (MQ=255)
                                                         tAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAgg                                            >  1:1671478/1‑100 (MQ=255)
                                                           gggCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGaa                                          <  1:2245413/100‑1 (MQ=255)
                                                            ggCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAac                                         <  2:1724403/100‑1 (MQ=255)
                                                             gCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTTCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAacc                                        >  2:1642854/1‑100 (MQ=255)
                                                             gCATTTTAAAACGTTGTTTGTAGCCGATTACTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAacc                                        <  2:734989/100‑1 (MQ=255)
                                                              cATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAac                                         <  1:285516/98‑1 (MQ=255)
                                                                ttttAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAaccacc                                     >  2:1889311/1‑100 (MQ=255)
                                                                     aaaCGTTGTTTGTAGCCGATTGCTCTATCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAAACACctc                                   >  2:1472434/1‑97 (MQ=255)
                                                                      aaCGGTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACGGTAGATATAGGAACCACCTCTGGc                               >  2:1068040/1‑100 (MQ=255)
                                                                       aCGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCt                              >  1:1970577/1‑100 (MQ=255)
                                                                        gatcttttgtAGCCGATTGCTCTTTCTTCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTg                             <  1:882389/95‑1 (MQ=255)
                                                                        cGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTg                             <  2:810154/100‑1 (MQ=255)
                                                                        cGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTg                             <  1:1923140/100‑1 (MQ=255)
                                                                          ttgtttgtAGCCGATTGCTCTTTCTTCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTa                           <  1:1494521/100‑1 (MQ=255)
                                                                           tgtttgtAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTaa                          <  2:1980807/100‑1 (MQ=255)
                                                                             tttgtAGCCGATTGCTCTTTCTTCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTgagatc                        >  2:882379/1‑95 (MQ=255)
                                                                             tttgtAGCCGATTGCTCTTTCTTCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAg                        >  1:1058567/1‑100 (MQ=255)
                                                                             tttgtAGCCGATTGCTCTTTCTTCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAg                        >  1:1058655/1‑100 (MQ=255)
                                                                                gtAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCta                     >  1:2164542/1‑100 (MQ=255)
                                                                                       atggCTCTTTCTCCATTCTGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGc              <  1:794302/97‑1 (MQ=255)
                                                                                          gCTCTTTCTTCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCttt           >  2:2032358/1‑100 (MQ=255)
                                                                                          gCTCTTTCTTCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCttt           >  2:1602386/1‑100 (MQ=255)
                                                                                          gCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCttt           >  1:2053242/1‑100 (MQ=255)
                                                                                            tctTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTcc         >  1:1839653/1‑100 (MQ=255)
                                                                                             ctTTCTTCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCt        >  2:2362408/1‑100 (MQ=255)
                                                                                              tttCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTg       >  1:1512332/1‑100 (MQ=255)
                                                                                              tttCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTg       >  1:1512313/1‑100 (MQ=255)
                                                                                                 ctcCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGata    >  2:1786958/1‑100 (MQ=255)
                                                                                                 ctcCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGata    >  1:1292233/1‑100 (MQ=255)
                                                                                                  tcCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGatat   >  1:1838898/1‑100 (MQ=255)
                                                                                                  tcCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGatat   <  2:1108706/100‑1 (MQ=255)
                                                                                                   ccATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATc  <  2:2375185/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATATAGGAACCACCTCTGGCTGTAAAGCTATAACGGCTTTCCTGATATC  >  NC_000913/584978‑585176

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: