Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,282,608 A→G G156G (GGT→GGC yjcF ← pentapeptide repeats protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,282,6080AG98.8% 900.0 / ‑1.6 254G156G (GGT→GGCyjcFpentapeptide repeats protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (131/120);  total (134/120)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 254 total aligned reads displayed.

AATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATACCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGACATTATTTATAAAATAT  >  NC_000913/4282511‑4282704
                                                                                                 |                                                                                                
aaTTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGcc                                                                                                <  2:1290574/100‑1 (MQ=255)
  ttCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCtt                                                                                              <  2:3985087/100‑1 (MQ=255)
  ttCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCAGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCtt                                                                                              <  1:4973819/100‑1 (MQ=255)
      ttttATACATATTACAGTCCTGAAAAATATCGTTTTGTAGTCTGCAACCTGCAAAATTCGCATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAg                                                                                          >  2:308249/1‑100 (MQ=255)
      ttttATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATATATCTTCATCTTTAATCACGATGCCTTTCAg                                                                                          <  1:1804118/100‑1 (MQ=255)
      ttttATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAg                                                                                          <  2:505480/100‑1 (MQ=255)
      ttttATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAAGTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAg                                                                                          <  2:3417449/100‑1 (MQ=255)
       tttATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAgg                                                                                         <  1:593963/100‑1 (MQ=255)
        ttATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGc                                                                                        <  1:298523/100‑1 (MQ=255)
        ttATACATATTACAGTCCTGAAAAATAGCGTTTTGTAATCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGc                                                                                        <  2:4573975/100‑1 (MQ=255)
         tataCATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCt                                                                                       >  1:97430/1‑100 (MQ=255)
          ataCATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATACCTTTCAGGCTa                                                                                      >  1:419360/1‑100 (MQ=255)
           taCATATTACAGTCCTGAAAAATAGCGTTTTGTATTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTaa                                                                                     >  2:3088110/1‑100 (MQ=255)
             cATATTACAGTGCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATTTTCATCTTTAATCACGATGCCTTTCAGGCTAAgg                                                                                   >  2:3585441/1‑100 (MQ=255)
             cATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCCAAATTCGAATTGGGTAAGTCTTCATCTTTAATCACGATGCCTTTCAGGCTAAgg                                                                                   <  2:4760819/100‑1 (MQ=255)
             cATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAgg                                                                                   <  1:673735/100‑1 (MQ=255)
             cATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAgg                                                                                   <  2:1428818/100‑1 (MQ=255)
             cATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAgg                                                                                   <  2:1438413/100‑1 (MQ=255)
              atatTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGgt                                                                                  <  2:1616547/100‑1 (MQ=255)
               tatTACGGTCCTGAAAAATGGCGTTTTGTAGTCTGCAACCTGCAACATTCGTATTGGATAAATCTTGATCTTGAATCACGATGCCTTTCAGGCTAAGgtg                                                                                 >  2:734502/1‑100 (MQ=255)
               tatTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGgtg                                                                                 >  2:1190053/1‑100 (MQ=255)
               tatTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGgtg                                                                                 >  2:1601520/1‑100 (MQ=255)
               tatTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGgtg                                                                                 <  2:3319119/100‑1 (MQ=255)
                atTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGgtgt                                                                                >  2:2655461/1‑100 (MQ=255)
                atTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGgtgt                                                                                <  2:3253602/100‑1 (MQ=255)
                  tACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGa                                                                              >  1:3578442/1‑100 (MQ=255)
                  tACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGa                                                                              <  2:1059780/100‑1 (MQ=255)
                   aCAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGaa                                                                             >  1:100936/1‑100 (MQ=255)
                    cAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGaaa                                                                            >  2:2332350/1‑100 (MQ=255)
                    atctCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGaaa                                                                            <  2:527296/97‑1 (MQ=255)
                       tCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAAt                                                                         <  1:3134299/100‑1 (MQ=255)
                       tCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAAt                                                                         <  1:1570936/100‑1 (MQ=255)
                         cTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCg                                                                       >  1:425112/1‑100 (MQ=255)
                           gAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTa                                                                     >  2:79901/1‑100 (MQ=255)
                           gAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTa                                                                     >  1:2456463/1‑100 (MQ=255)
                           gAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTa                                                                     >  1:2309143/1‑100 (MQ=255)
                             aaaaTAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAac                                                                   >  1:2028679/1‑100 (MQ=255)
                                aTAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAaccac                                                                <  2:2807204/100‑1 (MQ=255)
                                 tAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAaccacc                                                               <  2:664151/100‑1 (MQ=255)
                                 tAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAaccacc                                                               <  1:4765209/100‑1 (MQ=255)
                                  aGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAAccacca                                                              <  1:1211458/100‑1 (MQ=255)
                                    cGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCAtt                                                            <  2:3710773/100‑1 (MQ=255)
                                     gggtcgTAGGCTGCAAGCGGCAAAATGCGAATTGGATAAAGCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCAttt                                                           <  1:5152388/95‑1 (MQ=255)
                                     gTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCAttt                                                           <  2:3203631/100‑1 (MQ=255)
                                      ttttGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTa                                                          <  1:259975/100‑1 (MQ=255)
                                        ttGTAGTCTGCCACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTCTGAAAAATCGTAACCACCATTTAct                                                        <  2:4056663/100‑1 (MQ=255)
                                         tGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTActg                                                       <  2:1340936/100‑1 (MQ=255)
                                         tGTAGTCTGCAACCTGCAAAACTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTActg                                                       <  1:142053/100‑1 (MQ=255)
                                          gTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTActgt                                                      >  1:3794306/1‑100 (MQ=255)
                                          gTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTActgt                                                      <  1:2962495/100‑1 (MQ=255)
                                          gTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTActgt                                                      >  2:642768/1‑100 (MQ=255)
                                          gTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTActgt                                                      >  2:2946341/1‑100 (MQ=255)
                                           tAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTActgtc                                                     <  1:3427574/100‑1 (MQ=255)
                                           tAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTActgtc                                                     >  2:3500410/1‑100 (MQ=255)
                                            aGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTActgtct                                                    <  2:2975056/100‑1 (MQ=255)
                                             gTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTActgtctg                                                   >  1:4410909/1‑100 (MQ=255)
                                             gTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTActgtctg                                                   >  1:1969462/1‑100 (MQ=255)
                                              tCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTAAGGCTAAGGTGTGAAAAATCGTAACCACCATTTActgtctgt                                                  >  1:4506616/1‑100 (MQ=255)
                                               cTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCtgtg                                                 >  1:1287981/1‑100 (MQ=255)
                                               cTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCtgtg                                                 >  1:4315504/1‑100 (MQ=255)
                                                tGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTgg                                                >  1:17759/1‑100 (MQ=255)
                                                tGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTgg                                                >  1:2140725/1‑100 (MQ=255)
                                                 gCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTggg                                               >  1:2228664/1‑100 (MQ=255)
                                                      cTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTc                                          >  2:1359026/1‑100 (MQ=255)
                                                      cTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTc                                          <  1:3395526/100‑1 (MQ=255)
                                                         cAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGAt                                       >  1:321415/1‑100 (MQ=255)
                                                         cAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGAt                                       >  1:78556/1‑100 (MQ=255)
                                                             aTTCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGta                                   >  2:3917679/1‑100 (MQ=255)
                                                              ttCGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGtat                                  >  1:151119/1‑100 (MQ=255)
                                                                cGAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGtatta                                >  1:4072919/1‑100 (MQ=255)
                                                                 gAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGtattat                               <  2:4182976/100‑1 (MQ=255)
                                                                 gAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGtattat                               <  1:1043493/100‑1 (MQ=255)
                                                                 gAATTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGtattat                               <  1:3317255/100‑1 (MQ=255)
                                                                   aTTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACGACCATTTACTGTCTGTGGGATTTCGATAGtattatta                             >  2:4616138/1‑100 (MQ=255)
                                                                   aTTGGATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGtattatta                             <  1:517111/100‑1 (MQ=255)
                                                                      ggATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATgg                          <  2:2479601/100‑1 (MQ=255)
                                                                      ggATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATgg                          >  2:1056858/1‑100 (MQ=255)
                                                                      ggATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATCATgg                          <  2:4320340/100‑1 (MQ=255)
                                                                       gATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGa                         >  2:226169/1‑100 (MQ=255)
                                                                       gATAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGa                         >  1:89612/1‑100 (MQ=255)
                                                                        aTAAATCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGAc                        >  2:2400650/1‑100 (MQ=255)
                                                                          aaaTCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACgg                      <  1:4876662/100‑1 (MQ=255)
                                                                          aaaTCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACgg                      <  2:4477518/100‑1 (MQ=255)
                                                                           aaTCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGt                     >  2:1661916/1‑100 (MQ=255)
                                                                             tCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGa                   >  1:3730142/1‑100 (MQ=255)
                                                                             tCTTCATCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGa                   >  1:2812794/1‑100 (MQ=255)
                                                                                  aTCTTTAGTCACGATGCTTTTCAGGCTAAGGTCTGAAAAATCGTTACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGACatta              <  1:4733816/100‑1 (MQ=255)
                                                                                   tCTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGACagatc             >  2:4062346/1‑96 (MQ=255)
                                                                                   tCTTTAATCACGATACCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGACattat             >  1:2745108/1‑100 (MQ=255)
                                                                                    cTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGACattatt            >  2:752/1‑100 (MQ=255)
                                                                                    cTTTAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTCACattatt            <  1:632056/100‑1 (MQ=255)
                                                                                     tttAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGACAttattt           <  1:4600442/100‑1 (MQ=255)
                                                                                     tttAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGACAttattt           >  2:1764595/1‑100 (MQ=255)
                                                                                       tAATCACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGGGGGATTTCGATAGTATTATTATGGACGGTGACAttatttat         >  2:1016299/1‑100 (MQ=255)
                                                                                         aTCACGATGCCTTTCAGGCTAAGGTGTGAAAAATAGTAACCACCATTTACTGTATGTGGGATTTCGATAGTATTATTATGGACGGTGACATTATTTATaa       >  2:4691861/1‑100 (MQ=255)
                                                                                           cACGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGACATTATTTATaaaa     >  1:2655986/1‑100 (MQ=255)
                                                                                            aCGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGACATTATTTATAAAat    >  1:4420693/1‑100 (MQ=255)
                                                                                            aCGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGACATTATTTATAAAat    >  1:1385212/1‑100 (MQ=255)
                                                                                            aCGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGACATTATTTATAAAat    >  1:2351942/1‑100 (MQ=255)
                                                                                             cGATGCCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGACATTATTTATAAAata   <  2:1944720/100‑1 (MQ=255)
                                                                                                 |                                                                                                
AATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTGGATAAATCTTCATCTTTAATCACGATACCTTTCAGGCTAAGGTGTGAAAAATCGTAACCACCATTTACTGTCTGTGGGATTTCGATAGTATTATTATGGACGGTGACATTATTTATAAAATAT  >  NC_000913/4282511‑4282704

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: