Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,759,497 T→C 58.5% C171C (TGT→TGC yfjI → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,759,4970TC58.5% 109.0 / 191.8 178C171C (TGT→TGCyfjICP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (34/40);  new base C (44/60);  total (78/100)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 179 total aligned reads displayed.

GACTTAAAGAACTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGA  >  NC_000913/2759400‑2759595
                                                                                                 |                                                                                                  
gACTTAGAGAACTACTTACGATAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTgg                                                                                                  <  2:1778410/100‑1 (MQ=255)
gACTTAAAGAACTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTgg                                                                                                  <  2:2004442/100‑1 (MQ=255)
 aCTTAAAGAACTACTTACGACAAACCAGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGGTGCCTCACCTGCAGCTATTAACGCTCATCTCTGTggg                                                                                                 <  2:2805696/100‑1 (MQ=255)
  cTTAAAGAACTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGc                                                                                                <  1:1370870/100‑1 (MQ=255)
   ttAAAGAACTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCa                                                                                               >  2:877493/1‑100 (MQ=255)
    tAAAGAACTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCAc                                                                                              >  1:3429713/1‑100 (MQ=255)
    tAAAGAACTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCAc                                                                                              >  2:169766/1‑100 (MQ=255)
    tAAAGAACTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCAc                                                                                              >  2:3411/1‑100 (MQ=255)
       aGAACTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTgg                                                                                           >  1:1264290/1‑100 (MQ=255)
        gAACTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGc                                                                                          <  2:286952/100‑1 (MQ=255)
         aaCTACTTACGACAAACCCGACAGCTCCAGTGAGACTCAAATTTTTATTTAAAGATGCCACACCTGCCGCTATTACAGCTCATCTCTGTGGGCACTGGCg                                                                                         >  2:2211465/1‑100 (MQ=255)
         aaCTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCg                                                                                         <  1:2254386/100‑1 (MQ=255)
          actactTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGa                                                                                        <  1:3018075/100‑1 (MQ=255)
          actactTACGACAAACCCGAAAGCTCCAGCGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCAATTAAAGCTCATCTCTGCGGGCACTGGCg                                                                                         >  1:3210469/1‑99 (MQ=255)
           ctactTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTgg                                                                                           >  2:3299081/1‑96 (MQ=255)
            tactTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATc                                                                                      <  2:3883614/100‑1 (MQ=255)
            tactTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATc                                                                                      <  1:2823877/100‑1 (MQ=255)
              ctTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAg                                                                                    >  1:3815758/1‑100 (MQ=255)
                 aCGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCg                                                                                 <  1:649730/100‑1 (MQ=255)
                  cGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCgg                                                                                <  1:2957493/100‑1 (MQ=255)
                    aCAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGca                                                                              <  1:1616288/100‑1 (MQ=255)
                    aCAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGca                                                                              >  1:3256254/1‑100 (MQ=255)
                      aaaCCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGcatc                                                                            <  2:3084724/100‑1 (MQ=255)
                      aaaCCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGcatc                                                                            <  1:3685510/100‑1 (MQ=255)
                       aaCCCGAAAGCTCCATTGAGATTCAATTTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGcatca                                                                           <  2:1047268/100‑1 (MQ=255)
                       aaCCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGcatca                                                                           <  2:717020/100‑1 (MQ=255)
                       aaCCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGcatca                                                                           <  1:1578372/100‑1 (MQ=255)
                        aCCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGcatcat                                                                          >  2:2358978/1‑100 (MQ=255)
                        aCCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGcatcat                                                                          <  1:2150976/100‑1 (MQ=255)
                        aCCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGcatcat                                                                          <  2:3970069/100‑1 (MQ=255)
                         cccGAAAGCTATAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATg                                                                         <  2:1522713/100‑1 (MQ=255)
                            gAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTc                                                                       <  2:1954160/99‑1 (MQ=255)
                            gAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCt                                                                      <  2:3015169/100‑1 (MQ=255)
                             aaaGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCtg                                                                     <  1:218675/100‑1 (MQ=255)
                             aaaGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCtg                                                                     <  1:1551464/100‑1 (MQ=255)
                             aaaGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCtg                                                                     <  2:2099950/100‑1 (MQ=255)
                               aGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCtgat                                                                   <  2:514034/100‑1 (MQ=255)
                                gCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCtgatg                                                                  >  2:1667444/1‑100 (MQ=255)
                                 cTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCtgatga                                                                 >  2:1909049/1‑100 (MQ=255)
                                   ccAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAg                                                               >  2:412564/1‑100 (MQ=255)
                                    cAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGc                                                              >  2:1195018/1‑100 (MQ=255)
                                    cAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGc                                                              <  1:1107440/100‑1 (MQ=255)
                                          gaTTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAATCTCATCTCTGTGGACACCGGCGATCAGTCGGCCTCATGTCTGATGCAGCTGGGa                                                         >  2:1805389/1‑99 (MQ=255)
                                          gaTTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGAt                                                        <  1:1469721/100‑1 (MQ=255)
                                           aTTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATc                                                       >  1:701522/1‑100 (MQ=255)
                                            ttCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCa                                                      <  2:1856722/100‑1 (MQ=255)
                                             tCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCAt                                                     >  1:3480353/1‑100 (MQ=255)
                                              cAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCAtt                                                    >  1:3944232/1‑100 (MQ=255)
                                              cAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCAtt                                                    >  1:3143721/1‑100 (MQ=255)
                                              cAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCAtt                                                    >  1:630074/1‑100 (MQ=255)
                                                 aTTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCAttttt                                                 <  2:1578637/100‑1 (MQ=255)
                                                    tttatttaACCCTGACACACCTCCAGCTATTATACTTCATCTCTGTGGGCACTGGCTATCAGTCGGCATCATGTCTGATTAAGCTGGGATTATTTTTaa                                               <  1:3865361/99‑1 (MQ=255)
                                                     ttatttaACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATg                                             <  2:1792645/100‑1 (MQ=255)
                                                     ttatttaACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATg                                             <  1:544897/100‑1 (MQ=255)
                                                       atttaACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGt                                           >  1:3963014/1‑100 (MQ=255)
                                                        tttaACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGGGCGtttttgaggg                                            >  2:3387405/1‑93 (MQ=255)
                                                        tttaACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGtt                                          <  1:4050445/100‑1 (MQ=255)
                                                        tttaACGATGCCACACCTGCAGCTACTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGtt                                          <  1:2088001/100‑1 (MQ=255)
                                                         ttaACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTa                                         <  2:3628090/100‑1 (MQ=255)
                                                         ttaACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTa                                         >  1:1790341/1‑100 (MQ=255)
                                                         ttaACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTa                                         <  1:1031987/100‑1 (MQ=255)
                                                          taACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGtttc                                        >  2:2961551/1‑98 (MQ=255)
                                                          taACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTac                                        <  1:282215/100‑1 (MQ=255)
                                                           accGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTac                                        <  1:1220252/97‑1 (MQ=255)
                                                           aaCGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTaca                                       <  1:793083/100‑1 (MQ=255)
                                                            aCGATGCCAAACCGGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTacac                                      <  2:187171/100‑1 (MQ=255)
                                                             cGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTacaca                                     <  1:2248358/100‑1 (MQ=255)
                                                             cGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTacaca                                     >  1:1718389/1‑100 (MQ=255)
                                                              gATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTacacac                                    <  2:388798/100‑1 (MQ=255)
                                                              gATGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTacac                                      <  1:368032/98‑1 (MQ=255)
                                                               aTGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACt                                   <  1:1909358/100‑1 (MQ=255)
                                                               aTGCCACACCTGCAGCCAATAAAGCTCATCTCTGTGGAAACCGGCGATAAGTCCACCACAAGCCTGATACAGTTAAGCTCATCTTTAATGGTTacacac                                    <  2:2407305/99‑1 (MQ=255)
                                                                tGCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACtt                                  >  1:3358019/1‑100 (MQ=255)
                                                                 gCCACACCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTa                                 <  2:84447/100‑1 (MQ=255)
                                                                   cacaCCTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAAc                               >  1:2448448/1‑100 (MQ=255)
                                                                    gatcttGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACg                              <  2:2332650/95‑1 (MQ=255)
                                                                      acCTGCAGCTATTAAAGCTCATCTCTGTGGTCACTTGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAg                            <  1:3449193/100‑1 (MQ=255)
                                                                       ccTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGc                           >  1:1853464/1‑100 (MQ=255)
                                                                       ccTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGc                           >  1:564899/1‑100 (MQ=255)
                                                                        cTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCt                          <  1:745296/100‑1 (MQ=255)
                                                                        cTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCt                          >  1:2551846/1‑100 (MQ=255)
                                                                        cTGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCt                          <  1:357400/100‑1 (MQ=255)
                                                                         tGCAGCTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAgatc                         >  1:2332663/1‑97 (MQ=255)
                                                                          gCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGc                        >  1:815774/1‑100 (MQ=255)
                                                                             gCTATTAAAGCTCATCCCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATGTTTACTGGTTACACACTTAACGAGCTGcc                       >  2:1334780/1‑98 (MQ=255)
                                                                              cTATTAAAGCTCATCTCTGCGGGCACTGGCGATCATTTGGCATCATGTCTGATGAAGCTGCGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGtt                    <  2:3265182/100‑1 (MQ=255)
                                                                              cTATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGtt                    <  1:3506365/100‑1 (MQ=255)
                                                                               tATTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGttt                   <  1:800616/100‑1 (MQ=255)
                                                                                aTTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTa                  >  2:685753/1‑100 (MQ=255)
                                                                                aTTAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTa                  >  1:3159224/1‑100 (MQ=255)
                                                                                  tAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATc                <  1:976248/100‑1 (MQ=255)
                                                                                  tAAAGCTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATc                <  2:1264282/100‑1 (MQ=255)
                                                                                       cTCATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATaa           <  1:1437010/100‑1 (MQ=255)
                                                                                         cttctcTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGa         >  1:2793000/3‑100 (MQ=255)
                                                                                         cATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGa         <  2:2923877/100‑1 (MQ=255)
                                                                                         cATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGa         >  1:2871813/1‑100 (MQ=255)
                                                                                         cATCTCTGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGa         <  2:2464205/100‑1 (MQ=255)
                                                                                            ctctGCGGGCACTGGCGATCAGTCGGCATCATGTGTGATGAAGCTGGGATCATTTTTAATGATTACAGACGTAACGAGCTGCCGTTTATCAATAAGAt        <  2:687743/98‑1 (MQ=255)
                                                                                             tctGCGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGAtgtg     >  2:2399023/1‑100 (MQ=255)
                                                                                                gtgGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGa  >  2:366062/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
GACTTAAAGAACTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGA  >  NC_000913/2759400‑2759595

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: