Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,518,507 G→A 44.7% R14H (CGC→CAC)  yfeC → DUF1323 family putative DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,518,5070GA44.7% 77.9 / 226.8 181R14H (CGC→CAC) yfeCDUF1323 family putative DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (47/53);  new base A (39/42);  total (86/95)
Fisher's exact test for biased strand distribution p-value = 8.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 193 total aligned reads displayed.

AAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCAATGAACA  >  NC_000913/2518409‑2518606
                                                                                                  |                                                                                                   
 aaaTGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAgccgcc                                                                                                   <  2:185116/100‑1 (MQ=255)
 aaaTGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAgccgcc                                                                                                   >  1:2604217/1‑100 (MQ=255)
 aaaTGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAgccgcc                                                                                                   <  2:3241757/100‑1 (MQ=255)
 aaaTGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAgccacc                                                                                                   >  2:3523813/1‑100 (MQ=255)
 aaaTGAACTACTAATAGACCCACATACATGCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAgccacc                                                                                                   <  1:1455491/100‑1 (MQ=255)
    tGAACTAGTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGa                                                                                                >  2:2232902/1‑100 (MQ=255)
       actactAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCa                                                                                             >  1:60997/1‑100 (MQ=255)
        ctactaATAGACCCACATACATTCAGGGCATTGTTGTGATCAAGGAGTGGGTGACACCAGATGAACTTGCCAGACTGAACAGATATAGCCGCCAGACCAt                                                                                            >  1:2742495/1‑100 (MQ=255)
        ctactaATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCAt                                                                                            >  2:3096728/1‑100 (MQ=255)
        ctactaATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCAt                                                                                            <  2:606106/100‑1 (MQ=255)
          actaATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCagat                                                                                          >  2:4124302/1‑97 (MQ=255)
          actaATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCAttn                                                                                          <  1:3248405/100‑2 (MQ=255)
          actaATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATta                                                                                          >  1:578416/1‑100 (MQ=255)
          actaATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATta                                                                                          <  1:3532574/100‑1 (MQ=255)
           ctaATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATtaa                                                                                         >  1:1076946/1‑100 (MQ=255)
            taATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATtaat                                                                                        >  1:987468/1‑100 (MQ=255)
             aaTAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATtaata                                                                                       >  2:1966040/1‑100 (MQ=255)
             aaTAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATtaata                                                                                       >  2:1770083/1‑100 (MQ=255)
             aaTAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATtaata                                                                                       >  1:313693/1‑100 (MQ=255)
              aTAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATtaataa                                                                                      >  1:1216132/1‑100 (MQ=255)
               tAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATaaa                                                                                     >  1:2850855/1‑100 (MQ=255)
                aGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAAt                                                                                    >  2:3115290/1‑100 (MQ=255)
                 gACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATg                                                                                   <  1:1556391/100‑1 (MQ=255)
                 gACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATg                                                                                   >  1:3019955/1‑100 (MQ=255)
                   cccACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATggg                                                                                 <  1:2753774/100‑1 (MQ=255)
                   cccACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATggg                                                                                 <  1:3148926/100‑1 (MQ=255)
                   cccACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATggg                                                                                 <  2:1333331/100‑1 (MQ=255)
                   cccACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATggg                                                                                 <  1:161996/100‑1 (MQ=255)
                    ccACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATggg                                                                                 <  1:1230578/99‑1 (MQ=255)
                     cacaTACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTa                                                                               >  2:3872828/1‑100 (MQ=255)
                       catacatTCAGGGAATTGTTATGTTCAAGGAGCGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACg                                                                             >  2:4056502/1‑100 (MQ=255)
                       catacatTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACg                                                                             <  1:2500530/100‑1 (MQ=255)
                         tacatTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCa                                                                           <  2:301838/100‑1 (MQ=255)
                           catTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAAGAAATGGGTACGCaag                                                                         >  2:4206116/1‑100 (MQ=255)
                            atTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCaagg                                                                        >  1:2436280/1‑100 (MQ=255)
                                    aaTTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTgg                                                                <  2:947620/100‑1 (MQ=255)
                                    aaTTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTgg                                                                <  2:3371832/100‑1 (MQ=255)
                                    aaTTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTgg                                                                <  2:1133242/100‑1 (MQ=255)
                                     aTTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTgg                                                                <  2:4040388/99‑1 (MQ=255)
                                     aTTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGa                                                               <  2:2621874/100‑1 (MQ=255)
                                     aTTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGa                                                               <  1:734819/100‑1 (MQ=255)
                                       tGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACa                                                             <  1:1112375/100‑1 (MQ=255)
                                       tGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACa                                                             <  1:2054107/100‑1 (MQ=255)
                                        gTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACaa                                                            >  2:4194768/1‑100 (MQ=255)
                                        gTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAGGGAAGGCTGGACaa                                                            >  2:1657064/1‑100 (MQ=255)
                                             gTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCa                                                       >  1:443884/1‑100 (MQ=255)
                                             gTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCa                                                       >  2:375541/1‑100 (MQ=255)
                                             gTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCa                                                       >  2:2399985/1‑100 (MQ=255)
                                              ttCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCAc                                                      <  1:3342133/100‑1 (MQ=255)
                                                cAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCa                                                    <  1:4056480/100‑1 (MQ=255)
                                                 aaGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCaa                                                   >  1:1971528/1‑100 (MQ=255)
                                                 aaGGAGAGGATGACACCAGAGGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCaa                                                   <  1:754631/100‑1 (MQ=255)
                                                  aGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCaaa                                                  >  2:729828/1‑100 (MQ=255)
                                                  aGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCaaa                                                  >  1:1852721/1‑100 (MQ=255)
                                                     agagGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAAcc                                               >  2:3911966/1‑100 (MQ=255)
                                                        ggATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTgg                                            <  1:2897958/100‑1 (MQ=255)
                                                        ggATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTgg                                            >  2:3708503/1‑100 (MQ=255)
                                                        ggATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTgg                                            <  1:3746668/100‑1 (MQ=255)
                                                          aTGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCg                                          <  2:1503878/100‑1 (MQ=255)
                                                             acacCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGGAGGCTGGACAACGTCACCGAAACCTGGCGTcc                                       >  2:2789224/1‑100 (MQ=255)
                                                             acacCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTcc                                       <  2:2279397/100‑1 (MQ=255)
                                                             acacCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTc                                        >  1:941453/1‑99 (MQ=255)
                                                              cacCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCa                                      >  2:304727/1‑100 (MQ=255)
                                                               acCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAg                                     >  2:831569/1‑100 (MQ=255)
                                                                  aGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAgg                                    >  1:4219662/1‑98 (MQ=255)
                                                                   gATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTg                                 <  2:1216135/100‑1 (MQ=255)
                                                                   gATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTg                                 <  2:173007/100‑1 (MQ=255)
                                                                    aTGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTgg                                <  2:3792289/100‑1 (MQ=255)
                                                                     tGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGc                               <  1:3928285/100‑1 (MQ=255)
                                                                        aCTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCaaa                            >  1:2001272/1‑100 (MQ=255)
                                                                        aCTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAATGTCACCAAAACCTGGCGTCCAGGGTGGCaaa                            >  1:3582530/1‑100 (MQ=255)
                                                                        aCTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCaaa                            >  2:3645731/1‑100 (MQ=255)
                                                                        aCTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCaaa                            <  1:2256559/100‑1 (MQ=255)
                                                                        aCTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGATACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCaaa                            >  1:662013/1‑100 (MQ=255)
                                                                         cTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAg                           <  2:1532159/100‑1 (MQ=255)
                                                                         cTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGCACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGCCCAGGGTGGc                               >  1:1329258/1‑96 (MQ=255)
                                                                         cTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAGGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAg                           <  1:1657059/100‑1 (MQ=255)
                                                                         cTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAg                           >  2:2128766/1‑100 (MQ=255)
                                                                            gCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGccc                        >  2:2497000/1‑100 (MQ=255)
                                                                             ccAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCg                       <  2:3739358/100‑1 (MQ=255)
                                                                              cAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCGAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGa                      >  1:3470038/1‑100 (MQ=255)
                                                                              cAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGa                      >  2:1750140/1‑100 (MQ=255)
                                                                               aGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGAc                     <  1:2240702/100‑1 (MQ=255)
                                                                               aGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGAc                     <  2:2035188/100‑1 (MQ=255)
                                                                                 aCTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTg                   <  2:2807089/100‑1 (MQ=255)
                                                                                   tGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTgg                  >  1:1706055/1‑99 (MQ=255)
                                                                                    gACGGGCTAGAGCCGCAAGACCAGTAATGAATGGGTAGGCAAGGGAGGCTGGACAACGTGACGAAAAGCTTGCGTGCAGGATGGCAAGGCCCGACTgcat                >  1:491492/1‑97 (MQ=255)
                                                                                     aCCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTc               >  2:3414523/1‑100 (MQ=255)
                                                                                     aCCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTc               <  2:1794862/100‑1 (MQ=255)
                                                                                     aCCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTc               >  2:2918617/1‑100 (MQ=255)
                                                                                      ccGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCa              >  1:2532153/1‑100 (MQ=255)
                                                                                        ggCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACg            >  2:2706412/1‑100 (MQ=255)
                                                                                        ggCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACg            <  2:3337794/100‑1 (MQ=255)
                                                                                          cTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTc          <  2:463797/100‑1 (MQ=255)
                                                                                            ataGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCaa        <  1:3133214/100‑1 (MQ=255)
                                                                                             taGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCAAt       >  1:954671/1‑100 (MQ=255)
                                                                                             taGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCAAt       >  1:550202/1‑100 (MQ=255)
                                                                                             taGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCAAt       >  2:3619705/1‑100 (MQ=255)
                                                                                               gccaccAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCAATGa     >  2:1289316/1‑100 (MQ=255)
                                                                                                  gccAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCAATGAACa  <  1:3869054/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
AAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCAATGAACA  >  NC_000913/2518409‑2518606

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: