Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,589,489 A→G G171G (GGT→GGC ydeT ← fimbrial usher domain‑containing protein YdeT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,589,4890AG99.7% 1402.5 / ‑4.7 339G171G (GGT→GGCydeTfimbrial usher domain‑containing protein YdeT
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (162/176);  total (163/176)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 339 total aligned reads displayed.

CACCACTCATTCCGTAATAAATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAACCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCG  >  NC_000913/1589392‑1589586
                                                                                                 |                                                                                                 
caccacTCATTCCGTAATAAATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCa                                                                                                <  2:2889727/101‑1 (MQ=255)
 accacTCATTCCGTAATAAATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCAc                                                                                               <  2:2889783/101‑1 (MQ=255)
 accacTCATTCCGTAATAAATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCAc                                                                                               >  2:2890220/1‑101 (MQ=255)
   cacTCATTCCGTAATAAATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTg                                                                                             >  2:2890225/1‑101 (MQ=255)
      tCATTCCGTAATAAATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTg                                                                                          >  2:2890229/1‑101 (MQ=255)
        aTTCCGTAATAAATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGcc                                                                                        >  1:2890233/1‑101 (MQ=255)
          tCCGTAATAAATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAg                                                                                      <  2:2889841/101‑1 (MQ=255)
          tCCGTAATAAATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAg                                                                                      >  2:2890236/1‑101 (MQ=255)
          tCCGTAATAAATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAg                                                                                      <  2:2889578/101‑1 (MQ=255)
           ccGTAATAAATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGa                                                                                     >  2:2890240/1‑101 (MQ=255)
           ccGTAATAAATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGa                                                                                     <  1:2889680/101‑1 (MQ=255)
           ccGTAATAAATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGa                                                                                     <  1:2889722/101‑1 (MQ=255)
              taataaATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGa                                                                                  >  1:2890245/1‑101 (MQ=255)
               aataaATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGAt                                                                                 <  1:2889743/101‑1 (MQ=255)
               aataaATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGAt                                                                                 <  2:2889658/101‑1 (MQ=255)
                 taaATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGt                                                                               <  1:2889668/101‑1 (MQ=255)
                  aaaTCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTa                                                                              <  1:2889695/101‑1 (MQ=255)
                   aaTGTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTAt                                                                             >  2:2890248/1‑101 (MQ=255)
                    aTCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTAtt                                                                            >  1:2890250/1‑101 (MQ=255)
                     tCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTa                                                                           <  1:2889444/101‑1 (MQ=255)
                      cTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTAc                                                                          >  1:2890253/1‑101 (MQ=255)
                       tGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTAcc                                                                         >  1:2890259/1‑101 (MQ=255)
                        ggCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCt                                                                        >  1:2890262/1‑101 (MQ=255)
                        ggCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCt                                                                        <  1:2889894/101‑1 (MQ=255)
                        ggCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCt                                                                        <  1:2889623/101‑1 (MQ=255)
                           tgctgTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCg                                                                     >  1:2890267/1‑101 (MQ=255)
                           tgctgTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCg                                                                     >  2:2890268/1‑101 (MQ=255)
                             ctgTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTg                                                                   >  2:2890277/1‑101 (MQ=255)
                             ctgTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTg                                                                   <  1:2889883/101‑1 (MQ=255)
                             ctgTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTg                                                                   <  2:2890209/101‑1 (MQ=255)
                             ctgTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTg                                                                   >  1:2890274/1‑101 (MQ=255)
                              tgTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTgg                                                                  >  1:2890278/1‑101 (MQ=255)
                               gTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTggg                                                                 >  2:2890286/1‑101 (MQ=255)
                               gTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTggg                                                                 <  2:2889801/101‑1 (MQ=255)
                               gTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTggg                                                                 <  1:2889718/101‑1 (MQ=255)
                               gTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTggg                                                                 >  1:2890284/1‑101 (MQ=255)
                                 caccacTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCa                                                                                       <  2:2889830/77‑1 (MQ=255)
                                  accacTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGgtgt                                                              >  1:2890291/1‑101 (MQ=255)
                                    cacTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTa                                                            <  2:2889948/101‑1 (MQ=255)
                                    cacTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTa                                                            >  2:2890296/1‑101 (MQ=255)
                                    cacTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTa                                                            >  1:2890295/1‑101 (MQ=255)
                                    cacTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTa                                                            <  1:2889370/101‑1 (MQ=255)
                                      cTCCGGCTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTcc                                                                      <  2:2889237/89‑1 (MQ=255)
                                      cTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTAcc                                                          >  1:2890301/1‑101 (MQ=255)
                                       tCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCg                                                         >  2:2890304/1‑101 (MQ=255)
                                        ccGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGa                                                        <  1:2889602/101‑1 (MQ=255)
                                          gACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGAcc                                                      <  1:2889806/101‑1 (MQ=255)
                                           aCTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCt                                                     >  1:2890309/1‑101 (MQ=255)
                                           aCTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCt                                                     <  2:2890087/101‑1 (MQ=255)
                                            cTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTg                                                    <  2:2889692/101‑1 (MQ=255)
                                             tGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGa                                                   <  1:2890020/101‑1 (MQ=255)
                                               tAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAAc                                                 >  2:2890313/1‑101 (MQ=255)
                                                aaCCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACg                                                >  2:2890319/1‑101 (MQ=255)
                                                aaCCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACg                                                >  1:2890316/1‑101 (MQ=255)
                                                 aCCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGc                                               >  1:2890322/1‑101 (MQ=255)
                                                  ccGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGAc                                                                                    >  2:2890323/1‑63 (MQ=255)
                                                   cGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCGCTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACg                                                <  2:2889823/98‑1 (MQ=255)
                                                    gACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCtat                                            <  1:2889775/101‑1 (MQ=255)
                                                      cATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATaa                                          <  2:2889989/101‑1 (MQ=255)
                                                        ttAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAAtt                                        >  1:2890329/1‑101 (MQ=255)
                                                        ttAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAAtt                                        >  2:2890330/1‑101 (MQ=255)
                                                           gTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAg                                     <  1:2889858/101‑1 (MQ=255)
                                                            tATTACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAgg                                    >  1:2890335/1‑101 (MQ=255)
                                                              ttACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGtt                                  <  1:2889581/101‑1 (MQ=255)
                                                              ttACCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGtt                                  >  2:2890338/1‑101 (MQ=255)
                                                                aCCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGttat                                <  2:2889705/101‑1 (MQ=255)
                                                                aCCATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGttat                                <  1:2890237/101‑1 (MQ=255)
                                                                 ccATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGttatt                               <  2:2889535/101‑1 (MQ=255)
                                                                  cATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGttatta                              >  1:2890341/1‑101 (MQ=255)
                                                                  cATAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGttatta                              <  2:2889876/101‑1 (MQ=255)
                                                                   aTAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTtattat                             <  1:2889851/101‑1 (MQ=255)
                                                                   aTAAGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTtattat                             >  1:2890345/1‑101 (MQ=255)
                                                                     aaGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATcc                           >  1:2890349/1‑101 (MQ=255)
                                                                     aaGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATcc                           <  1:2889925/101‑1 (MQ=255)
                                                                      aGCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTAtt                                                                            <  2:2889597/51‑1 (MQ=255)
                                                                       gCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCgg                         <  2:2889647/101‑1 (MQ=255)
                                                                       gCTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCgg                         >  2:2890352/1‑101 (MQ=255)
                                                                        cTCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGgc                        >  1:2890354/1‑101 (MQ=255)
                                                                         tCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATc                            >  2:2890360/1‑96 (MQ=255)
                                                                         tCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGgca                       <  2:2889742/101‑1 (MQ=255)
                                                                         tCCACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACAATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGgca                       <  2:2889813/101‑1 (MQ=255)
                                                                          ccACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGgcag                      >  1:2890361/1‑101 (MQ=255)
                                                                          ccACGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGgcag                      <  1:2889751/101‑1 (MQ=255)
                                                                             cGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGcagcag                   <  1:2890089/101‑1 (MQ=255)
                                                                             cGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGcagcag                   >  2:2890364/1‑101 (MQ=255)
                                                                             cGATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGcagcag                   <  1:2889835/101‑1 (MQ=255)
                                                                              gATAATTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTa                                                                           <  1:2889659/44‑1 (MQ=255)
                                                                               aTAATTAAGAGAACTGTAGCCACTGATGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCagag                 >  1:2890370/1‑101 (MQ=255)
                                                                                  aTTAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGc              >  2:2890372/1‑101 (MQ=255)
                                                                                    tAAGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCa            <  1:2890217/101‑1 (MQ=255)
                                                                                     aaGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCAt           >  2:2890376/1‑101 (MQ=255)
                                                                                     aaGAGAACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCAt           <  2:2889624/101‑1 (MQ=255)
                                                                                       gagaACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATaa         <  1:2889760/101‑1 (MQ=255)
                                                                                       gagaACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATaa         <  1:2889498/101‑1 (MQ=255)
                                                                                        agaACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATaaa        <  1:2889943/101‑1 (MQ=255)
                                                                                         gaACTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAAc       <  2:2889764/101‑1 (MQ=255)
                                                                                          aaCTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAAcc      <  1:2889701/101‑1 (MQ=255)
                                                                                           aCTGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAAccc     <  1:2889797/101‑1 (MQ=255)
                                                                                             tGTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAAccccc   >  1:2890381/1‑101 (MQ=255)
                                                                                              gTAGCCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTCACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCg  <  2:2889671/101‑1 (MQ=255)
                                                                                                 |                                                                                                 
CACCACTCATTCCGTAATAAATCTGGCTGCTGTCACCACTCCGACTGTAACCGACATTAGTATTACCATAAGCTCCACGATAATTAAGAGAACTGTAACCACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCG  >  NC_000913/1589392‑1589586

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