Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,059,704 A→T 20.5% intergenic (‑3/‑49) yqfE ← / → argP pseudogene, LysR family/transcriptional regulator for arginine transport and DNA replication genes; replication initiation inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,059,7040AT20.5% 398.2 / 105.7 191intergenic (‑3/‑49)yqfE/argPpseudogene, LysR family/transcriptional regulator for arginine transport and DNA replication genes; replication initiation inhibitor
Reads supporting (aligned to +/- strand):  ref base A (77/74);  new base T (19/20);  total (96/95)
Fisher's exact test for biased strand distribution p-value = 8.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-01

Only 100 of 195 total aligned reads displayed.

AATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACG  >  NC_000913/3059605‑3059802
                                                                                                   |                                                                                                  
aaTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACaaaa                                                                                                    <  1:100144/100‑1 (MQ=255)
 aTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAAt                                                                                                   <  2:1405345/100‑1 (MQ=255)
   gCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtg                                                                                                 <  1:298903/100‑1 (MQ=255)
    cATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgt                                                                                                >  2:3737520/1‑100 (MQ=255)
    cATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgt                                                                                                >  2:3256310/1‑100 (MQ=255)
     aTTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                               <  1:2682861/100‑1 (MQ=255)
     aTTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                               <  1:4133173/100‑1 (MQ=255)
     aTTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                               <  2:2261177/100‑1 (MQ=255)
      ttAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCTCGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTgtgt                                                                                              <  1:3987269/100‑1 (MQ=255)
       tAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTc                                                                                             <  1:262933/100‑1 (MQ=255)
         gcgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCag                                                                                           <  1:1992973/100‑2 (MQ=255)
         gcgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCag                                                                                           <  1:1225215/100‑2 (MQ=255)
         gcgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAt                                                                                           >  2:2360903/1‑100 (MQ=255)
         gcgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAt                                                                                           <  2:817705/100‑1 (MQ=255)
          ctcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCagg                                                                                          <  2:3226309/98‑3 (MQ=255)
          cgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGTATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCagg                                                                                          <  2:957875/100‑3 (MQ=255)
          cgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGTTTCGGTATATGCAACCTGACACAATTTTGTGTCagg                                                                                          <  1:2962812/100‑3 (MQ=255)
          cgcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATa                                                                                          >  1:1671313/1‑100 (MQ=255)
            cAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTgtgtaaggag                                                                                        >  1:1029713/1‑94 (MQ=255)
            cAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggtt                                                                                        <  2:1100196/100‑5 (MQ=255)
            cAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggag                                                                                        >  1:3226323/1‑96 (MQ=255)
             aaaTCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggt                                                                                         <  2:2532207/98‑4 (MQ=255)
               aTCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCa                                                                                     <  2:307835/100‑1 (MQ=255)
                 cccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAgg                                                                                   >  2:646705/1‑100 (MQ=255)
                 cccTGATCTTTTAAGCTGCGCTTTTAACCACGGAGAGAAATGTTACGCGCCGACGGCTTCGGTATATGCAAACTGACACAAAATTGTGTCATAGTGCAgg                                                                                   >  2:1055838/1‑100 (MQ=255)
                  ccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGa                                                                                  <  2:1902941/100‑1 (MQ=255)
                  ccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGa                                                                                  >  2:1951385/1‑100 (MQ=255)
                     gATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaaa                                                                               >  1:2189541/1‑100 (MQ=255)
                     gATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaaa                                                                               >  2:3819764/1‑100 (MQ=255)
                       tCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAg                                                                             <  2:1827002/100‑1 (MQ=255)
                       tCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAg                                                                             <  1:472760/100‑1 (MQ=255)
                       tCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAg                                                                             <  1:2277836/100‑1 (MQ=255)
                        cTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGc                                                                            >  1:591535/1‑100 (MQ=255)
                        cTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATCGTGTCATAGTGCAGGAAAAAGc                                                                            >  2:2019710/1‑100 (MQ=255)
                           ttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAtt                                                                         <  2:917366/100‑1 (MQ=255)
                            tAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAttt                                                                        >  1:3043832/1‑100 (MQ=255)
                            tAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAttt                                                                        >  1:512227/1‑100 (MQ=255)
                             aaGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTa                                                                       <  1:3167111/100‑1 (MQ=255)
                                cTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                    <  1:3738436/100‑1 (MQ=255)
                                 tGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAg                                                                   <  2:1640891/100‑1 (MQ=255)
                                 tGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGAAAAAAAAAAAAAAAAAGCATTTACCAg                                                                   <  2:4132712/100‑1 (MQ=255)
                                  gcgcTTTTACCCACGGATACAAATGTTCCCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGTAAAAAGCATTTACCAgg                                                                  <  2:1357128/100‑1 (MQ=255)
                                  gcgcTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAgg                                                                  >  2:2737753/1‑100 (MQ=255)
                                     cTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                               >  1:1745868/1‑100 (MQ=255)
                                     cTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                               >  1:2885276/1‑100 (MQ=255)
                                     cTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                               <  1:2590913/100‑1 (MQ=255)
                                       tttAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                               <  1:1271728/98‑1 (MQ=255)
                                         tAACGACGGATACAAATGTTGCCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTATAAAAAAAAAAAAAAAAGCATTTACCAGGAGCAGac                                                           <  2:2780685/100‑1 (MQ=255)
                                         tAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATCTGCAACCTGACACACAATAGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAg                                                             >  2:2327832/1‑98 (MQ=255)
                                           aCCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaa                                                         >  1:959782/1‑100 (MQ=255)
                                            ccACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaac                                                        >  2:2702752/1‑100 (MQ=255)
                                            ccACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaac                                                        >  1:2620184/1‑100 (MQ=255)
                                               cGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGc                                                     >  1:1493617/1‑100 (MQ=255)
                                                ggATACAAATGTTACTCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCa                                                    <  1:4137194/100‑1 (MQ=255)
                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                               <  2:3145329/100‑1 (MQ=255)
                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                               >  1:2532603/1‑100 (MQ=255)
                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAGTTACCAGGAGGCGACAACAGCAATGaa                                               >  2:1232119/1‑100 (MQ=255)
                                                      aaaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaaa                                              <  1:1144858/100‑1 (MQ=255)
                                                      aaaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAAAAAAAAAAAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaaa                                              <  2:2112394/100‑1 (MQ=255)
                                                         tGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACg                                            <  1:2732494/99‑1 (MQ=255)
                                                          gTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGcc                                          >  2:3473802/1‑100 (MQ=255)
                                                             aCCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCgg                                       <  2:1555365/100‑1 (MQ=255)
                                                              cccGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGa                                      >  1:3352693/1‑100 (MQ=255)
                                                               ccgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGAc                                     <  2:313122/100‑1 (MQ=255)
                                                               ccgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGAc                                     >  1:659298/1‑100 (MQ=255)
                                                                   cgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACa                                 <  2:1760536/100‑1 (MQ=255)
                                                                        gCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACa                            >  1:162455/1‑100 (MQ=255)
                                                                           tCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTa                         >  1:3689049/1‑100 (MQ=255)
                                                                             ggTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGTCAACAGCAATGAAACGCCCGGACTACAGAACATTACa                       <  2:2234426/100‑1 (MQ=255)
                                                                             ggTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACa                       >  1:2238731/1‑100 (MQ=255)
                                                                             ggTATATGCAACCTGACACAAAATTGTGTCATAGTGCAAGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACa                       <  2:1694794/100‑1 (MQ=255)
                                                                             ggTATATGCAACCCGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACa                       >  1:4097083/1‑100 (MQ=255)
                                                                              gTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAg                      <  2:696508/100‑1 (MQ=255)
                                                                              gTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAg                      <  2:39647/100‑1 (MQ=255)
                                                                                tcctgcactatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGc                    >  1:3577434/11‑100 (MQ=38)
                                                                                 tatGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCa                   >  2:4113992/1‑100 (MQ=255)
                                                                                 tatGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCa                   >  1:3695004/1‑100 (MQ=255)
                                                                                 cctgcactatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCa                   >  2:603139/10‑100 (MQ=38)
                                                                                   tgcactatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACt                 >  1:2424871/8‑100 (MQ=255)
                                                                                   tGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACt                 <  2:2610644/100‑1 (MQ=255)
                                                                                     cactatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACTAGGAGCAGACAACAGCAATGAAACGACCGAACTAAAGAACATTACAGGCACggg               >  2:4196597/6‑100 (MQ=255)
                                                                                      actatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTg                >  2:1932420/5‑98 (MQ=255)
                                                                                      actatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGAACTGGa              >  2:1030275/5‑100 (MQ=255)
                                                                                      aaCCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGa              <  2:1465114/100‑1 (MQ=255)
                                                                                       ctatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGa              >  1:2679398/4‑99 (MQ=255)
                                                                                         ctaacacAAAATTGTGTCATAGTGGAGGAAAAAGCATTTACCAGGAGCAGACAACAGCGATGAAACGCCCGGACTACATAACATTACAGGCACTGGATGc           <  2:276305/97‑1 (MQ=255)
                                                                                         cTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGc           <  2:1493614/100‑1 (MQ=255)
                                                                                         atGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGc           >  1:3119733/2‑100 (MQ=255)
                                                                                         atGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGc           >  1:474380/2‑100 (MQ=255)
                                                                                           gACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCgg         >  2:3442329/1‑100 (MQ=255)
                                                                                           gACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCgg         >  1:4068476/1‑100 (MQ=255)
                                                                                            tctcaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGt        <  2:2403018/97‑1 (MQ=255)
                                                                                             cacaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg       >  1:3469766/1‑100 (MQ=255)
                                                                                             cacaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg       >  1:3741816/1‑100 (MQ=255)
                                                                                              acaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGCTGa      >  2:1261110/1‑100 (MQ=255)
                                                                                               caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGtaga     >  1:2403023/1‑97 (MQ=255)
                                                                                               caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGAt     >  1:1378475/1‑100 (MQ=255)
                                                                                                aaaaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATa    <  1:801391/100‑1 (MQ=255)
                                                                                                 aaaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATa    <  2:3658576/99‑1 (MQ=255)
                                                                                                  aaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACg  >  2:3182715/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
AATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACG  >  NC_000913/3059605‑3059802

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: