Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 585,035 T→C Y200C (TAC→TGC)  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,0350TC98.9% 983.8 / ‑1.7 282Y200C (TAC→TGC) ompTDLP12 prophage; outer membrane protease VII; outer membrane protein 3b
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (137/141);  total (139/143)
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.90e-01

Only 100 of 282 total aligned reads displayed.

GTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATAT  >  NC_000913/584938‑585131
                                                                                                 |                                                                                                
  gtTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAggg                                                                                              >  2:669048/1‑100 (MQ=255)
  gtTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAggg                                                                                              >  2:3574926/1‑100 (MQ=255)
  gtTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAggg                                                                                              >  1:1999208/1‑100 (MQ=255)
   ttCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGc                                                                                             <  1:4343516/100‑1 (MQ=255)
     cATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCAt                                                                                           >  2:1644816/1‑100 (MQ=255)
       tCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCCATGCAGGGCAttt                                                                                         <  1:4911125/100‑1 (MQ=255)
        cGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCAtttt                                                                                        >  2:2902172/1‑100 (MQ=255)
           tATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTaaa                                                                                     <  1:5410714/100‑1 (MQ=255)
            aTCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTaaaa                                                                                    >  1:944795/1‑100 (MQ=255)
             tCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAAc                                                                                   >  1:2925161/1‑100 (MQ=255)
             tCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAAc                                                                                   <  2:1062133/100‑1 (MQ=255)
               aGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACgt                                                                                 <  1:2250078/100‑1 (MQ=255)
                gatgatTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACgtt                                                                                >  1:3824155/1‑100 (MQ=255)
                gatgatTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACgtt                                                                                >  1:2556748/1‑100 (MQ=255)
                 atgatTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACgttg                                                                               <  1:4406805/100‑1 (MQ=255)
                     ttCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGttgtttg                                                                           <  2:1194799/100‑1 (MQ=255)
                      tCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGttgtttga                                                                          <  2:3330602/100‑2 (MQ=255)
                          cccAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGcc                                                                      <  1:1664579/100‑1 (MQ=255)
                             aGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGAt                                                                   <  1:1302163/100‑1 (MQ=255)
                               ccGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTg                                                                 >  2:373405/1‑100 (MQ=255)
                                cGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGc                                                                <  1:1540838/100‑1 (MQ=255)
                                cGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGc                                                                <  1:1873719/100‑1 (MQ=255)
                                cGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGc                                                                <  1:3267498/100‑1 (MQ=255)
                                cGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACCTTGTTTGTAGCCGATTGc                                                                <  2:768003/100‑1 (MQ=255)
                                 gCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGct                                                               <  1:148052/100‑1 (MQ=255)
                                 gCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGct                                                               >  1:1094786/1‑100 (MQ=255)
                                 gCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGct                                                               >  1:4184106/1‑100 (MQ=255)
                                  cTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGctc                                                              <  1:3314495/100‑1 (MQ=255)
                                      aTTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTtc                                                          >  1:82003/1‑100 (MQ=255)
                                      aTTTAAATGTGCCACCGAGTTCAAAATCGTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTtc                                                          <  2:2749106/100‑1 (MQ=255)
                                       tttAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTtct                                                         >  1:2747812/1‑100 (MQ=255)
                                         tAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTcc                                                       <  1:123747/100‑1 (MQ=255)
                                         tAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTcc                                                       <  1:1257998/100‑1 (MQ=255)
                                         tAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTcc                                                       <  2:2057076/100‑1 (MQ=255)
                                         tAAAGGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCCGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTcc                                                       <  1:241016/100‑1 (MQ=255)
                                            aTGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCAtt                                                    >  2:92634/1‑100 (MQ=255)
                                              gtgCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCg                                                  >  1:4082287/1‑100 (MQ=255)
                                              gtgCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCg                                                  <  1:3601530/100‑1 (MQ=255)
                                               tgCCACCGAGTTCAGAATCTTCATAAGGATAACTTCCAGTCAAGCCAGTGCAGGCCATTTTAAAACGTTGTTTGTAGCCGATTGCTATTTCTCCGTTCgg                                                 >  2:5129324/1‑100 (MQ=255)
                                               tgCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCgg                                                 <  2:3059778/100‑1 (MQ=255)
                                               tgCCACCGAGTTCAAAATCCTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCACGTCAAAACCATGTTTCTAGCCGATCGCTCTTTCTCCATTAgg                                                 >  2:4742035/1‑100 (MQ=255)
                                                tccACCGAGATCAAAATCTTCATAACGATAGCTTCGAGTCAAGCCAATGCAGGACATTTTAAAACTTTTTTTTTAGTTGATTGCTCTTTTTCCTTTCggt                                                >  2:1352857/2‑99 (MQ=255)
                                                gCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTTTTTGTAGCCGATTGCTCTTTCTCCATTCggg                                                >  2:4375594/1‑100 (MQ=255)
                                                 ccACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGa                                               <  1:3192852/100‑1 (MQ=255)
                                                  cACCTAGTTCAAAATCTTCATAACGATAACTTCCAGTTACGCAAATGCAGGGCACTTTAATCCGTGGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGaa                                              <  2:1835581/100‑1 (MQ=255)
                                                  cACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGaa                                              >  2:2403925/1‑100 (MQ=255)
                                                   aCCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAg                                             <  2:1770864/100‑1 (MQ=255)
                                                   aCCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAg                                             <  2:2119523/100‑1 (MQ=255)
                                                    ccGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAgg                                            >  1:1754303/1‑100 (MQ=255)
                                                    ccGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGACATTTTAAAACGATGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAga                                            >  2:3415681/1‑99 (MQ=255)
                                                          ttaaaaTCTTCAGAACGATAACTTAGAAACAAGCCAACGGAGGGAAATTTAAAACGTTGTTTATAGACGGTTGATATTTCTCCATTCGGGAAGAAGCCGa                                      >  2:2611677/3‑100 (MQ=255)
                                                          tCAAAATCTTCATAACGATAACTTCCCGTCAAGCCAATGCAGGGCCTTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGa                                      >  2:3992554/1‑100 (MQ=255)
                                                           cAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCGATTCGGGAAGGAGCCGat                                     <  1:4636373/100‑1 (MQ=255)
                                                           cAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGat                                     >  1:812516/1‑100 (MQ=255)
                                                           cAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGat                                     <  1:104164/100‑1 (MQ=255)
                                                           cAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGat                                     <  2:376446/100‑1 (MQ=255)
                                                           cAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGat                                     >  2:2209729/1‑100 (MQ=255)
                                                           cAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGat                                     <  2:2574114/100‑1 (MQ=255)
                                                            aaaaTCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGata                                    >  1:1967865/1‑100 (MQ=255)
                                                            aaaaTCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGata                                    <  1:4208562/100‑1 (MQ=255)
                                                             aaaTCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGatat                                   >  1:3551317/1‑100 (MQ=255)
                                                              aaTCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATatc                                  <  2:5248596/100‑1 (MQ=255)
                                                               aTCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATatca                                 >  1:4086681/1‑100 (MQ=255)
                                                               aTCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATatca                                 >  2:1122737/1‑100 (MQ=255)
                                                                tCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATatcat                                <  1:395443/100‑1 (MQ=255)
                                                                tCTTCATAACGATAACTTCCAGTCAAGCCAATGCAGCGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGTAAGGAGCAGatgttat                                >  2:1220797/1‑95 (MQ=255)
                                                                    cATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATctct                            >  1:3791012/1‑100 (MQ=255)
                                                                    cATAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATctct                            <  1:527571/100‑1 (MQ=255)
                                                                      tAACGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGa                          >  2:3238749/1‑100 (MQ=255)
                                                                       aaCGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGaa                         <  1:5217913/100‑1 (MQ=255)
                                                                        aCGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAAt                        <  1:2152760/100‑1 (MQ=255)
                                                                        aCGATAACTTCCAGTCAAGCCAATGCAGGGCATTGTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAAt                        <  1:3098484/100‑1 (MQ=255)
                                                                         cGATAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATc                       <  2:4015453/100‑1 (MQ=255)
                                                                           aTAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTGGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATccc                     <  1:2452126/100‑1 (MQ=255)
                                                                           aTAACTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATccc                     <  2:5332901/100‑1 (MQ=255)
                                                                              aCTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCcctcc                  >  1:5486583/1‑100 (MQ=255)
                                                                              aCTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCcctcc                  <  2:149714/100‑1 (MQ=255)
                                                                              aCTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCcctcc                  >  2:419300/1‑100 (MQ=255)
                                                                               cTTCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCcctcct                 <  2:4919700/100‑1 (MQ=255)
                                                                                ttCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCctcctc                <  1:756863/100‑1 (MQ=255)
                                                                                ttCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCctcctc                >  1:170680/1‑100 (MQ=255)
                                                                                 tCCAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCa               >  1:3466034/1‑100 (MQ=255)
                                                                                   cAGTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGa             >  1:2395466/1‑100 (MQ=255)
                                                                                     gTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAAc           <  2:4627014/100‑1 (MQ=255)
                                                                                     gTCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAAc           <  1:1138414/100‑1 (MQ=255)
                                                                                      tCAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACt          <  2:2340926/100‑1 (MQ=255)
                                                                                       cAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTg         <  2:569940/100‑1 (MQ=255)
                                                                                       cAAGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTg         >  2:300345/1‑100 (MQ=255)
                                                                                         aGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTa       <  1:2617642/100‑1 (MQ=255)
                                                                                         aGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTa       >  2:4294001/1‑100 (MQ=255)
                                                                                         aGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTa       >  2:1976744/1‑100 (MQ=255)
                                                                                         aGCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTa       <  1:911671/100‑1 (MQ=255)
                                                                                          gCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAg      <  1:5147435/100‑1 (MQ=255)
                                                                                          gCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAg      <  1:603199/100‑1 (MQ=255)
                                                                                          gCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAg      <  2:1934712/100‑1 (MQ=255)
                                                                                          gCCAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTGCTCCATTCGGGAAGGAGCTGATATCATCTCTGAATCCCTCCTCAGAACTGTAg      <  2:5080359/100‑1 (MQ=255)
                                                                                           ccAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGa     <  2:5383343/100‑1 (MQ=255)
                                                                                            cAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGat    >  1:4776248/1‑100 (MQ=255)
                                                                                            cAATGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGat    <  2:2217270/100‑1 (MQ=255)
                                                                                              aTGCAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGatat  >  2:1598799/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGATATCATCTCTGAATCCCTCCTCAGAACTGTAGATAT  >  NC_000913/584938‑585131

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: